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^ FAULT-TOLERAi^f TIME-OUT COMMUNICATION PROTOCOL 
AND SENSOR APPARATUS FOR USING SAME 

CROSS-REFERENCE TO RELATED APPLICATION 
The present application claims the benefit of provisional patent application No. 
60/450,000, filed February 25, 2003, which application is hereby incorporated by reference. 

TECHNICAL FIELD 
The present invention relates to sensors, e.g., sensor arrays and networks. More 
particularly, the present invention pertains to communication methods for use with such sensors. 

SUMMARY OF THE INVENTION 
The present invention provides various embodiments of communication protocols for 
use, for example, in communication of sensor information in system, sensor network, sensor 
array, etc. 

In one embodiment of a method of communicating sensor information, the method 
includes providing a sensor architecture comprising a base station and a plurality of sensors 
associated with the base station. Each sensor is operable to communicate sensor information to 
the base station and backup information to one or more other sensors of the plurality of sensors. 
Failed communication of sensor information from at least one sensor of the plurality of sensors 
to the base station is detected, e.g., a time-out occurs. Task information representative of the 
failed communication is broadcast to the one or more other sensors receiving backup sensor 
information from the at least one sensor to request priority information therefrom. Backup 
sensor information is requested to be sent to the base station from one of the one or more other 
sensors receiving such backup sensor information based on the priority information received 
from each of the one or more other sensors. 

In another embodiment of a method of communicating sensor information, the 
method includes providing a sensor architecture comprising a base station and a plurality of 
sensors associated with the base station. Each sensor is operable to communicate sensor 
information to the base station and backup information to one or more other sensors of the 
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plurality of sensors. Failed commumcation of sensor inf<Hiiia!ion from at least one sensor of 
tt^ plurality of sensors to the base station is detected. Upon <tetecti<n)» an evaluation is 
p^onned of the potential ability of each of the one or mcxe oQxer sensors in providing 
backup s^isor infonnation to the base station. The backup sensor information is requested to 
be sent to the base station from one of the one or more other sensors based on the evaluation. 
For example^ the evaluation maybe p^ormed by considmng tihe time needed for each sensor 
to transmit the backup sensor information to the base station and/or the number of tasks 
offered to each sensor. 

Hie present invention contemplates the use of the sensor communication described 
herein in various sensor networics, apparatus, or systems. Hardware, programs and/cnr 
algorithms fc^ use in providing the described above features accordmg to the piesent 
invention are also contemplated. 

Hie above smnmaiy of the present invention is nc^ int»ded to describe each 
embodiment or every implementation of the present invention. Advantages^ together with a 
more complete und^tanding of the invention, win become s^aient and appreciated by 
referring to the following detailed description and claims. Various otter embodiments of the 
present invention are provided by tte telow Description. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 

In the following detailed description of the embodiments, reference is made to 
^ific embodiments in which the invention may be practiced. It is to be understood that 
other embodiments may be utilized and {processing step/structural changes may te macte 
without departing from tte scope of tte present invention. Further, the word •*backup," wten 
used, for example, in association with a sensor or node, may te interchangeable with tte tenn 
**sibling.** 

Tte following TV. Detailed Description of tte Invention*' provides various 
embodiments of the present invention. It will te recognized that various components or 
elements of tte present invention and/or method steps carried out as descrited therein may te 
included and/or may te optional according to the present invention. 



IV. Detailed Description of the loTention 

Abstract- 
Flow measurenant is an in^rtant technology used in nearly every industry. Distributed micro flow-sensor 
arrays and networics (DMFSA/N), buflt from coUcctions of spatiaDy scattered, cooperating intelligent 
micro flow-sensor nodes, can improve the accuracy and reliabiUty of the system. However, it is unrealistic 
to expect ail the sensor nodes and communication links in the system to be fault-free all the time. A fauh- 
Udeiant time^mt Piotocol (FTTP) with two alternatives . one widi blhid station, and die other with smart 
base station was proposed in our research. Experimem results diow that compared widi FTTP with bUnd 
base station* although FTTP widi smart base station requires mm complex structures, it invohres less 
communicatim. 

L Distribnted Micro Flow<Sensor Arrays and Networls (DMFSA/N) Structure 
Before FTTP is presented in detail, it is necessary to briefly describe the DMFSA/N structure that suppwts 
tfje communication among die nodes in die system. Rgure 1 shows the clustered DMFSA/N structure. In 
diis structure, DMFSA/N are divided into sensor cluster tmiis (SCUs). each of which consists of a set of 
tntdUgent micro flow-sensor nodes and a base staticm. The inieUigem sensor nodes widiin the same SCU 
are also called sibling nodes. An intdligem micro flow-senscM* node consists of a dedicated processing unit 
and an associated nucio flow-sensor diat measures die variable of interest A more powerful procesring 
unit is used as die base station, which acts as die control node of die respective SCU, controlling die signal 
transmission of die sensor node, mana^g information rerouting m case of link fatlores and integrating die 
information from aU die sensor nodes in die SCU. In order to make die system tolc^ 
sensor node sends information not only to die base station, bm abo to one or more odier sibling nodes diat 
arc also called its backup nodes. Two communication jffotocols can be used among die sibUng nodes: (1) 
Broadcasting Protocol, in which each senses- node broadcasts its information to aU its sibling nodes; (2) 
Bvenfy wide gossiping protocol (EWGP), a revised Gossiping Ph)tocol (GP), in which instead of sendmg 
information to only one odier sibling node by random selecticm according to GP» die baclcup nodes are 
diosen in die way diat die n>]es backup nodes are evenfy distribmed a^ 

not hai^ that some nodes get too many 'duties*. whUe odiers are 'starved*. These communicatioii 
protocols can also be applied in die communication among die base sution. The external processor b die 
commander of die entire nctworic, in which infwraation fiom all die SCUs is integrated for decision makiqg 
as required. The clusters can be dqiloyed aloi^ die fltow padi, dqpending on die flow situation and 
qsplications (Figure 1). 
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IL Fault-tolerant TiiDe-<mt Conummkation Protocol 

The design of &alt-tc^eraQt timeout protocd (FTTP) is motivated by the comlunation of time-out and ta^ 
coofdtnalioD protocols and the need for fattlt*tolerant integration. At the beginning of each data 
transmission cycle^ the base station broadcasts the data trananission request to all the sensor nodes in the 
SC!U. Recdving the reqi^ each sensor node sends its data to the buffer inside the base station. Each 
sensor node has a corresponding buffer in the base station to stcwre its data so that the base station can trwc 
the source of the data. Then the data will be retrieved into the processing unit In the base station. According 
to FTTP, the base station will stop waiting for the information irom the sensor node if a obtain amount of 
time. 7» has passed, where r is the pre-determined tBne<)m period 

Base station will then announce the rerouting task message to its siblmg nodes. Two alternatives exist for 
vteie the base statim should send the task annrancemem: (I) with a Blind Base ^ftadoiv Ae message is 
sort to aU the sii^ nodes; (2) with a Smart Base Statioo, the inessage is o^ 
can finish the task carHcst, based on the cuirart status <^ the notom the system 
FTTP with Blind Base Station 

In this protocol, when the link of a sensor node, say node v^, fails, the base station will broadcast the 
reroutmg task message (RTM)f<» v^^ to aO its ailing nodes. IUx:dving RTM» its backup iiode £^ 
backup node of v^) generates a bid value frvii and sends it to the base station. 6ai€ ^A^t, where JSJIT^ is the 
set of backup nodes of . M^i, can be calculated as follows: 



«*ere «, Is Ae transmission Tale of node 6« (time ^ ta tt^ 

of tasks that have been offered to node Becanse it is unrealistic for the intelligent sensor node to keep 
track of its transmission rate over time due to its computation limit, », is assumed to be constant for node 
6«. The base station then evahiates the priorities of the backup nodes based on dieir bid values (die lower 
a» hid value, the higher dte priority), and sends the task offer message (TOM) to the backup node widi the 
Wjhesi piiorily. b^,^ Receiving the task offer message. sends dns required backup data message 
(BM4) to the base statloa The stnictiire$ of the base station and inlelHg^ 
FITP witti blind base statioB are shown in Rgure 2 and 3 respectively. 

FITP M>ith Smart Base Statbm 

In this protocol, in order to limit die unnecessary communication and reduce die communication traffic, 
instead of broadcastmg RTM for v». whose link fidls. to aD its sibling nodes, the base station first 

evahates die priorities of die its backup nodes and then sends lOM only to the node tfiat can fi^ 
dieearUest The time when die backup mide finishes die task. can be cafcadated as (assummg die 
transmission Starts at time 0): 

tu=imu+l)ru (2) 
where m^ls die number of tasks diat have been offered to die bachqi node and y,, is die 
transmissten rate of *!, (time spert ta transnritting one message). The base statiOT 
node whose is die smallest In order to find out die base station needs to have a kro 

base which keeps track of die Inforaiation for each node - its backup nodes, transmission ra^ 
number of tasks diat have been offered Kgure 4 shows die structures of base st^ 

nodes for die FTTP widi smart base station. TlK modd of die knowledge base is shown in Rgnre 5. i^ 
which « is die munber of sensor nodes m die chister and die backup nodes are hypodieti^ 
only. Figure 6 shows die kigic chart of die FTTP widi smart base station. 

Compared widi die FTTP widi blind base station, die FTTP widi smart base station requires less 
communication, and has die advantage of being able to keep track of die updated transmission rate of each 
sensor node, but it requires a base station widi more complex structure (widi an additional knowledge 
base). Therefore, a tradeoff has to be ccmsidcred m order to determine when to apply each protocoL 




Fi9ire2. Structures ofthe BasdSfiatfon and Ihtdligeat Sensor Nodes for the FTTP wiA Blind Base 
Station (Dasb line means the link may or may not e^dst dq)^^ 
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Figure 3, Lo^c CbM ofthe FTTP with Blind Base Station 
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Fiffm 4. Stm^ves of the B ase Nation aood hUteHignit Sensor Nodes for ^ FTTP with Sman 
Base Station (Dash ttm means the link may or may not exist depending on the protocol used) 
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Rgiire 5. Logic Oiait of the ITn> widi Smart Base Station 
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The following ^'Vlt Infonnation** provides infonnation legardiBg one or more 
applications for the present invention and problems that the present invention may be used to 
resolve* 



: Infonii^kMB 



Distributed micro flow-seroor arrays and oetwoilcs (DMFSA/N) will open an origiiial directioii In tiie 
af^licattons of flow measurement and le^ detection. The atnlity to place tl^ micro flow-sensor nodes 
when and where diey are needed provides better control of process. However, It is unrealistic to expect all 
die mkro sensor nodes and communicatioa links along die path to be &tilt*firee all die time because of the 
constraints of mao sensors and com|^ flow enviromnent In order to make die system link fiallure 
tolerant and satisfy the teal-time control requirement, fault-tolerant time*<nit communicatii»i protocol 
(FTTP), which has two alternatives, is proposed • FTTP integrates theccmventional time-out 

protocol with ftutt-tc^eraot inlbnnadon tnt^ration characteristics* 



i Lit addresses die feolt-toleranceissiK in die design of DMFSA/DMFSN 

2. This i»otocol is easy to understand and implement 

3. It can easily be revised to cater to the applicadons similar to the DMi^ 

It addresses the issues in the ctesign of distributed msao flow-sensor anays and networic 
(DMFSA/DMFSN) or odier similar appUcations. Due to die physical ccmstraints of die micro sensors and 
effects of die environment diey are measuring, it is unrealistic to expect all die micro sensor nodes and 
communication links along the path to be fault-free all the time. In ord& to make the system link failure 
tolerant and satisfy the real-time control requiremem, fault-tolerant time-out proiood* which has two 
alternatives^ is proposed in our research. 

In the past, similar problems were raised in the design of sensor n^woiks. The most common apfanoach was 
to increase the number of the connections among die sonsor nodes or the number of sensor nodes. 




The following materials entitled 'l>istriboted Micro Flow-Sensor Network (DMFSN) 
Design and Modeling^ provide generalized informaticHi about a potential communication 
protocol. 
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The following article entitled 'THstributed Micro Flow-Sensor Network (DMFSN) 
Design and Modeling'* provides generalized information about a potential communication 
protocol. 



Distributed Micro FloW'-Seiisor Netwaric/Array Design and Modeling 



Abstract This ^apex describes the design issues of die distributed micro flow-sensor networks/arrays 
O^MFSNypMFSA) and sensor netwcnfc aichitecturcs and infonnadon conmmnicatioQ jHotocols diat can 
be used in DMRSN/DMFSA. A sensor networlc/afray Emulator and evaluator, c^led Teamworic 
Integration Evaluator <TIE)/MEMS, was developed to evaluaie die eifkiency of DMFSN in terms of its 
comnnmication cost Initial simulation results diovs duu die duster sessor network provides better 
eCRciency than-hierafchicid and committee sensor netvorlcs* 

Keywords: MEMS, Networik Aicbitectnres, Plroiocols* Fault-Tolerant. 11nie-0ut» Teamw(»k Ihtegnuion 
Evaluator 



L INTRODUCTION 

The measurement of fluid flows is an important technology used in nearly every industry. Flow sensors are 
used to measure dte movement of the fluid flow (liquid flow or gas flow). They are widely used in a variety 
of ai^lications, such as automotive applications, envimnmental monitoring, biomedical applications, etc. 

Before die advent of micro omiiimtzation, senses systems tend to be bulky aiui es^iensive, so tsag]& sensor 
systems were used to save space and cost The advant^ of single sensor systems is that diey are relatively 
easy to construct and analyze. However, single sensor systems have limited applications and usages. Fdr 
instance, if we want to measure several variables, say, tenperature, pressure and flow density, at d» same 
time, single sensor systems cannot be used. In addition, a single soKsor cannot guarantee to deliver 
accurate inform^n att die tioie, because It is me^taUy cfaaUenged by iioi^ 

Microelectromechankal system (MEMS) technology deals widi designing and falmcating entire electrical 
and mechanical systems, usually on a single silicon dap. Micro fbw-sensors are one of the most popular 
MEMS devices. Due to dieir small size, the advanta^ of micro flow-sensors are: 

• Interfere less wxdi die enviroimient they are measuring 

• Lower manufacturing cost (less materials) 

• Can be uppUed in narrow spaces, sudi as innde die living organisms, small {npes, automobile 
engines, etc. 

• Multiple sensors can be used widi redundanqr to improve infonnation accuracy and fault tolerance 
of die system 

Distributed nncro fk>w-sra$or networts/arrays (DMFSN/DMFSA) are built from collections of qiatially 
scattered, intelligent micro flow sensor nodes. Each node has die ability to measure die local flow widiin its 
accuracy limits, processes die raw sensor data, and cooperates with its neigliboring nodes, which can 
improve die accuracy and reliability of the information. Sensors incorporated with dedicated signal 
processing functions are called intelligent sensors, or smart sensors. The main roies of dedicated signal 
processing are to eiduaK:e design flexibility and realize new sensing functions. Additional roles are to 
reduce loads on central processing units and signal transmission lines by distributing information 
processing tt> die lower layers of die system (Yamasati, 1995). 
Developing effective DMFSN, however, also faces great challenges, some of which are: 



« Scaling issues. The differences in dimensions between the micro woHd and macio wmid cause 
differences in balance between magnitudes of forces. 

• Suitable micramachining technique. To achieve the realization of nucix) flow sensor stmcture^ the 
stnicnue has to be released using nucromacfainrag techniques CToslw and Fumihito. 1999X 

• Appropriate network architecture. Since micro sensor nodes ait distributed spatially, aj^priate 
Sensor network architecture has to be found to support efficient, undisturbed commimicadon 
among die sensor nodes. 

• Optimal placement and deployment of micro flaw sensors. Flow meastirein»t is a con^lex 
process because of the numerous variables diat can affect accuracy, such as flow stabtU^* nm* 
homogeneous flow, temperature variation along die flow path* etc. Hence, die placement and 
depk^rment of mkro sensors m die flow space (e.g., posttkining widiin the cross*section of a pipe 
and along the axes of flow) is hnpcmant inorto to nitnimize the enx»r due to uncertain foctors. 

• Fault-tolerant informa^on commmication protoa^ and integra^an algorithm. Flow measurenaent 
is a complex process, so it is imrealistic to expect all die sensor nodes and communxcatkm 
along the path to be fault-fiee alt die dme. 

• Low eommurdcation cost. Co mmnnicatfo a cost is often related ID communication time and power 
con$uii^<m. Minimizing communicatton cost is in^oitant in online control applications and 
energy-constrained appficaiions. 

The objectiveofdiis pi^ is to discuss die sensor network architectures and communication protocols diat 
can be used for DMFSN/DMFSA, and die development of a sensor network simulation tool, TIE 
CTeamwi^ Integration EvaluatorVMEMS. The purpose of TGBAIEMS is to simulate and evaluate die 
efUctency of die architecmres and protooob. In diis paper, sone recent proposed mkro flow seas<»' 
designs are presented in section 11; Different ^pes €i archhectures and communicaticxi protocols diat can 
be used for DMFSN are discussed in section m and IV, respectively; Section V describes TIE/MENfS 
structure and preliminaiy simulation results. Finaify, conclusions and further research arc discwssed xq 
sectbnVI. 



n. MICRO FLOW-SENSOR DESIGN REVIEW 

Realization of micro flow-sensors can use several, different opmting princi{^ which inchide 
aneoraieter functions, physical deflections, tempeiauire^qimieiit layers, time^-fli^ coriolb^fbfce, 
etc 

Rgure 1 illustrates an expkided view of a silioon diaphragm capacitive sensor, called Piandd micro flow 
sensor (PMFSX for flow-velocity measurement (tesigned by Boberig et oL (1998). Like die well-known 
classical prandd tube, the sensor realized fiow-velodty detection by irvasuring the piessme diffetcace 
betwera die stagnant fluid pressure in firont of die sensor and static pstssmt in die flow around die sensor. 
This difference results in a deflection of a siliom diaphragm suspended boss, which servers as die comner 
dectrode of an mtegtated cq>acittff diat is direcdy exposed to die fluid 10 be measured 
figure 2 shows an exaiKq>le of the mass fk>w sensor designed by Nguyen and Kiehnsch^ (1995). This 
diermal flow sensor contains polysilicon heaters and measuremBnt resistance on a dii^^gm. The fluid 
flows in an anisotropically etched micro-channel widi a cross*section of O.Gvaa? and a lengdi of IQnm 
The diannel is covered widi a Pyrex glass plate by anodic bonding. A small thermal response time fiom 1 
to 4ms is reached with a di^thragm ttnckn^ from 10 to 30ura 




1. Pofysilicon resistance 

2. Aluminum conductor 

3. Silicon chip 

4. Pyrex plate 

5. Flow input and output 



Figure 1. Expkxied view of die Prandd micro fiow-sensor developed by Berberig et aL (1998) 
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m DiSTRuinrED nocro sensor network ARCorrEcruRE 

A distributed micio sensor network is a collection of nncro intelligent sensor nodes that are cK^nbuted 
s{>atxafiy. Its d)straction is a Gr^ G=(V. E)» whete V re^tsents^tfte inteHigeat micro sensor iKxIe. This 
node consists of a micropfocessor unit and an associated mao sensor. £ represents the comnamicaiioii 
links connecting the sensor nodes. 

The two most common^ conventional architectures for distributed s^oscar networks are die Committee (or 
Anardm) and the Hlemrdiicai architecnires, first proposed by Wesson et oL (1981). In die 0»nmittee 
stmcture. each node in die network is autonomous and can send information to all other nodes m the 
netwwk. Thus if there are N sensors in the network, it needs interconnections that are of die onfer N^ The 
committee architecture is not suitable for a large distributed micro sensor network, because it is too 
expensive with respect to the communicaticm cost and not easily extensible (Figure 3(a)). In the 
Hierarchical structure, sensor nodes are connected to form strict hierarchies. Each node (except die lowest- 
level leaf node) receives information fiora the lower level node; integrates the infiMrmadon received 
according to its posidon in the hierarchy, and tten sends infonDadon tp die nock at tiie upper leveL The 
root uode can be a more powerful jEnocessor, ct a nonM sensor node la wlddi di^ final icsuh is generated. 
The H^rarchical aichitechwe has s^end advantages, sudi as constant node degree and easy extensibility^ 
iHit it may not produce accurate inf<vmation at hi£^ leveb because the errors cm aocamnlate ^n/hea the 
tnform^n goes up the hierarchy. It is also not tolerant to link ^ures (Figure 3(b)). 
Quster structure overcomes the majcH* limitations of both the Committee and die Hierarchical architectuies 
by limiting the number of oMiununicadon channels. At the same time, the cluster architcctute tolerates die 
sensor node and link failure. £u the chist^r an^tecture, die s^asor network is divkletf into Mit^ 
wtits(SCUh each of which oxisists of a set of intelligent sensor nodes. An im^igent seasor nodeooiBisis 
of a processhig unit and an assodated nscio sensor dial measmts the variables of mleicst A sinknode & 
needed fctr eadi cluster to mtegrate information from all die sensor nodes in the duster. Considering die 
complexity <^ co m pt itati op and enogy constraint of the mcso sensor node, a ram powerful p rocessor is 
chosen as the base station for each cluster. Figure 3(c) shows the chister architecture for die distributed 
micio sensor network, which is die revision of die network tq)ology proposed by lyengaor er dL (1994). hk 
dus stmcture, the base station acts as d^ control node die respective SCU, msms^^g mfkumadon 
rerouting in case of link failures and combming the mformadcm from all the intelligeiit sensor nodes in the 
SCU. The hitelligent sensor nodes widtm one SCU (also called sibUng nodes) can be fhliy connected (using 
broadcasting coimmmicatlon protocol) or connected with only one or several selected sensOT iKxIe (using 
gossiping comnnmtcation protocol), d^nding on the link failure probability of the system. Broadcasting 
and gossiping communicadon protocob are described in more detail in part IV. Base stations in die 
network can also be fully or partially interconnected, for the same reason as connections among sensor 
nodes. The external processor is die commander of die distributed micro sensor network, in which 
information from all the SCUs is integrated to produce a final result, or make decisions as requiredL The 
clusters can be deployed along the flow padi, depending on the flow situation and applicadons. 
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IV. INFORMATION COMMUNICATKW PROTOCOLS 



1. PoUa-tQ-Point comnumicaihn pnaocal (FTPh Id this piotoool, seom nodes A and B 
coummnicate exclusively with each other widwiil interfering widi otli«^ nodes. (Ezhildidvan, 
1989» Henozebnan ei al., 1999) 

2. Broadcasting c<mtmunication prctocoL En dus protocol, the tnformaticm ont by a smsor node 
is leoeivcd by aU die sensor nodes widiin a ceitatn lange of the send^ 

piolocol is dial vidien a ^gle node sends kfcMuutf^ 

nei^b(»ing nodes using only one transmission (Hedetmoni. 1988, Pete, 1996). 

3. Gosling communication protocol In diis protocol, instead of indiscriimnately sending 
infonnation to aO its nei^UxMring nodes, each s^isor iKxte only foiwaids d» data to (me of its 
neighboring nodes by randomly selecting or accc»ding to sc«ne criteria, such as die distance 
between those nodes (Hedetnieim, 1988» Pieic, 1996). 

4* Pauit-toterant time-out protocol (FTTP). This protocol is an In^vemeot of die tradidonal time- 
out protocol widi fault-tolerant propedies. Aconding to the time-out fvotocol, the receive node 
will stqp waiting for die infonnadon finom die sender node if a certain amount of time, T, has 
passed, if the committee ardutecture or cluster architecture is used for the sensor netwoik, the 
receiver node can announce die reroming task to the sensor nodes, which have the bacfci^ 
infonnation of die sensor node widi die fialed liidL &i order to reduce the message trafRc and limit 
the unnecessaiy communicadott, instead of broadcasting the task announcement message to all 
odier nodes, the receiver node can only sokI the message to ^se most suitable nodes based on 
their previous performance. The receiver nodes have to have 'iGDowIedge** about die 
communication abilify of those nodes in order to make such 'Intdligenf* <fecision. Theiefbie» n 
database that stores sudi knowledge Meds to be buflt in the receiver node. In the chister sensn* 
network, the base stations and die external processor are examples of such receiver nodes. For the 
committee sensor networic, however, it b not economic for each sensor node to have a knowledge 
base in its dedicated processing unit Figure 4 (a) stows die foult-tolerant time-out jnotocol for the 
cluster sensor networic, in which St, Sj.. . .,Sb are die micro sensor nodes, n is die nun^r of sensor 
nodes in die chister; Bi, B2, are die buffers for sensor nodes St, S3, Sy,,,Sn respectively. 

Figure 4(b) shows the knowledge base structure, which consists of three attributes, the sensor node 
ID, backi^ nodes which have the redundant informadon of die sensor node, and die transmission 
delay time for the sensor node. Transmission delay dme is the measure of ^e performance of the 
sensor node, which is equal to die difference between die time when die base station receives die 
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ittfonnanon from the sensor node and the time when the base station issues the task announcement 
to the sensor node. ^ 
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Figure 4. (a) Fault^tolmnt time-out protocol (FTTP) nx>del 
(b) Knowledge base structure 



V. TIECrEAMWORKINTERGRATION£VALUATOR)/A^^ 

The objective of TIE/MEMS, an extended version of die TIE CTeamwork Lite^ation Evaluatkm) SmaaSmt 
devek^ied by PRISM (Khanna and Nof» 1994, Huang and Nof, 1996, Anussomnitisam and Nof» 2000), is 
to slnnilate and evaluate the architecture and information integration and comnmnication protocol 
developed for die distributed micro sensor netwoik. TEE/MEMS is written using MPI (Message Passing 
IntdfaceXSnir et at, 1996). MPI can be used to simulate die dsfiferem conmumicatimi schemes among die 
sensor nodes, such as point-to-point c<»iimuiiication» collective cornmimicatkm, group anmnunicatioo, and 
so on* In addiiioii, implementations <^ MPI on top of standard Unix inter-piooessor communication provide 
portalnlity to multiple processors, woikstation chisters. and heterogeneous n^oiks of workstations 
(NOW). 



5./. Stmcture ofTW/MEMS 

The overall structure of TIE/MEMS is depicted in Figure 5. Distributed micro sensor network architecture, 
information communicadon protocol, sensor network parameters, sudi as the number of intelligent micio 
sensor nodes, and constraints unique to MEMS sensors (e.g., physical and chemical constraints can»ng 
micro-stress, proximity noise, and other disturbances to measurement and communication) are the inputs of 
TIE/MEMS; the outputs are the communication costs, which are represented by die simulation time and 
energy consumption. If the communication costs are too high, we need to revise the network architecture 
and/or die information communication protocol. Table 1 lists some examples of constraints related to 
MEMS sensor operations. 
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Figwre 5. Structure of TIE^MEMS 



Table I. Examples of constraints related to MEMS seosois 



Scaling issues 


Viscosity force; Surface effect 


Ptoblons in nucnmiachining 


Squeeze film damping; F^cular contanunation; Static stiction 


Response 


E>eciease of lespcmse time diie to pn>tective ^Is» coatings or 
diaphragms. 


Reliability of oiicioelectionic device 


Breakdown of a wire-bond to a junction; Defective encapsulation; 
Defective setmcondnctor suiterial; Theimal nmaway 


Ufctkot 


Aging; Contamination of die device Titr^Hf^f pow^ supply 


Noise^tofeiences 


Degraded signal due to die variance of temperature, humidity, €tc 



5*2. PreUmimuy simulation result 

Plreliminary simulations were nin on Origin 2000 workstations at Con^uter Sciaixe in Purdue Universt^. 
Table 2 shows the results of normalized communicadon time fiom 90 simulation runs for the sensor 
netwofk architectures and communicati<M protocob described in section IV. The communication time of 
die hierarchical network is used as the baseline. Table 3 shows die ANOVA results (oM).95) of die data in 
Table 2, From Table 3, it can be seen diat when die nuinber d sensor nodes is smaU, the communication 
cost of all network architectures and communication protocols is similar. Therefore, committee architecture 
or cluster architecture might be preferred in order to in]{^ve the network fault-tolerance. When the nunto 
of sensor nodes increases, the performance of committee network, especially the committee netwoiic with 
point-to-point communication protocol deteriorates dramatically, so die hierarchical netwoik and cluster 
network may be preferred. In eidier case, cluster network is a recommended choice. 

Two important points need to be mentioned: 

1) The communicadon cost of the hierarchical sensor networic is greatly influenced by the height and fan- 
out of die tree structure. The less height and larger fan-out the structure has, the smaller die communkadon 
cost is, but die more vulnerable it is to die link and/or node faOures. 



2) The communication cost of the cluster sensor network is greatly influenced by the number of sensor 
nodes in each cluster and the number of clusters. The less the number of sensor nodes in each cluster, and 
Che more clusters the sensor network has, the less communication time needed within each cluster, but morn 
communication time may be needed among the dusters, so we need to calculate the optimal nunniber of 
clusters and number of sensor nodes in each cluster to minimize the communication cost 



Tab!e2. TIE/MEMS analysis results: Normalized communicadon tin» vs. No. of sensor nodes 
(as ratio of the communicadon time in the Iffierarchical architecture) 





Sensor network alternatives 


No. of sensor nodes 


» 


Architecture 


Communication 
Pkotocot 


4sensor 
nodes 


8 sensor 
nodes 


15 sensor 
nodes 


21 sensor 
nodes 


30 sensor 
nodes 


1 


Committee 


FTP 


3.1783 


1.1027 


3^083 


6.7054 




2 


Comminee 


Broadcasting 


13852 


1^857 


1.8919 


1.8335 


WA* 


3 


Hiei^ichical 


PTP 


1 


1 


1 


I 


1 


4 


Cluster 


Broadcasting 


U274 


0.2968 


0.1145 


0.Q5Q2 


0.0159 


5 


Ouster 


Gossiping 


1.1651 


0.2811 


0.1027 


0.0419 


0.0145 



*Data not available because the number of messages exchanged exceeded the maximum 
limit stq^rted by Origin 2000 



Table3. AVONA results (Ges95%) of communicadon time (cost) vs. No of sensor nodes 
(1,23.4,5 correspond tolf 'm Table 2 respecdvely; i>i means die cost of i is larger than 
that of j) 



No. of sensor 
nodes 


ANOVA results 


4 sensor 
nodes 


no significanUy different among all network architectures and communication 
protocols 


8 sensor 
nodes 


1» 2 and 3 are not significantly diiferent, 4 and 5 are not significandy 
different* but 1,23 are significanUy different from 4 and 5 and 1,23>4>S 


IS sensor 
nodes 


4JS are not significandy different, but 4JS are significandy differaat from 1,23 
and 1,23 are also sign^kandy different 1>2>3>43 


21 sensor 
nodes 


44 are not significanUy different, but 4,5 are significandy difteient from 1,23 
and 1,23 are also signiHcandy different 1>2>3>43 


30 sensor 
nodes 


4,5 are not significandy different, but 4,5 are significantfy different from 3 
and3>43 



VL CONCLUSKmS AND FUTURE RESEARCH 

Alternative sensor network architectures and communication (m>tocob are discussed and comi>ared in the 
paper. TIE/MEMS is a modeling tool developed to simulate and evaluate the cost effectiveness of the 
network architectures and communicadon protocob. 

Several other issues, such as the optimal deployment of micro flow sensors, and fault tolerant algorithms, 
are also considered in diis research. These issues depend on die flow environment, sensor physical 
constraints like data accuracy, time'<tependent behavicMr, and energy constraints (for wireless 
communication}. Therefore, in the next step, dM>se constraints will be investigated to find out die number of 
sensor nodes needed to achieve the required data accuracy. TDE/MEMS will be applied to fmd the 
reconunended sensor architecture and communication protocols to minimize the communication cost and 
yield die best sensory information. 
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The foDowing materials enUtled "DISTRIBUTED MICRO FLOW-SENSOR 
ARRAYS AND NETWORKS: DESIGN OF ARCHITECrURE AND FAULT-TOLERANT 
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ABSTRACT 



TTie measurement of fluid flow is an important technology used in neariy every industry. 
Flow-sensors measure the movement of fluid flow (liquid or gas flow). Conventionally, 
single sensor systems are used to save space and cost The development of MEMS 
enables production of large amounts of micro sensors at low cost. Distributed micro flow- 
sensor arrays and networks (DMFSA/DMFSN) are built from collections of spatially 
scattered, intelligent micro flow-sensor nodes. Each node is able to measure the local 
flow within its accuracy limits, process the raw sensor data, and cooperate with its 
neighboring nodes. Multi-sensor architectures can improve the reliability and fault- 
tolerance of the system. The objective of this research is to address the design of the 
network architectures, communication protocols, and sensor information integration 
algorithm, for DMFSA/DMFSN, by developing a MEMS sensor network simulation tool, 
called TEE (Teamwork Integration Evaluator)/MEMS to analyze the efficiency of the 
sensor network architectures and conamunication protocols. New cluster network 
architecture, two versions of fault-tolerant time-out communication protocol (FTTP) - 
one with blind base stations, and one with smart base stations, and a fault-tolerant sensor 
integration algorithm (FTSIA) are proposed and analyzed in this thesis. TDE/MEMS is 
applied to evaluate five combinations of network architectures and communication 
protocols and compare the two versions of FTTP. The experiment results show that under 
certain conditions the proposed cluster architecture can support efficient information 
conmiunication and integration in DMFSA/DMFSN because it requires relatively low 
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communication cost and tolerates possible sensor and link failures. The results also diow 
that although the FTTP with smart base station requires an additional knowledge base in 
its structure compared with the one with blind station, it significantly reduces the 
communication time. In addition, in order to evaluate the effectiveness of the proposed 
FTSIA, three cases of simulation experiments were conducted to compare the output 
results from FTSIA and statistical calculations. The experiment results indicate tiiat the 
FTSIA yields more accurate measurement results tiian non-FTSIA statistical calculation, 
when some of tfie sensors yield faulty information. An empirical experiment with an 
array of eight pressure sensors was developed to illustrate FTSIA. General design 
guidelines are developed based on the experiments. 
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CHAPTER 1 INTRODUCTION TO THE PROBLEM 



The measurement of fluid flow is an important technology used in nearly every industry. 
Flow-sensors are used to measure the movement of fluid flow (liquid flow or gas flow). 
They are widely used in a variety of applications, such as automotive applications, 
environmental monitoring, biomedical applications, etc. 



Before the advent of micro minimization, sensor systems tended to be bulky and 
expensive, so single sensor systems were used to save space and cost. AIthoug|i single 
sensor systems are relatively easy to construct and analyze, they have at least two 
disadvantages: 

1, They have limited applications and usages. For instance, if several variables, say, 
temperature, pressure and flow density, need to be measured at the same time, 
single sensor systems cannot be used. 

2. In addition, a sinjgle sensor cannot guarantee delivery of accurate information all 
the time, because it is inevitably affected by noise or other uncertain disruptions. 

Thus, single sensor systems are not suitable for critical applications. 



Microelectromechanical system (MEMS) technology deals with design and fabrication of 
entire electrical and mechanical systems, usually on a single silicon chip. The 
development of MEMS enables production of large amounts of micro sensors at low cost. 
Micro flow-sensors are one of the most common MEMS devices. 



LI. Disadvantages pf Single SensQr gy^tems 



1.2. MEMS Flow-Sensors in Distributed Arrays and Networks 
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Due to their small sizes, the advantages of micro flow-sensors are: 

1. Interfere less with the environment they are measuring 

2. Lower manufacturing cost 

3. Can be ^plied in narrow spaces, such as inside living organisms, small pipes, 
automobile engines, etc. 

4. Multiple sens(MS can be used with redundancy to improve information accuracy 
and fault tolerance of the system. 

Distributed micro flow-sensor arrays and networks (DMFSA/DMFSN) are built from 
collections of spatially scattered, intelUgent micro flow-sensor nodes. Each node has die 
ability to measure Ae local flow within its accuracy Umits, process the raw sensor data, 
and cooperate with its neighboring nodes. Sensors incorporated with dedicated signal 
processing functions are called intelligent sensors, or smart sensors. The main roles of 
dedicated signal processing units are to enhance design flexibiUty and realize new 
sensing functions. Additional roles are to reduce loads on cratral processing units and 
signal transmission lines by distributing informatitm processing to the lower layers of die 
system (Yamasaid, 1995). In this diesis, the micro flow-sensor array refers to a group of 
the same type of micro flow-sensors which are placed close together to measure the same 
variable of interest. Different groups of micro flow-sensor arrays (either of the same type 
or different types) can be distiibuted widely in the flow environment to form a micro 
flow-sensor network. 



1.3. Problem Definition 
Developing effective DMFSA/DMFSN faces great challenges. Five major design issues 
are: 

1. Suitable micromachining technique. To achieve die realization of micro flow- 
sensor structure, the structure has to be released using suitable micromachining 
techniques (Toshio and Fumihito, 1999). 
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2. Appropriate network architecture, NA. Since micro sensor nodes are distributed 
spatially, appropriate sensor network architecture has to be found to support 
efficient conununication among the sensor nodes. 

3. Cost-efficient communication protocol, CP. Communication cost C is often 
related to communication time, c„ and energy consunn|>tion, c„ of the system. 
Minimizing communication cost is important in the online control and energy- 
constraint aiq>lications. 

4. Fault-tolerant conununication protocol and information integration algorithm, lA. 
It is unrealistic to expect all the micro sensor nodes and communication links 
along the path to be fault-free all the time because of the constraints of micro 
sensors and complex flow environment. 

5. Optimal placement and deployment of micro flow-sensors. Flow measurement is 
a complex process because of the numerous factors that can affect accuracy, such 
as flow stability, non-homogeneous flow, etc. Hence, the placement and 
deployment of micro sensors in the flow space is important in order to minimize 
the error due to the uncertain factors. 



1.4. Scope and Objectives 
The objective of this research is to address the issues No. 2, 3 and 4 in section 1.3 in the 
design of the network architectures, NA, communication protocols, CP, and sensor 
information integration algorithm, lA, forDMFSA/DMFSN, Issues No. 1 and 5 are out of 
the scope of this thesis. 

In Chapter 2, some recent research is reviewed on the design of micro flow-sensors, 
sensor network architectures, communication protocols, and multi-sensor fusion and 
integration. The proposed cluster architecture and fault-tolerant time-out communication 
protocol (FITP) for DMFSAA5MFSN are presented in Chapter 3. In order to evaluate 
the sensor network architectures and protocols developed for DMFSA/DMFSN, 
TIE(Teamwork Integration Evaluator) /MEMS is developed. In Chapter 4, the 
background and structure of TIE/MEMS are described. In addition, two sets of simulation 
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experiments were conducted to iUustrate the application of HE/MEMS: one is to evaluate 
the efficiency of five combinations of sensor network architectures and communication 
protocols based on the communication time and energy consumption; and the other is to 
evaluate and compare two alternatives of FTTP, FTTP with bhnd base station and with 
smart base station. Those experiments were conducted on SGI Origin 2000 workstation at 
the Department of Computer Science in Purdue University. In Chapter 5, the proposed 
fault-tolerant sensor integration algorithm (FTSIA) is presented. Three cases of 
simulation experiments were conducted to test the effectiveness of this algorithm. In 
these experiments, the results from FTSIA were compared with the statistical mean of the 
same simulated measurement values. In chapter 6, a case study of measuring the airflow 
using a pressiire sensor array is described to illustrate the application of FTSIA in real 
fiow measurement. The experiment setup and results are also described. Finally, the 
design recommendations of DMFSA/DMFSN are summarized, and the conclusions and 
future research are discussed in Chapter 7! 



1.5. Contributions and Sipnifica nce of die Reseair^ 
This research is being undertaken in response to the real problems encountered in the 
flow measurement area. Given the widespread use and importance of flow-sensore, the 
limitations of the current flow measurement technology are suriMising. DMFSA/DMFSN 
opens original new directions for this area. Some of the main contributions of this 
research are: 

1. Ditelligent Sensors manufactured using advanced IC (Integrated Circuit) 
technology will fundamentally change the nature of sensing and control systems. 

2. The ability to place the micro flow-sensors when and where they are needed 
enables better control of the process. 

3. Redundancy enables more reliable monitoring and control because occurrence of 
failures of some sensors will not inhibit the sensing ability of the system. 
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CHAPTER 2 LITERATURE REVIEW 

This Chapter reviews the recent research on the design of niicro flow-sensors, sensor 
network architectures, communication protocols, and multi-sensor fusion and integration. 

2. L Micro Plow-Sensor Design 
Realization of micro flow-sensors can rely on several, different operating principles. The 
most conunon ones include: 

L Anemometer 

A thermal anemometer refers to sensors that measure total heat loss. The heat loss is due 
to some heating element (such as a polysilicon resistor) heating the fluid. The heal loss 
depends on the flow rate of the fluid and increases or decreases with the increase of flow, 
depending on the operating mode for the smicture (Rasmussen and Zaghloul, 1998). 
Anemometer is the most commonly used type of micro flow-sensor (For instance. 
Nguyen and Kiehnscherf, 1995; Mayer et al., 1997.) Figure 2.1 shows the mass flow- 
sensor designed by Nguyen and Kiehnscherf (1995). This thermal flow-sensor contains 
polysilicon heaters and measurement resistance on a diaphragm. The fluid flows in an 
anisoiropically etched micro-channel with a cross-section of 0.6mm^ and a length of 
10mm. The channel is covered with a Pyrex glass plate by anodic bonding. A small 
thermal response time from 1 to 4ms is reached with a diaphragm thickness from 10 to 
30|im. 
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Detection 
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Fluid 
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Contact 
Pad^ 




Cut-out Membrane '^'^^ 
Figure 2. 1 Structure of the Mass Flow-sensor Designed by Nguyen and Kiehnscherf (1995) 



2. Time-of-Flight 

A good method to measure flow would be by marking the flowing medium and detecting 
the movement of the mark. For example, a mark can be a floating object which can move 
with the medium while being stationary with respect to the medium. The time that it 
would take the object to move with the flow from one position to another could be used 
for the calculation of the flow rate (For example, Shofi and Esashi, 1994; Van Kuijk et al., 
1994.) 

Both anemometer and sensors using the thermal time-of-flight principle belong to thermal 
transport flow-sensors. They are more sensitive than other types of flow-sensors and have 
a broad dynamic range. They can be employed to measure very minute gas or liquid 
displacements, as well as fast and strong currents. Major advantages of these sensors are 
the absence of moving components and the ability to measure very low flow rates. 

3. Pressure Gradient 

A fundamental equation in fluid mechanics is the Bernoulli Equation, which is strictly 
applicable only to the steady flow of nonviscous, incompressible medium. With the 
Bernoulli Equation the fluid velocity can be found by measuring pressures along the flow. 
The pressure gradient flow measurement requires introducing flow resistance. The 
concept is analogous to Ohm*s law: voltage (pressure) across a fixed resistor is 
proportional to current (flow). The advantages of the pressure gradient method include 
the absence of moving components and use of the standard pressure sensors which arc 
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readily available. One of its disadvantages, however, is the restriction of flow by resistive 
devices. An example of the pressure flow-sensor is shown in Berberig et a]. (1998). 

4. Ultrasonic Sensors 

Flow can be measured by employing ultrasonic waves. The main idea behind this 
principle is the detection of frequency or phase shift caused by the flowing medium. One 
possible implementation is based on the Doppler Effect. According to the Doppler Effect, 
when the target moves toward or away from the transmitting antenna, the frequency of the 
reflected radiation will change. An alternative method relies on the detection of the 
increase or decrease of the effective ultrasound velocity in the medium. The effective 
velocity of sound in a moving medium is equal to the velocity of sound relative to the 
medium plus the velocity of the medium with respect to the source of the sound. The 
difference between the velocity of the sound wave propagating upstream and the one 
propagating downstream is exactly twice the velocity of the medium, which allows one to 
determine the velocity of the flow. Several researchers have implemented the ultrasonic 
flow-sensors (For example, von Jena and Magori, 1995; Ishikawa el al., 2000.) An 
advantage of the ultrasonic flow-sensor is its ability to measure flow without direct 
interference with the fluid. 

5. Physical Deflection 

One example of a micro flow-sensor using the physical deflection principle is the device 
designed by Gass et al. (1993). It measures the drag force exerted by the fluid onto an 
obstacle - a cantilever beam having overall dimensions of 0.03mm x Inun x 3mm. This 
sensor can detect very low flow (down to nanoliter resolution). 

2.2. Sensor Network Architecture 
The abstraction of a sensor network can be represented as AW = (V. £) , where V is the set 
of sensor nodes and E is the set of communication links connecting the sensor nodes, 
V = ](i = 1,2,,.., N^), where v. is the sensor node i and Ns is the number of sensor nodes 
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in the networic.E:={e^}(i «U,„.,Ar,;;=lA...,iV,), where e^is the communication link 
between v. and Vj . 

Two most conventional network architectures that can be used for sensor networks are the 
committee architecture (CO A) and the hierarchical architecture (HIA), first discussed by 
Wesson et al. (1981). In the COA (Figure 2.2(a)), each node in the network is autonomous 
and can broadcast information to all other nodes in the network. Thus, if M is equal to n, 
it needs interconnections that are of the order 0(n^). With such a high number of 
interconnections, this architecture is robust relative to link failures because the messages 
from a source node to a destination node can be easily rerouted by bypassing the failed 
links- The COA, however, is not suitable for a sensor network with a large value of iV„ 
because it is too expensive with respect to communication cost, C, and not easily 
extensible. In the HIA, nodes are connected to form strict hierarchies. Each node (except 
the lowest-level leaf nodes) receives information from the next lower level node, 
integrates the information received according to its position in the hierarchy, and then 
sends information to the node at the next upper level. The root node is the node at the 
highest level, which is responsible for interpretation of data at the global level and for 
decision-making. For a network with equal to /i. HIA requires Oin) interconnections. A 
special case of HIA is a complete binary tree (Figure 2.2(b)). HIA has several advantages, 
including constant node degree and easy extensibility, but it may not produce accurate 
information at higher levels because the errors can accumulate when the information is 
transmitted up the hierarchy. In addition, HIA is not tolerant to link failures. 
Two other sensor network architectures, proposed by Iyengar et al. (1990) and Iyengar et 
al. (1994). are the Flat Tree architecture (FTA) and Multilevel binary de Bruijn 
architecture (MBDA). In FTA, the network is partitioned into clusters and each cluster is 
internally organized as a binary tree. There are three types of nodes in each cluster leaf 
nodes, intermediate nodes and root nodes. Leaf nodes are at the first level of the netwoik 
and the parent of a node is at one level higher than its children. All the root nodes in the 
network are fully interconntected (Figure 2.2(c)). FTA is not robust to link failures, 
MBDA is a multilevel network (MLN) with the top level fully connected and with each of 
the other levels interconnected as a de Bruijn network. MLN is a network in which each 
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node of the network can be associated with a level number. Assuming that a node i is at 
level m, then the neighbors of node i are the set of nodes at the same level, that is m, to 
which node i is connected. The parents of node / are the set of nodes at level m-l to which 
node I is connected. The children of node i are the set of nodes at level m+l to which i is 
connected. If each node in the network can have at most one parent and r children, the 
network is called r-ary multilevel network. Definitions and more details about the de 
Bruijn graph and network can be found in Leighton (1992). MBDA is not suitable for the 
DMFSA/DMFSN because of its complex structure. 





010* 




Level 2 



110 



111 

(c) (d) 
Figure 2.2 Network Architectures 

(a) Committee Architecture. COA (b) Hierarchical Architecture. HIA 
(c) Flat Tree Architecture, FTA (d) Multilevel Binary de Bruijn Architecture, MBDA 



23. Communication Protocol 
In recent years, there has been increasing interest in the development of communication 
protocols for the distributed micro sensor network, especially the wireless micro sensor 
network (for example, Utete and Durrant-Whyte. 1994; Gossink et al., 1998; Agre et al.. 
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1999; Chandrakasan et al., 1999; Heinzelman et al., 1999; Bhushan and Rengaswamy, 
2000.) The sensor network can improve the accuracy of infonnation obtained via 
collaboration among sensor nodes. 

The simplest communication protocol is the point-to-point protocol (PIP). In this 
protocol, sensor nodes A and B communicate exclusively with each other without 
interfering with other nodes (Figure 2.3(a)). 

Another commonly used communication protocol is the broadcasting protocol (BP). In 
this protocol, the information sent out by the sender node is received by all the other 
nodes within a certain range of the sender node (Figure 23(b)). The advantage of this 
protocol is that when a single node sends information out to a broadcast address, it can 
reach aJi its neighboring nodes using a single transmission. BP is. therefore, ^ cheaper 
than FTP protocol in one to-many communication, where n is the number of neighboring 
nodes for the sender node. BP is a common type of the classic flooding protocol, which 
means that upon receiving a piece of data, each node then stores and sends a copy of the 
data to all of its neighbors. It disseminates data quickly in a network where bandwidth is 
not scarce and links are not loss-prone. However, since a node always sends data to its 
neighbors, regardless of whether or not the neighbor has already received the data from 
another source, it leads to the implosion problem and wastes resources by sending 
duplicated copies of data to the same node. 

Gossiping protocol (GP) is an alternative to the classic flooding protocol, in which 
instead of indiscriminately sending information to all its neighboring nodes, each sensor 
node only forwards the data to one randomly selected neighbor (Heinzelman et a/., 1999) 
(see Figure 2.3(c)). While the GP distributes information more slowly than BP. it 
dissipates resources, such as energy, at a lower rate as well. In addition, it is not as robust 
relative to link failures as BP because a node can only rely on one other node to re-send 
the information for it, in case a link failure happens. 
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(a) (b) (c) 

Figure 2.3 Communication Protocols 
(a) Point to Point Protocol, FTP (b) Broadcasting Protocol, BP 

(c) Gossiping Protocol, GP 



In order to solve the problem of implosion, a sensor protocol for Information via 
Negotiation (SPIN) was proposed by Ifeinzelman et al* (1999). SPIN nodes negotiate with 
each other before transmitting data, which helps ensure that only useful information will 
be transferred. SPIN nodes use three types of messages to communicate: ADV (new data 
advertisement), REQ (request for data) and DATA (data message). SPIN protocol works 
in three stages (ADV-REQ-DATA). The protocol starts when a node advertises the new 
data that is ready to be disseminated. It advertises by sending an ADV message to its 
neighbors, naming the new data (ADV stage). Upon receiving an ADV, the neighboring 
node checks to see whether it has already received or requested the advertised data. If not, 
it responds by sending an REQ message for the missing data back to the sender (REQ 
stage). The protocol completes when the initiator of the protocol responds to the REQ 
with a DATA message, containing the missing data (DATA stage). Figure 2A shows an 
example of this protocol. 

In a relatively large sensor network, there is typically no central base station to which all 
nodes can directly communicate, and thus a low-power routing protocol is needed to route 
the data from the sensors to the high-powered base-station. An energy efficient protocol 
for the networked sensors was proposed in Chandrakasan et al. (1999). According to this 
protocol, each node maintains a routing table where die cost of each path is the power 
dissipated in transmitting information along diat path. Therefore, each node can transmit 
data to the high-powered node using a minimum amount of total energy. Furthermore, 
most nodes in the sensor network "go to sleep" except for a few "monitor-site" sensors 
that are chosen randomly and reassigned periodically. These monitor sensors run simple 



algorithms to detect the presence of an "event". When an event is detected, each monitor- 
site sensor "wakes up" its neighbors, which in turn 'Vake up* their neighbors. 




(1) (2) (3) 




(4) (5) (6) 



Figure 2.4 An Example of the SPIN Protocol (Heinzelman et al., 1999) 
(1) node A starts by advertising its data to node B (2) node B responds by 
sending a request to node A; (3) node A sends data to node B; (4) upon 
receiving the data, node B sends out advertisement to its neighbors; (5) node 
B receives the requests ftom its neighbors and sends data to them 



2.4. Multi-Senftor Fusion an d Integration 
A sensor system must interact with its environment. Measurement errors (such as noise) 
and sensor failures of single devices are unavoidable. Multi-sensor Fusion and Integration 
is the process of combining information from different sensors. Using readings from 
several independent sensors reduces the vulnerability of a system to the failure of a single 
component. Combining readings from several different kinds of sensors in a system can 
result in more accurate information than otherwise possible. Possible applications of 
multi-sensor fusion and integration include virtually all systems involving signal 
processing, such as in aeronautics (Brown and Olson, 2000), manufacturing (Kuo. 1998), 
remote sensing (Paul, 2001), medical application (Kubler et al., 2000), etc. 



40 



13 



There are several types of multi-sensor fusion and integration depending on the types of 
the sensors and their deployment (Iyengar et al., 1995): 

1. Competitive Integration: Each sensor ideally measures identical information, but 
in reality, is subject to noise and other disruptions. In this type of integration, even 
if some sensors are faulty, the infomiation from the erroneous sensors can be 
discarded during the integration process if a good integration algorithm is used. 

2. Complementary Integration: Partial and overlapping (complementary) information 
is available from each sensor. The partial information must be integrated into 
complete information. If some sensors are faulty, the result of the integration can 
be erroneous. | 

3. Cooperative Integration: Partial and non-overlapping (cooperative) information is 
available from each sensor. The partial information must be integrated into 
complete information. No sensor can be faulty; otherwise, the results would be 
faulty, 

4. Independent Integration: Unrelated information is available from the sensors. 
Integrating them simply means adding more fields to a data record. 

5. Temporal Integration: It requires integration of information over a given time. The 
data, sampled at time intervals, must be processed to re-create the notion of time 
sequence. 



6. Spatial Integration: It requires integration of information over a given space. 
EHfferent mathematical tools have been used in the research on multi-sensor fusion and 
integration. The most conomonly used are (Brooks and Iyengar, 1998): 

L Algorithm: An algorithm may be regarded as a finite sequence of instructions to 
accomplish a given task. In addition to usually having input and output, an 
algorithm has three properties: each instruction is definite, i.e., its meaning is 
unambiguous; each instruction is effective, i.e., it can be easily performed by a 
machine or a person; the sequence of instructions is finite, i.e., it terminates after a 
finite number of steps. 
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2. Linear Algebra: It is the branch of mathematics concerned with solving 
simultaneous linear equations. It provides tools for grouping together and solving 
sets of related equations. 

3. Probability: It is a methodology that has been developed for representing and 
quantitatively measuring random events. It is the source of most basic methods for 
dealing with uncertainty. 

4. Meta-Heuristics: It is a term used to characterize a number of methods that have 
been found to be practical algorithms (meaning, not guaranteeing optimality) for 
solving non-linear problems. It includes tabu search, genetic algorithms, simulated 
annealing, and artificial neural networks. 

In this thesis, it is assumed that all sensors are of the same type and measure the same 
variable of interest, therefore, competitive integration is the research interest here, A 
sensor integration algorithm is developed for DMFSA/DMFSN, because it is easily 
implemented in the processing unit with limited computation power. 
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CHAPTER 3 DMFSA/DMFSN ARCHITECTURE AND COMMUNICATION 

PROTOCOL 

In this chapter, the proposed cluster network architecture and fault-tolerant time-out 
communication protocol for DMFSA/DMFSN are presented. 

3. 1 > Cluster Network Architecture for DMFSA/DMFSN 
As mentioned in Chapter I, suitable network architecture is important to support efficient 
information conununication and integration in DMFSA/DMFSN. Figure 3.1 shows die 
proposed network architecture for DMFSA/DMFSN. Cluster architecture (CLA) is 
chosen because it overcomes the major limitations of both COA and HIA» by limiting die 
number of communication channels while at the same time tolerating possible sensor 
node and link failures. DMFSA/DMFSN are divided into sensor cluster units (SCUs), 
each of which consists of a set of intelligent micro flow-sensor nodes, and a base station. 
The intelligent sensor nodes within the same SCU are also called sibling nodes. Each 
intelligent micro flow-sensor node consists of a processing unit and an associated micro 
flow-sensor that measures the variables of interest. A sink node is needed for each cluster 
to integrate information from all the sensor nodes in the cluster. Considering the 
complexity of computation and/or energy constraint of the micro sensor node, a more 
powerful, dedicated processor is chosen as the base station for each cluster. The base 
station acts as the control node of the respective SCU, controlling the signal transmission 
of the sensor node, managing information rerouting in case of link failures, and 
integrating the information from all the sensor nodes in the SCU. In order to increase 
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the system tolerance to link failures, the sensor node sends information not only to the 
base station, but also to one or more other sibling nodes that are also called its backup 
nodes. Two communication protocols can be used among the sibling nodes: 

(1) BP. in which each sensor node broadcasts its information to all its sibling nodes; 

(2) Evenly Wide Gossiping Protocol (EWGP), in which instead of sending information to 
only one other sibling node by random selection (Heinzelman et al.. 1999)» the backup 
nodes are chosen in a way that the roles of backup nodes are evenly distributed among 
the sibling nodes. The intention is that no nodes are assigned too many 'duties', while 
others are *starved' (Liu et al., 2001c). 

The selection of protocol type and the number of backup nodes, //^, in EWGP depends on 
the link failure rate, r/, of the system. If n is relatively high, BP and EWGP with larger 
value of Ni, are preferred, because they have better fault tolerance. Otherwise, EWGP 
with smaller value of Ni, is better considered, because less communication is needed. The 
comparison of the conununication cost of BP and EWGP with Ni, equal to one is shown 
in section 4.2.1 in chapter four. The same protocol works for the conununication among 
the base stations. The external processor is the commander of the whole network, in 
which information from all the SCXJs is integrated for decision making as required. The 
clusters can be deployed along the flow path, depending on the flow situation and 
applications. 

Intelligent 
Sensor 



Sensor Cluster Unit 
(SCU) 



Intelligent 
Sensor 




Base 

^^Y^>Station 




External 
Processor 



Figure 3.1 Proposed Cluster Architecture for the DMFSA/DMFSN 
(Dash line means the link may or may not exist depending on the protocol used) 
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3.2. Fault-tolerant Time-out Compunic ation Protocol (FTTP) 
The design of this new protocol is motivated by the combination of time-out and task 
coordination protocols (Esfarjani and Nof, 1998), and the need for fault-tolerant 
measurement integration. First, the base station broadcasts the data transmission request 
to all the sensor nodes in the SCU. Receiving the request, each sensor node sends its data 
to the buffer inside the base station. Each sensor node has a conresponding buffer in the 
base station to store its data so that the base station can trace the source of the data. Then 
the data will be retrieved into the processing unit in the base station. Fault-tolerant Time- 
out Communics^on Protocol (FTTP) is an improvement over the conventional time-out 
protocol by incorporating it with the fault-tolerant property (Liu and iNof, 2001b). 
According to the time-out protocol, the base station will stop waiting for the information 
from the sensor node if a certain amount of time, T, has passed, where T is the pre- 
determined time-out period and its value depends on the application. Base station will 
then announce the rerouting task message to its sibling nodes. Two alternatives exist for 
where the base station should send the task announcement: (1) with a Blind Base Station, 
the message is sent to all the sibling nodes; (2) with a Smart Base Station, the message is 
only sent to the sibling node that can finish the task earliest, based on the cwrent status of 
the nodes in the system. 



3.2.1. FTTP with Blind Base Station 
In this protocol, when the link of a sensor node, say node , fails, the base station will 
broadcast the rerouting task message (RTM) for v^to all its sibling nodes. Receiving 
RTM, its backup node (the tth backup node of ) generates a bid value dv«. bki € BKt 
, where BKk is the set of backup nodes of Vj .bv^cmbc calculated as follows: 

bv^^{m^^l)rfti (1) 
where ^is the transmission rate of node bid (time spent in transmitting one message) and 
ntki is the number of tasks that have been offered to node bti. Because of the computation 
limit of the intelligent sensor node, it is unrealistic for it to keep track of its u^smission 
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rate over time, }fcis assumed to be constant for node bu. The base statim then evaluates 
the priorities of the backup nodes based oa their bid values (the lower die bid value, the 
higher the priority), and sends die task offer message (TOM) to die backup node wiUi die 
highest priority, btfi). Receiving die task offer message, bta, sends die required backup 
data message (BDM) to the base stadon. 

The structures of die base station and intelligent nodes, and die logic chart of die FTTP 
widi blind base station are shown in Figure 3.2 and Figure 3.3 respectively. 




B|, B2 . Bj are the buffers for the corres|x»iding 
intelligent sensor nodes 



Rgure 3.2 Structures of die Base Station and Intelligent Sensor Nodes for die FTIP widi 
Blind Base Station 

(Dash line means die link may or may not exist depending on die protocol used) 
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Base Station 

Broadcasts rerouting task message 
(RTM) for to all its sibling 
sensor nodes 



Receives tv^i 



Determines the priorities of 
sensor nodes 



Sends task offer message (TOM) 
to the sensor luxie with the 
highest priority 



Receives BDM 



RTM 



TOM 



BDM 



Intelligent Sensor Node 
Receives RTM 




I Yes 

Generates bid value, i^v^ 
Sends bv^ to base station 




Yes 
Receives TOM 



Sends the required backup data 
Message (BDM) to base station 



Figure 3.3 Logic Caiart of the FTTP with Blind Base Station 



3.2.2. FTTP with Smart Base Station 
In this protocol, in order to limit the unnecessary communication and reduce the 
communication traffic, instead of broadcasting RTM for , whose link fails, to all its 
sibling nodes, the base station first evaluates die priorities of the its backup nodes and , 
then sends TOM only to the node that can finish the task the earliest. The time when the 
backup node bju finishes the task, tu, can be calculated as: 

(2) 

where m^. is the number of tasks that have been offered to the backup node fc^ and Yu is 
the transmission rale of bjti (time spent in trmsmitting one message). The base station 

(.1 
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then sends TOM to the node Bki, whose is the smallest In oider to find out , 
however, the base staticm needs to have a knowledge base which Iceeps track of the 
information for each node - its backup nodes, transmission rate, and the number of tasks 
diat have been offered. Figure 3.4 shows die structures of base station and intelligent 
sensor nodes for the FTTP with smart base station. The model of the knowledge base is 
shown in Figure 3.5, in which n is the number of sensor nodes in the cluster and the 
backup nodes are hypothetical for illustration only. Figure 3.6 shows the logic chart of 
the FTTP with smart base station. 

In comparison, relative to the FTTP witti blind base station, the FTTP with smart base 
station requires less communication, and has Ae advantage of being able to keep track of 
the updated tr^ismission rate of each sensor node, but it aiso re<joiites a base station with 
more complex structure (with an additional knowledge base). Therefore, a tradeoff has to 
be considered in order to determine when to q)ply each protocol. The experimental 
results of the comparison of C for the FTTP with bUnd base station versus the FTTP with 
smart base station are analyzed in section 4.2.2 in chapter four. 




Intelligent 
Sensor Node 



B|. B2, Bj are tt» buffers for the correspondiqg 
inteUig»M sensor nodes 



Figure 3.4 Structures of the Base Station and Intelligent Sensor Nodes for the FTTP with 
Smart Base Station 

(Dash line means the link may or may not exist depending on the protocol used.) 
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Fipitc 3.5 Modes of the Knowlaite* ®^^ for the FTTP with Smart Base Station 



Base Station 

Retrieves infcmnation firom 
Knowledge Base and calculates 
for node ba 



Sends task offer message (TOM) to sensor 
node bt(if whose r^;; is the smallest 



TOM 



Receives BDM 



BDM 



Intelligent Sensor Node 




Receives TOM 



Sends the required backup data 
information message (BDM) 



Rguie 3*6 Logic Chan of the FTTP with Smart Base Station 
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CHAPTER 4 TIE (TEAMWORK INTEGRATION EVALUATOR)/ MEMS 

TIE (Teamwork Integration EvaIuator)/MEMS is developed in this thesis as a new 
version of the TIE simulators developed by the PRISM group (Khanna and Nof, 1994; 
Huang and Nof, 1996; Anussomnitisam and Nof, 2000.) The objective of TIB/MEMS Is 
to simulate and evaluate the architectures and communication protocols developed fw the 
distributed micro sensor network. Tffi/MEMS is programmed with MPl (Message 
Passing Interface). MPI is not a new programming language; rather it is a Hbrary of 
subprograms that can be called from C and Rwtran 77 programs. Since its completion in 
June 1994, MPI has received widespread acceptance because it is based on message 
passing, one of the powerful and widely used paradigms of programming parallel systems 
(Snir et al., 1996). MPI can be used to simulate different communication schemes among 
the sensor nodes, such as point-to-point communication, collective communication, group 
communicaHon, and so on. In addition, implementations of MPI on top of the standard 
Unix inter-processor communicaUon enable portability to multiple processors, 
woricstation clusters, and heterogeneous networks of workstations (NOW) (Liu and Nof, 
2(X)la). 



4.1. Structure of TIE/MEMS 

The overall structure of TIE^MS is depicted in Figure 4.1. The inputs of TEE/MEMS 
are: 

(1) micro sensor network architecnire, NA; 

(2) communicaticm protocol, CP; 

(3) sensor network parameters (such as the number of micro sensor nodes, iV„ the number 
of sensor noctes in each cluster, N^, the number of failed links, Ni, etc.); 

76 
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(4) constraints unique to the given MEMS sensors, e.g.. physical and chemical constraints 
causing micro-stress, proximity noise, and other disturbances to measiuement and 
communication. 

TIE/MEMS consists of a collection of libraries, each of which is built to simulate one 
type of sensor networic architecture and communication protocol. After the requued 
inputs are entered, the corresponding library will be called to execute the program. The 
output of TEE/MEMS is the communication cost, C, which is represented by the 
communication time C|, or/and energy consumption If C is too high, NA, or CP, or 
both need to be revised. TTE/MEMS can be extended by adding libraries for the newly 
developed NA and CP. Table 4.1 lists some examples of constraints related to MEMS 
sensor operations. 




t 

Conununication 
Cost C(c^c J 



Yes, revise NA 




Yes, revise CP 



Effective W4 
and CP 



Figure 4. 1 Structure of TIE/MEMS 
(Dashed line - future development) 
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Table 4.1 Examples of Constraints Related to MEMS Sensors 







Scaling issues 


Viscosity force; Surface effect 


Problems in 
micromachining 


Squeeze film damping; Particular contamination; Static stfction 


Response 


Decrease of response time due to protective gels, coatings or 
diaphragms 


Reliability of 
microelectronic 
device 


Brteltdown of a wire-bond to a junction; Defective encapsufatlon; 
Defective semiconductor material; Thermal runaway 


Ufetime 


Aging; Contamination of the device; Limited power supply 


Noise/interferences 


Degraded signal due to the variance of temperature, humidity, etc. 



4.2. Simulation Ex periments and Results 
In order to illustrate how TIE/MEMS is used to evaluate the efficiency of the network 
architectures and communication protocols discussed in the previous chapters, two sets of 
experiments were conducted. The experiments were conducted on the SGI Origin 2000 
workstation at the Department of Computer Science in Purdue Univereity. The flow 
diagram of TIE/MEMS programming structure is shown in Figure A.1 in Appendix A. 
The data results of the two experiment sets are shown in more detail in Tables A.1 - A.3 
in Appendix A. 



4.2.1. First Experiment Set 
In this experiment set, five combinations of network architectures and communicatiai 
protocols were simulated and compared They are HIA and FTP, COA and FTP, COA 
and BP. CLA and BP and CLA and EWGP with the number of backup nodes. Nb, equal 
to one. Thirty experiments were run f<w each combination and the outputs of the 
experiment are c, and c, per communication round. The results of c, vs. Ns and the 
normaUzed communicaUon time, c«, vs. Ns (normalization by the ratio to c, of HIA and 
FTP) are shown in Table 4.2 and Table 4.3, respectively. Figure 4.2 shows the results of 
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Cni VS. Ns graphically. Table 4.4 and 4,5 show the results of vs. Ns (assuming the 
transmitting cost is 600mv per message transmission and the receiving cost is 200mw per 
message reception) and the normalized energy consumption, c^, vs. Ng (nonnalization by 
the ratio to of the leaf node of in HIA) respectively. 



Table 4.2 Communication Time (Seconds) of Rve Combinations of Network 
Architectures and Connmunication Protocols vs. No. of Sensor Nodes 



AM 


CP 


N, 


3 


8 


15 


20 


30 


HIA 


PTP 


0.000102 


0.000441 


0.000782 


0.000902 


0.001674 


COA 


PTP 


0.000725 


0.006153 


0.022920 


0.040894 


0.093200 


COA 


Br 


0.000471 


0.002233 


0.004835 


0.009217 


0.020555 


CLA 


BP 


0.000138* 


0.000161* 


0.000192^ 


0.000205* 


0.000242* 


CLA 


EWGP 


0.000112* 


0.0001 18* 


0.000115' 


0.000115* 


0.000121** 



*Nc=3; 2Nc=4; ^Nc=5; 'Nc=5; ^Nc=6 



Table 4.3 Normalized Communication Time of Five Combinations of Network 
Architectures and Conmiunication Protocols vs. No. of Sensor Nodes 
(Normdization by the ratio to Communication Time of HIA and PTP) 



AM 


CP 


At 


3 


8 


15 


20 


30 


HiA 


PTP 


1 


1 


1 


1 


1 


COA 


PTP 


7.143842 


13.946657 


29.307476 


45.319050 


55.668379 


COA 


BP 


4.637767 


5.061201 


6.182039 


10.214362 


12.277506 


CLA 


BP 


1.355993' 


0.364110* 


0.245972^ 


0.227731* 


0.144386* 


CLA 


EWGP(Aft=1) 


1.103448' 


0.267473* 


0.147046' 


0.127443* 


0.07234' 
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3 8 15 20 30 

N. 



Figure 4.2 Normalized Communication Time of Five CombinaUons of Network 
Architectures and Communication Protocols vs. No. of Sensor Nodes 



Table 4.4 Energy Consumption of Five ComWaations of Networic Architectures and 
Communication Protocols (mw) vs. No. of Sensor Nodes 
(* Ce for the leaf nodes ; ** for the intermediate nodes ; *** c, for the root node; 
no intermediate nodes for HIA with when N, is equal to 3) 



AM 


CP 




3 


8 


15 


20 


30 


HIA 


FTP 


600* 


600- 


600* 


600* 


600* 


N/A 


800" 


1200" 


1400** 


1400- 


400*** 


800*" 


800*" 


800*** 


1600*** 


COA 


PTP 


4800 


44800 


168000 


304000 


696000 


COA 


BP 


2400 


15400 


50400 


87400 


191400 


CLA 


BP 


1000' 


1200^ 


1400' 


1400* 


1600* 


CLA 


EWGP (W4=1) 


800' 


800^ 


800' 


800* 


800' 
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Table 4.5 Nonnalized Energy Consumption of Five Combinations of Network Architectures 
and Communication Protocols vs. No* of Sensor Nodes 
(Normalization by the ratio to the energy consumption of the leaf node in HtA and FTP) 



AM 


CP 


N, 


3 


8 


15 


20 


30 


HIA 


FTP 


1* 


r 


r 


1* 


r 


N/A 


1.333333" 


2.000000" 


2.333333" 


2.333333** 


0.666667*" 


1.333333*" 


1.333333*" 


1.333333*" 


2.666667*" 


COA 


PTP 


8.000000 


74.666667 


280.000000 


506.666667 


1160.000000 


COA 


BP 


4.000000 


25.666667 


84.000000 


145.666667 


319.000000 


CLA 


BP 


1.666667' 


2.000000^ 


2.333339* 


2.333333* 


2.666667* 


CLA 


EWQP (N^\) 


1.333333' 


1.333333^ 


1.333333" 


1.333333* 


1.333333^ 



The eatperiment pesults indicate that when Hs is small, there is namuch difference in C 
(both C| and among the five combinations of NAs and CPs. In this case, COA and 
CLA are preferred because of their fault-tolerant characteristics. When Ns increases, 
however, C of the COA network, especially, the COA with PTP, increases dramatically. 
Therefore, HIA and CLA are favored for large sensor networks, CLA is a good choice in 
either case. 



4.2.2. Second Experiment Set 
In diis experiment set, FTTP with blind base station and FITP with smart base station, 
were compared. Nc and Ni are inputs of the experin^nt. Thirty experiments were run and 
the output is Ct per communication round. Table 4.6 summarizes the experiment results, 
which show, as expected, that the FTTP with blind base station requires less evaluating 
time, c^;, than the FTTP with smart base station. The reason for this is that in the FTTP 
with blind base station, the bid values are calculated by the intelligent sensor nodes 
locally before they are sent to the base station, and the base station only evaluates the bid 
values and chooses the sensor node with the minimal bid value. Much more infoimation 
transmission, including sending and receiving information between die sensor nodes and 
the base station, however, is involved in the FTTP with blind base station. Because the 
information transmission time, c^, is much longer than c^r which requires only 

7r 



28 

computation, c,, which is the sum of Cei and c„, of the FTTP with blind base station is 
longer than that of the FTTP widi smart base station. The difference of c, becomes man 
evictent when Ne and ^ increase. On the other hand, the strwture of the smart base 
station is more complex than the blind base station, as shown in Figure 3.2 and Hgure 
3.3. Figure 4,3 shows the ratio of c, of the FTTP with blind base station to the FTTP with 
smart base station, n», vs. Ne when Ni equal to 1, 2 and 3. 



Table 4.6 Communication Time, c„ of the FTTP with BHnd Base Station and FTTP 

with Smart Base Station vs. No. of Sensor Nodes in the Cluster, Ne , and No. 
of Failed Links, M 



FTTP 


N, 








^' 


FTTP with Blind 
Base Station 


1 


4 


0.00000217 


0.00014277 


0.00014493 


8 


0.00000270 


0.00026113 


0.00026383 


15 


0.0000029O 


0.00047737 


0.00048027 


21 


0.00000393 


0.00077307 


0.00077700 


2 


4 


0.0(»00313 


0.00034660 


0.00034973 


8 


0.00O0042O 


0.00058190 


0.00058610 


15 


0.00000587 


0.00114650 


0.00115237 


21 


0.00000740 


0.00170613 


0.00171353 


3 


4 


0.00000463 


0.00055680 


0.00056143 


8 


0.00000565 


0.00082713 


0.00083277 


15 


0.00000770 


0.00172153 


0.00172923 


21 


0.00001045 


0.00258752 


0.00259797 


FTTP with Smart 
Base Station 


1 


4 


0.00000433 


0.000(^43 


0.00002677 


8 


0.00000564 


0.00002891 


0.000)3455 


15 


0.00000716 


0.00003306 


0.00004023 


21 


0.00000817 


0.00004049 


0.00004866 


2 


4 


0.00000790 


0.00004533 


0.00005323 


8 


0.00000989 


0.00005448 


0.00006437 


15 


0.00001390 


0.00005850 


0.00007240 


21 


0.00001848 


0.00006068 


0.00)07914 


3 


4 


0.00001060 


0.00006213 


0.00007273 


8 


0.00001330 


0.00006696 


0.00008026 


15 


0.00001683 


0.00007069 


0.00008751 


21 


0.00002386 


0.00(W7291 


0.00009677 
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Figure 4.3 R^o of Ccmimunication Time of the FTTP with Blind Base Station to the 
FTTP with Smart Base Statim vs. No. of Sensor Nodes in the Qoster when 
Uie Number of Failed links Equal to 1» 2 and 3 
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CHAFIER 5 FAULT-TOLERANT SENSOR INTEGRATION ALGORITHM 

A system is said to be fauli-tol^t if its performance is not affected by faults in the 
system. A new fault-tolerant sensor integration algorithm (FTSL\) is developed in this 
thesis to combine information from multiple, same type of sensors in the cluster and 
produce ndiable results even if some sensors are faulty. FTSIA is appUed in the base 
node of a sensor architecture, such as the b^ station in the cluster network or the root 
node in the hierarchical architecture. Its purpose is to combine information from the 
sensor nodes in the respective cluster. The algorithm developed in this research is based 
on the ideas originated by Marzullo (1990) and later improved by Jayasimha (1996). 



S i. Concrete and Abstracf.^yps^j^ 
A concrete sensor (CS), is a device that measures the physical variable of interest in the 
environment, therefore, it can be also called a physical sensor. An intelligent sensor node 
is a CS. An abstract sensor (AS), is a piecewise continuous function that maps a physical 
variable into a dense interval [low. high] of physical values. TOs interval contains the 
value of the physical variable read by CS. The reason for using AS is that it is more 
reasonable to consider the reading from a sensor as a continuous set of values instead of a 
point value, because of the variation of die data due to noise, manufacturing, and other 
uncertainty issues. 

AS is derived from the point value received fix)m CS. If a CS. say v, , reads a value, say 
r, , and its maximum tolerable variation can be ±£,, then AS of v, can be defined as the 
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interval [Ij,t/, ]=[r, -^m^i which is caUed the interval estimate of the reading 
from V, , where U is the lower bound of v, , and Ui is the higher bound of v, . 

5.2. Some Definitions andTenninologdes 
i. An AS is correct if its interval estimate contains the true value of the measured 
physical variable, otherwise, it is faulty. If both AS Ai and A2 are coirect. they 
must intersect and the intersection contains the true value of the physical variable. 
An AS is said to be tamely faulty if it is not correct^ but its interval estimate 
overly with tfiat of a correct AS. A faulty AS that is not tamely faulty is said to 
\)tmdely faulty. 

n. A connected interval sequence, 5, is a sequence that has no overlaps or gaps, that 
is, S is made up of several intervals C/, C2, .... C*, and u^^ =/c^^^(l<i <*-!), 
where u^^ is the upper bound of interval C, and l^^ ^ is the lower bound of interval 
Cm . 

iii. Interval h completely overlaps h if Ia ^ h and ua > ub. 

iv, A distinct m-interval is an interval in which no more than m intervals intersect 

5,3. The Logic of PTSTA 
faput: Ns ASs {U Ui\ l^i^Ns, with at most/(either tamely faulty/, or widely faulty /h,) 
of them are faulty. / is a parameter determined by the manufacturer of the CSs from 
which the ASs were derived, based on their constraints. 
Output: A final output interval and a list of all the possibly faulty sensors 
Assumption: f <\N^ 

Pseudo PiDgram of FTSIA: 

I. Import Ns ASs, whose lower bound set is LB ^< Lff.^LBj,..., LB.v^., LB^ > and higher 

boundis HB =< HB^.HB, TO^ HB ^ >; 
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2. Sort the ASs by their lower bounds, and this sorted list Is t ^ l,.L,^ > and the 

coiteq»B<«ng upper bound set is « „<// I, • //•«•• 

3- Sort the first (Ns -/-I) intervals of 1/ and this sorted h'st is U ; 

4. Initialize the variables 

/ * ; /♦ ImtiaJize die interval set. in which each interval has three attributes, lower bound /, higher 
bound wand level of interval TO */ » • 

last :« 0 ; /♦ Initialize the number of intervals in / */ 

5. Main part of die algorithm 

J f ^own to 0 do 

5.1Insert(/^ ^ initsapiMropriatcplaceinf/ to ke^ £/ sorted ; 
52 Find the intersection of ( iv, - / ) intervals ^ 

^^//,-/ ^^0, ^^<o »s first clemftrt in t; */ 

5.3 Insert n := [L^,., . C/ j which is a ( A^, ^ / ^interval, into I and make appropriate changes to I 

5.3.1Locate the connected interval sequence 5 :.[C..C,_ C J in I in which each C, intends 
widi/V; 

5.3.2 Condition 1: If S is empty, insert )V to die tall / 
if(5 = 

tor := 1; 
else 

5.3.3 Condition 2: If diere is only one mterval d is S. find the intersection of N and C, 
if(S = [C,] ) 

if(i/(,j = )dien 

m,^ :=m^^ +1; 
else /♦if (t/^,) >«e,)*/ 

J ; 

:« + 1 ; 

y :=(/uiV,uiVj)-C,; 
end if 

else /♦if (/^^ < L„ ,^)V 



86 



^Af, AT,-/* 

/asr last + 1 ; 
else /*if > ii,^ */ 

/ojr :* Uut + 2 : 
end if 
end if 
end if 

else 

53.4. Condition 3: intervals Ca, Cp.i coinpletely overlaps N 
for Jfc 2 to p - 1 do 

end for 

else /*if <L^^^,Split[i^ ,11^ ] into two intervals*/ 

/fwf := tef/ + 1 ; 
end if; 

if )then 

'"at, := + 1 ; 
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tost last + 1 ; 
else /•if £/,„ ^ u^^ */ 

if({/,„ >«, )tfien 

end if 
end if 
end if 
end if 

5.35 Ronove Ae AS t ] from the list of ASs 

endftw 

6. Find out aU the possibly faulty sensors, print out the faulty sensor set and remove it from the lower 
bound set and higher bound set 

for I := 0 up to - 1 do 
for >:s Oopto AT, do 

if((«F/>= u, }ac{LB ^ ^ / » /♦senswy is possibty faulty*/ 

CLB := Lff - { CUB := UB - (UB , J; l*CLb and CVB are lower bound and higher bound 

sets that consist of absolutely correct sensors */ 

pnnt: sensor j is possibly faulty in interval i; 
end if 
end for 
endfor 

7. Find out the final output interval 

sort C£J? and flieswled list is CL =< CLy.CLt^,CL^ >: 
sort CffB and the sorted list is Cff =< CH ,,CH ^^..CH „ >; 
/, B[CL.,Cff J; is the final data interval estimate*/ 

The algorithm proposed by Maraullo (1990) only provides a single interval from all the 
sensors without fault detection. The algorithm proposed by Jayasimha (1996) provides all 
the (yV5-i)-interval (0<i:Sy). but in the flow measurement applications, single output results 
are preferred. The iinprovement of the FTSIA proposed in this thesis over the two 
algorithms described above is that it not only detects the possibly faulty sensors and 
surely faulty sensors but also generates a final data interval estimate from the correct 
sensors after removing readings of those faulty sensors. Examples of execution of FTSIA 

§2. 
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in (Ufferent cases are shown in Appendix B. Application of FTSIA in a real flow 
measurement is shown in Chapter 6. 

5.4. Simulation Experi^ nents and Resuttx 
Three experiment cases were simulated to test FTSIA. The mean from FTSIA, which is 
the average of the higher bound and lower bound of the final interval generated from 
FTSIA, w«c compared with the statistical mean of the same measurement values. Hie 
data acquired from the FTSIA and statistical calculations are also affected by the 
distribution of the errors of sensor readings. In order to evaluate the reliability of both 
methods, the worst situation was simulated in all three cases as described below. All the 
senscff nodes generate readings above Oie true value of the .measiued f^ysical variaWe, 
i.e., At, = - r* > 0 ( / = 1,2,..., N, ), where t * is the true value of the measured physical 
variable, r, is the reading received from the sensor node v, , and Ar, is the enor of the 
reading from v,. . The sensor readings were simulated using the random number generator 
function rand( ) in Microsoft Excel. In all three cases, assuming t* = 0.5 , and the 
maximum tolerable variation of the sensor reading, £ . is the same for all the sensor nodes 
ande = ±0.0)5. 

Shift of values from 0.5025 to 0.505 needs to be detected because the conect sensors 
should give readings no greater than 0.505. Therefore, the hypotheses are: 
Ho: M).505 
H,: n<0.505 

Based on initial experiments, variances of both mean from FTSIA and the statistical 
mean of iV, sensor readings are about 3.6x10'*. Assuming the type I risk is 0.025 and type 
n risk is 0.01, the sample size required to detect the shift is no less than one. If sample 
size is set to be 12 and type H risk is kept to be 0.01, then when the average of those 12 
samples is less than 0.502903. reject H©; otherwise, accept Ho (Hicks and Turner, 1999) 
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5.4.1 Case I 

In this case, all the sensors are assumed to be conect, geno-ating readings with mois 
witiiin ±1% of T* and their interval estimates contain t*. Rve experiment sets were 
conducted in this case: is equal to 8, 15, 21, 27 and 36. In all the experiments of this 
case. Ho is rejected for both mean from FTSIA and statistical mean of the sensor . 
readings. Kgure 5.1 shows the results finm FTSIA and die statistical mean vs. Nj. 




Kgure 5.1 Output Interval Estimate from FTSIA and Statistical Mean vs. No. of Sensor 
Nodes (7^.5, e=±0.005) 



5.4.2 Case 2 

In this case, some sensors are assumed to be tamely faulty, i.e., they generate readings 
with errors beyond ±1% but within ±2% of r*. Since the tamely faulty sensors intersect 
with some correct sensors, they can form ( iV, - / )-intervaIs that do not contain die true 
value. Because it is impossible to know in advance which interval is the correct niterval, 
all the possibly faulty sensors have to be deleted, and only those absolutely comet 
sensors will be kept. The absdutely correct sensors arc those overlapping with all the 
{N, -/ )-intervals. Following this procedure, however, some correct sensors will be 
considered as possibly faulty sensors and thus deleted. In ord» for die algoriUun to be 
effective in this case,/, must satisfy tiiat/, < because Acre will be no correct sensor 
left after deleting all die possibly faulty sensors when f,^-^N,, 
Twelve experiment sets were conducted in this case: AT, is equal to 8 with^ equal to 1, 2 
and 3; Ns is equal to 15 with/, equal to 1, 3, 5 and 7; and Ns is equal to 21 with/, equal to 
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1, 3, 5, 7 and 10. In this case, when is equal to 8 and// is equal to 1 and 2, Ho is 
rejected for the mean of FTSIA, but accepted for the statistic^ mean. When Nj is equal to 
8 and/, equal to 3, Ho is accepted for both the mean of FTSIA and statistical mean. When 
Ns is equal to 15 and^ is equal to 1, H© is rejected for both the mean of FTSIA and 
statistical mean. When Ns is equal to 15 and/, is equal to 2 and 3, Ho is rejected for the 
mean of FTSIA, but accepted for statistical mtan. When Ns is equal to 15 and f, is equal 
to 5 and 7, Ho is accepted for both the mean of FTSIA aQd statistical mean. When Ns is 
equal to 21 and/, is equal to 1 and 3, Hq is rejected for both the mean of FTSIA and 
statistical mean. When N, is equal to 21 and/i is equal to 5, Ho is rejected for the mean of 
FTSIA, but accepted for statistical mean. When Ns is equal to 21 and/, t$ equal to 7 and 
10, Ho is accepted for both the mean of FTSIA and statistical mean. The output interval 
estimates from FTSIA and statistical mean vs. /, when Ns is equal to 8, 15 and 21 are 
shown in Figures 5.2 -5.4 respectively. 
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Figure 5.2 Output Interval Estimates from FTSIA and Statistical Meaji vs. No. of Tamely 
Faulty Sensors (rM)-5, tefcO,005, iV^S) 
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Figure 5.3 Output Interval Estimates firom FTSIA and Statistical Mean vs. No. of Tamely 
Faulty Sensors (rM).5. £s=±0.005, N^15) 
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Kgure 5.4 Output biterval Estimates from FTSIA and Statistical Mean vs. No. of Tamely 
Faulty Sensors (rM).5, e^iOXtOS, lff^21) 



In this case, some sensors are assumed to be widely faulty, i.e., they generate readings 
with errors beyond ±2% of r*. In Older to control the experiments, the errors of the 
readings from the widely faulty sensors are assumed to be within ±10% of r*. Because 
when the readings from the widely faulty sensors are so close that their interval estimates 
intersect, they will form (N,-f )-intervaIs that deviate the results from the correct one. 

In order to prevent this situation, the number of widely faulty sensors, needs to satisfy 



Twelve experiment sets were conducted in this case: Ns is equal to 8 with/w equal to 1, 2 
and 3; Ns is equal to 15 with/w equal to 1, 3, 5 and 7; and Ns is equal to 21 vnthfy, equal 
to 1, 3, 5, 7 Md 10. 

In all the experiments in this case, Ho is rejected for the mean of FTSIA, but accepted for 
statistical mean. 

The output interval estimates from FTSIA and statistical mean vs.fw when Ns is equal to 
8, 15 and 21 are shown in Figures 5.5 -5.7 respectively. 



5 A3. Case 3 
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Figure 5.5 Output Interval Estimates ftom FTSIA and Statislical Mean vs. No. of Wictely 
Fault> Sensors (rM).5. e=±0.005. Nr=8) 
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Figure 5.6 Output Interval Estimates from FTSIA and Statistical Mean vs. No. of Widely 
Faulty Sensors (rM).5, £J=±0.005, NM5) 
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Figure 5.7 Output Interval Estimates from the FTSIA and Statistical Mean vs. No. of 
Widely Faulty Sensors(rM).5. e=±0.005, iV^21) 



81 



9 




5.5. Conclusions of FTSIA Experiments 



Analysis of Variance (ANOVA) was used to further analyze the results from the 
experiments in section 5.4. The mean from FTSIA and statistical mean were compaied. 
Table 5.1 summarizes the analysis results, which show that when fsO or /, is small 
compared with A/^ there is no significant difference between tfie mean from FTSIA and 
the mean from the statistical calculations. However, wben/r increases or /w>0, Ae mean 
fnmi FTSIA is significantly more accurate than the statistical mean. 
Table 5.2 shows how Ns and)? affect the mean from FTSIA in case 2. The reason for this 
result is that tamely faulty sensors can intersect wth some correct sensors and fwm 
( / >interval5 that do not contain the trae value, and the correct sensors that do not 
int^ect with the incorrect intervals will also be deleted as possibly faulty sensors. The 
correct sensors left after deletion usually have relatively higher errors because they must 
not only intersect with the correct sensors, but also with the tamely faulty sensors. For 
example, in case 2, when ^5=15 and/ =7, and when Nr=21 and/, =10. only one correct 
sensor is left after deleting all the possibly faulty sensors and this sensor has the highest 
sensor reading, i.e., the largest error, although still within the tolerable limits. Therefore, 
the performance of FTSIA will degrade gracefully when increases, but the final data 
interval estimate it generates will still contain the true value. Table 5.2 also indicates that 
in order to enable FTSIA generate more accurate results, the number of correct sensors 
left after deletion, /i^ should be at least half of Ns, i.e.. ^^N^ . Based on the logic of 
the algorithm. /tcCan be calculated as 



Therefore, when all the faulty sensors are tamely faulty in the system, in order to generate 
more accurate result using the FTSIA^iV, should satisfy 



For example, when/r=l and Ns > 4, wheny; =3 and Ns >12. and when /r=5 and Ns > 20, 
the FTSIA will yield more accurate results. 

The widely faulty sensors (case 3), however, do not cause die same problem as tamely 
faulty sensors. They do not intersect with die correct sensors and are simpler to detect and 



/i,=iV,~2Vr. 



(3) 



iV,>4V, 



(4) 
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delete, without affecting the conect sensors (see Table 5.3). Therefore, Ns and/, do not 
have significant effect on the results from FTSIA, as long as flie/n, k less than half of Af^ 
i.e., < . Therefore, when all the faulty soisois are wi(tely faulty in the system, N, 
should satisfy 

N,>2*f^ (5) 

Since Equations (4) and (5) are linearly independent, it can be concluded that when both 
tamely faulty and widely faulty sensors exist in the system, in order for FTSIA to 
generate more accurate results, should s^fy 

Ar,>4*/,+2*/. (6) 

The outputs of FTSIA also infer which sensw nocfes are possibly tamely faulty, and 
which are widely faulty. If some sensor node exists in the ail the possibly faulty sensor 
sets, it is widely faulty. On the other hand, if some sensor node appears in some of the 
possibly faulty sensor sets, but not all of them, it is possibly tamely faulty. 



5,6. FTSIA Algorithmic Characteristics 
The simulation experiment results show that the proposed FTSIA is effective because 
although when /is large relative to A^^, and the performance of the FTSIA may degrade 
gracefully, it always determines the closed interval that contains the true value. This 
algorithm has a complexity of 0{Ns\o%{Ns) because of the sorting in step 2 and 3 shown 
in the logic of the algorithm (the complexity of computing statistical memi is of 0(iV,)), 
This algorithm is simple and easy to implen^t (section 5.3). The results of three 
experimental cases are shown in more detail in Tables C.l -C.3 in Appendix C. 
Summary of the 29 experiments (with 12 treatments each) in all three cases is shown in 
Figure C.1 in Appendix C, 

Given the failure rates of the sensors and the required accuracy of the results, can be 
determined. 



Table 5.1 ANOVA Results (ass95%) of Effects of No. of Total Sensor Nodes and No. of 
Faulty Sensors on the Mean fifom FTSIA and Statistical Mean (Froin duee 
experimental cases) 
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Table 5.2 ANOVA Results (as4>5%) of Effects of No. of Total Sensor Nodes and No. 
Tamely Faulty Sensors on the Meai ftom FTSIA (Rom Case 2) 
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Table 5.3 ANOVA Results (0=^95%) of Effects of No- of Total Sensor Nodes and No. of 
Widely Faulty Sensors on the Mean from the FTSlA(Fn>m Case 3) 
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CHAPTER 6 CASE STUDY-FLOW MEASUREMENT WITH A PRESSURE SENSOR 

ARRAY 

This section illustrates the ai>plication of the FTSIA described in Chapter 5 through a 
case study - measuring the flow using a pressure sensor array. 

I 

6.1. Decision on the Number of Sensors 
Before the experiment was conducted, had to be decided to generate more accurate 
result using the FTSIA. The initial experiments showed that among the 8 available 
pressure sensors. /is equal to 2. Therefore, according to Equation (6). all 8 available 
pressure sensors had to be used to form a pressure sensor array in the experiment. 

6.2. Description of the Experim«>in 
This experiment was conducted in the PRISM lab at Pindue University. The experiment 
setup is shown in Figure 6.1. It consisted of a vacuum cleaner, which was used to suck air 
to go through the pipe connected to it, and a pressure sensor array, which had 8 pressure 
sensors. 8 holes were drilled on die pipe and signals were then input into the pressure 
sensors through die 1/8 Tygon Tubing. Because the input signals were very low Oess than 
0.01 V). differential inputs were used to reduce the noise effect. Data Acquisition (DAQ) 
board was used to change the analog signals into the digital signals which were input into 
the computer for data analysis. In this experiment, the DAQ used was KPCI-3107, 
manufactured by Keithley Instruments Inc. KPaOlO? offers some features that satisfy 
the experiment, such as maximum sample rate up to lOOKS/s, 16 single-ended or 8 
differential inputs. 12 gain ranges (1, 2, 4, 8. 10, 20, 40, 80, 100, 200, 

72. 
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400, 800) and compatible with LabView, etc. LabView is a graphical prograimning 
language that has been widely adopted throughout industry, academia. and govomnent 
labs as the standard for data acquisiticm, instniment control software, and analysis 
software (Bishop, 2001). Finally, the data retrieved from LabView were input into the 
FTSIA to generate a final data interval estimate ou^ut. 



Air 




DAQ 



imim 



LabView 



mm 



FTSIA 



Final Data Interval 
Estimate Output 

Figure 6. 1 Experiment Setup of Flow Measurement with A Pressure Sensor Array 
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^3. g^periment Results 
Table 6.1 shows an example of the data, t^, received from the pressure sensor array and 
the lower bounds, L,, and higher bounds, Uh of the data intervals, assuming £^=±0.0005. 
After the L, and i/, are input into the FSTIA, the outputs generated are as shown in Table 
6.2. The output results show that there is an 8-interval, which means the data from all 
sensors intersect, two 7-intervals and one 6-intervaL When there are tamely faulty sensors 
in the system, and if the lower bounds (higher bounds) of their data intervals are lower 
(higher) than the maximum (minimum) of the higher bouncte (lower bounds) of the 
conect data intervals, it is possible that the data from all the sensors intersect, so it cannot 
be concluded that there is no faulty sensor in the system. Therefore, the intervals of three 
possibly faulty sensors^ 4, 5, and 7 have to be removed fiom the interval set, and the final 
interval (the intersection of all the absolutely correct sensors) is [0.006747, 0,007685]. In 
addition, the results show that sensors 4, 5 and 7 are possibly tamely faulty because none 
of them exist in all of the possibly faulty sets and thus can be removed from the system 
for re-calibration or replacement. 



Table 6.1 An Example of the Data Received from 8 Pressure Sensors (volt) and the 

Lower Bounds and Higher Bounds of the Data Interval (Assuming te±0.0005) 
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Table 6.2 Summary of the Outputs from FTSIA (Data from Table 6.1 are inputs) 
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CHAPTER 7 CONCLUSIONS AND FUTURE RESEARCH 



In this thesis, three issues of design of the distributed micro flow-sensor aiiays and 
networks (DMFSA/DMFSN) were discussed, including the sensor network axchitectuie, 
communication protocol and sensor information integration algorithm. TIE/MEMS was 
developed to evaluate the efficiency of the network architectures and communication 
protocols. 



7.1. Design Recom mendations of DMFSA/DMFSN 
Figure 7.1 shows the overall logic of the design recommendations of DMFSA/DMFSN. 
In order to design DMFSA/DMFSN, first of all, the physical constraints of the micro 
flow-sensors used in the system have to be examined, such as the deterioration rate, 
link failure rate, 17, maximum number of failed sensors, f„ etc. Based on the 
information, the number of sensors required in the system, AT,, to provide more accurate 
results can be decided using Equation (6). 

Next, the type of network architecture, NA, and conununication protocols, CP, for 
DMFSA/DMFSN need to be chosen. If Ns is relatively small (depending on the response 
time and/or energy constraints in the applicaUon). COA or CLA is preferred. Otherwise, 
HIA or CLA is preferred. However, in both cases, if the flow environment is not steady 
along the path, which is always the case in the real flow measurement, CLA is better than 
the other alternatives because the clusters can be deployed along the path. ITie choices at 
CP depend on the sensor physical constraints. If rj and n in system are relatively hi^ 
(depending on the applications of DMFSA/DMFSN). BP and EWGP with larger number 
of backup nodes, Nb, are preferred. Otherwise, the EWGP with less Nt, is preferred. FTTP 
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with smart base station is preferred to FTTP with blind base station when the 
communication time, Cf, is critical in the application. Otherwise, if die base station widi 
simpler structure is required, FTTP with blind base station is preferred. These statements 
are based on the experimental results ftom section 4.2 in chapter 4. 
Finally, after the data from all die sensor nodes are received, FTSIA is applied in the base 
node and the final data interval estimate will be generated. Li addition, the FTSIA also 
outputs the list of all of possibly faulty sensor no<tes. 



£xaiiime sensor physical 
constraints (r^, n./^/^ etc,) 




Choose FTTP with 
Smait base station 



Choose FTTP with 
Blind base statton 



Input data from ail sensor 
nodes into the FTSIA 
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Final data interval 
estimate and list of 
possibly faulty sensors 



Figure 7.1 Overall Logic of Design Recommendations of DMFSA/DMFSN 
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7,2. Conclasions of Results 
MEMS technology is the merger of the IC worid with the mechanical worid. It allows the 
integration of the sensors and actuators with IC circuitry at an unpiecedented level. 
MEMS flow-sensor can be used in the applications in which the conventional sensor 
systems cannot or are not suitable to be used. DMFSA/DMFSN, built from a coUection 
of cooperating intelUgent micro-flow sensor nodes, can improve the reliability and fault- 
tolerance of the systeoL 

The appropriate sensor architecture for DMPSA^MFSN depends on the number of 
sensor nodes needed in the system to provide accurate information. When the number of 
sensor nodes in the system is relatively small. COA and CLA are better architectures 
because of thfeir fault-tolerant characteristics. Otherwise, HIA and CLA are preferred in 
order to limit the communication cost. For example. Table 4.3 shows that when AT, is 
greater than 15. C, of COA and BP is more than 10 times of C, of MA and FTP. and C, of 
COA and FTP is more than 30 times of C, of HIA and PTP. The proposed CLA can 
support efficient information communication and integration in the system, because it not 
only requires less communication cost, but also tolerates a certain level of sensor and Hnk 
failures. For example. Table 4.3 shows that when Ns is greater than 15. C, of CLA and BP 
is less than 1/4 of C, of HIA and FTP. and Q of CLA and EWGP (wiUi Ni, equal to 1) is 
less than 1/8 of C, of HIA and PTP. 

The best communication protocols depend on the number of sensor nodes in the system, 
sensor constiraints and requirement of response time of the system. Although BP has 
better fault tolerance than EWGP. it requires more communication. For example. Table 
4.3 shows tiiat when Ns is equal tolS (N^ is equal to5), Q of CLA and BP is abnost twice 
Cr of CLA and EWGP. and the difference increases when increases. Therefore. BP is 
suitable for tiie system in which tiie sensors and other IC components are not very robust. 
Compared witi, Uie FTTP with blind base station, the FTTP witii smart base station 
requires less communication time, but it demands the base station to have an additional 
knowledge base tiiat keeps track of the information of tite sensor nodes. Therefore, FITP 
with smart base station is preferred in die on-line control applications, in which response 
time is critical. 
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In order to evaluate the efficiency of micro sensor network architectures and 
conununication protocols, TDE/MEMS has been developed using MPI. TIE/MEMS can 
be expanded by adding libraries for the newly developed network architectures and 
communication protocols. 

Information from different sensw nodes needs to be integrated. In oider to d&bect the 
faulty sensors and minimize their effects on the system, a FTSIA was proposed in this 
. research. The results from simulation experiments show that when f, is large relative to 
Ns, and when/w >0, the FTSIA gives more accurate results than statistical calculations. 
For example. Table 5.1 shows that when/, is equal to 1 and Ns is equal to 8, when is 
equal to 3 and M is equal to 15, and when/, is greater than 0, the mean from FTSIA is 
significantly more accurate than the statistical mean. The mean from FTSIA is also 
affected by/ and N,. For example. Table 5.2 shows that when Ns is equal to 15. the mean 
from FTSIA when / is equal to 1 or 3 is significantly more accurate than that when / is 
equal to 5 and the mean from FTSIA when / is equal to 5 is significantly more accurate 
than that when / is equal to 7. The experiment results indicate that Ns should satisfy 
Equation (6) in order for FTSIA to give more accurate results. 



7.3. Future Research 

Three issues that have not been addressed in this thesis yet important in the design of 
DMFSA/DMFSNarc: 

1. Optimal selection and placemmt of sensor arrays: Optimal selection and 
placement of sensor arrays is not a trivial issue in the flow measurement because 
of several factors, such as the complex flow environment, possible interferences 
of the sensors, etc. Although the optimal placement depends on the flow 
environment, there might be some algorithms (functions) that can be applied in 
planning placement in the general flow envut)nment. 

2. Multi-dimensional FTSIA: A multidimensional sensor is one that returns data in 
D dimensions where D>1. In the flow environment, the velocity can be a 3- 
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dimensional vector. In such case, multi-dimensional fault-tolerant sensor 
integration algorithm needs to be developed. 
3. Further development of TIE/MEMS: MEMS flow sensor physical constraints (as 
shown in Table 5.1) need to be quantified and incorporated into TDE/MEMS to 
evaluate the sensor network architecture and communication protocols. 
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Figure A. 1 Flow Diagram of TIE/MEMS Programming Structure 
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0.000215 




12 


0.00000! 


0.000163 


0.000164 




12 


0.000002 


0.000263 




13 


0.000003 


0.000125 


0.000128 




13 


0.000004 


O.00O234 


A tftwyvi 


14 


0.000002 


0.00013 


0-000132 




14 


0.000002 


O.00O227 




15 


0.000002 


0.000154 


0000156 






0.000002 


0.000246 


A tVWiAS 


16 


0.000002 


0.000146 


0.000148 




16 


0:000002 


0.000258 


0.00026 


17 


0.000002 


0.000123 


0.000125 




17 


O.000002 


0.00026 


0.000262 


18 


0.000002 


O.O0OI3I 


0.000133 




18 


0.000002 


0.00021! 


0.000213 


19 


0.000002 


O.000152 


0.000154 




19 


0.000004 


0.000329 


0.000333 


20 


0.000002 


0.000162 


0.000164 




20 


0.000003 


O.000307 


0.00(»l 


21 


0.000002 


0.000164 


0.000166 




21 


0.000003 


0.000287 


0.00029 


22 


0.000003 


0.000164 


0.000167 




22 


o:ooooo2 


0.00027 


0.000272 


23 


0.000002 


0.000136 


0.000138 




23 


0.000003 


0.000299 


a0003Q2 


24 


0.000002 


0.000126 


0.000128 




24 


0.000002 


0.000245 


0.000247 


25 


O.0000O2 


0,00012 


0.000122 




25 


0.000004 


0.000304 


0.000308 


26 


0.000002 


0-000175 


0.000177 




26 


O.0000O3 


0.000241 


0.000244 


27 


0.000002 


0.000123 


0.000125 




27 


0.000002 


0.000244 


0.000246 


28 


0.000002 


0.000128 


0.00013 




28 


0.000002 


0-000278 


0.00028 


29 


0.000002 


0.000124 


0.000126 




29 


O.O000O2 


0.00029 


0.000292 


30 


0.000003 


0.000162 


0.000165 




30 


O.0O0001 


0.000269 


0.00027 



(Com.) 
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c. //^ =15 


.N^ =sl 




i 


1 =21.^, =1 




Treatment 




r 


Ci 




1 reoimciji 






c, 


1 


0.000002 




0.000420 




1 


0.000004 


0000663 


0.000667 


2 


0 000002 




0.000481 




2 


0.000003 


0 000854 


0.000857 


*/ 


0 000002 




0,000442 




3 


<MXX)004 




0.000642 


A 


0000001 




0.000498 




4 


0.000004 


0 000784 


0.000788 








P.0P0453 




5 


0.000004 


0.0008^ 


OOOC^ 


V 


0000003 


n omvio'> 


0.000495 




6 


O.O00006 


0.000794 


0.000800 


7 


0 00000$ 




0.000489 




7 


0.000004 


0000946 


0;000950 


fi 
O 






0.000472 




8 


0.000003 




0.000772 


9 






0.000464 




9 


O.000006 


0.000751 


0.000757 




0000004 




&j0065O4 




10 


0.000003 


0000688 


0.000691 


11 


& 000003 


n nofu3Q 


0.000442 




11 


0.000004 


0.000734 


0.000738 


12 


00000O2 




0.000454 




12 


0.000004 




0.000843 




0 000003 




0,000512 




13 


0.000003 


0000601 


0.000694 


Id 






0.000460 




14 


0 000003 


0 oooomt 


0.000911 


15 


0 000002 


n nofu*70 


0.000481 




IS 


0.000002 


0000906 


0.000908 




0 000002 




0.000542 




16 


0.000004 


0.000995 


0.000999 


17 


0.000002 


0000177 


0.000474 




17 


0.000004 


0.000734 


0.000738 


18 


0.000002 


0.000730 


0.000732 




18 


0.000004 


0.000674 


0000678 


19 


0.000005 


0.000497 


0.000502 




19 


0.000004 


0.000764 


0.000768 


20 


0.000002 


0.000532 


0,000534 




20 


0.000004 


0.000984 


0.000988 


21 


0.000004 


0.000430 


0.000434 




21 


0.000004 


0.000595 


0,000599 


22 


0.000002 


0.000456 


0.000458 




22 


0.000004 


0.000866 


0.000870 


23 


0.000002 


0.000433 


0.000435 




23 


0.000004 


0.000762 


0.000766 


24 


0.000003 


0.000494 


0.000497 




24 


0.000003 


0.000683 


0.000686 


25 


0.000004 


0.000464 


0.000468 




25 


0.000004 


0.000744 


O.000748 


26 


0.000002 


0.00(M24 


0.000426 




26 


0.000004 


0.000747 


0.000751 


27 


0.000003 


0.000441 


0.000444 




27 


0.000005 


0.000716 


O.00O721 


28 


0.000002 


0.000476 


0.000478 




28 


0.000005 


0.000840 


0.000845 


29 


0.000005 


0.000472 


0.000477 




29 


0.000003 


0.000575 


0.000578 


30 


0.000003 


0.000437 


0.000440 




30 


0.000005 


0.000722 


O.O0O727 



(Cont.) 





e. A^^ = 4 


sW, =2 






f. iV^=8,Af, =2 




Treatment 


Cet 


Cti 






Treatment 


Cct 


Co 


c, 


1 


0.000003 


0.000362 






1 


0.000006 


0.000577 


0.000583 


2 


0.000005 


0.000339 


U.UUUJ44 




2 


0.000004 


0.000652 


0.000656 


3 


0.0000O4 


O.00O32O 


V.UUU344 




3 


0.000004 


0.000581 


0.00(»85 


4 


0.000004 


0.000358 


U.UUU362 




4 


0.000005 


0.000651 


0.000656 


5 


0.000004 


0.000342 


A AAA'S 




5 


0.000006 


0.000592 


0.000598 


6 


0.000002 


O.0OO39O 


A AAA^OI 




6 


aoQooos 


0.000513 


0.000518 


7 


0.000003 


0.000369 


A AAAl-7-9 




7 


0.000005 


0.000595 


0.000600 


8 


0.000003 


0.000376 


AnnAt70 




8 


0.000005 


0.000583 


0.000588 


9 


0.000004 


0.000370 


A AAAt^it 




9 


0.000003 


0.000557 


0.000560 


10 


0.000002 


0.000344 


A Afy|T4iC 
U.UUIU40 




10 


0.000004 


0.000SS6 


0.000590 


11 


0.000003 


0.000346 






11 


0.000004 


0.000584 


0.000588 


12 


0.000004 


0.000346 


A AAA^</% 

U.UUU350 




12 


0.000003 


O.00O583 


0.000586 


13 


0.000002 


0.000313 






13 


0.000005 


0.000504 


0.000509 


14 


0.000002 


0.000308 


A AAA')* A 

0.000310 




14 


0.000004 


0.000541 


O.O0O545 


15 


0.000003 


0.000319 


A AAA'^'V) 




15 


0.000005 


0.00O56I 


0.000566 


16 


0.000004 


0-000384 


0.000388 




16 


0.000004 


0000512 


0.000516 


17 


0.000004 


0.000338 


0.000342 




17 


0.000006 


0.000588 


0.000594 


18 


0.000005 




0.000296 




18 


0.000003 


0.000623 


0.000626 


19 


0.000003 


0.000382 


0.000385 




19 


0.000003 


0.000595 


0.000598 


20 


0.000002 


0.000353 


0.000355 




20 


0.000004 


0.000602 


0.000606 


21 


0.000005 


0.000301 


0.000306 




21 


0.000003 


0.000623 


0.000626 


22 


0.000002 


0.000374 


0.000376 




22 


0.000006 


0.000550 


0.000556 


23 


0.000004 


0.000330 


0.000334 




23 


0.000003 


0.000656 


0.000659 


24 


0.000002 


0.000339 


0.000341 




24 


0.000004 


0.000593 


0.0)0597 


25 


0.000003 


0.000373 


0-000376 




25 


O.OO00O4 


O.O0O571 


0.000575 


26 


O.O00OD2 


0.000310 


0.000312 




26 


0.000004 


0.000541 


0.000545 


27 


0.000003 


0.000298 


0.000301 




27 


0.000003 


0.000520 


0XX)Q523 


28 


0.000003 


0.000371 


0.000374 




28 


0.000004 


0.000630 


0.000634 


29 


0-000002 


0.000374 


0.000376 




29 


0.000003 


0.000568 


0.000571 


30 i 


0.000002 


0.000378 


0.000380 




30 


0.000004 


0.000625 


0.000629 



63 



(Cont.) 



g- A^,=lS,Ar, =2 h. = 21.^, = 2 



Treatment 




Ctt 


Ci 




Treatment 




Cct 


Ct 


1 


0.000005 


0.00I05& 


0.0O1063 




1 


0.000008 


0.001726 


0.001734 


2 


0.000006 


0.001044 


0.001050 




2 


0.000007 


0.001637 


OJ001644 


3 


0.000003 


0.001193 


0.001196 




3 


O.0OO0O6 


0.001834 


0.001840 


4 


0.000006 


0.001255 


0.001261 




4 


0.000008 


0.001882 


0.001890 


5 


O.0O0O08 


0.00JI74 


0.001182 




5 


0.000006 


0.001507 


0.001513 


6 


0.000006 


0.001050 


0.001056 




6 


O.0000O6 


0.001637 


0.001643 


7 


0.000006 


0.001164 


O.OOII70 




7 


0.000007 


0.001719 


0.001726 


8 


0.000006 


0.001226 


0.001232 




8 


0^)00006 


0.001458 


0.001464 


9 


0.000004 


0.001230 


0.001234 




9 


0.000006 


0,001649 


0.001655 


1& 


O.O0D0O7 


0.001298 


0001305 




10 


^.000609 


0.001944 


0.001953 


11 


0.000004 


0.001044 


0.001048 




11 


0.000006 


0.001906 


0.001912 


12 


0.000006 


0.001302 


0.001308 




12 


0.000006 


0.001725 


0.001731 


13 


0.000005 


0.000966 


0.000971 




13 


0.000007 


0.001499 


0.001506 


14 


0.000010 


0.0012S4 


0.001294 




14 


0.000008 


0.001798 


0.001806 


15 


0.000007 


0.001142 


0.001149 




13 


0.000006 


0.001685 


0.001691 


16 


0.000006 


0.001214 


0.001220 




16 


0.000008 


a.001891 


0.001899 


17 


0.000009 


0.00 109O 


0.001099 




17 


0.000006 


0.001625 


0.001631 


18 


0-000004 


0.001038 


0.001042 




18 


0.000006 


0.001635 


0.001641 


19 


0.000005 


0.001236 


0.001241 




19 


0.000005 


0.001706 


0.001711 


20 


0.000006 


0.001113 


0.00II19 




20 


0.000010 


0.(X)1793 


0.001803 


21 


0.000006 


O.00U54 


0.001160 




21 


0.000008 


0.001799 


0.001807 


22 


0.000006 


0.001182 


0.001188 




22 


0.000009 


0.001S36 


0.001545 


23 


0.000005 


0.000948 


0.000953 




23 


0.000007 


0.001734 


0.001741 


24 


0.000006 


0.001212 


0.001218 




24 


0.000006 


0.001703 


0,001709 


25 


O.0000O4 


0.001021 


0.0O1025 




25 


O.OOOOlO 


0.001579 


0.001589 


26 


0.000006 


0.001194 


0.001200 




26 


0.000010 


0.001566 


0.001576 


27 


0.000006 


0-0OJ282 


O.0O1288 




27 


0.000009 


0.001880 


0.001889 


28 


0.000009 


0.001354 


0.001363 




28 


0.000009 


0^1894 


0.001903 


29 


0.000006 


0.000778 


0.000784 




29 


0.000009 


0.001615 


0.001624 


30 


0.000003 


0.001149 


0.001152 




30 


0.000008 


0.001622 


0.001630 



(Cont.) 











iV,«8,iV, =3 




Treatment 


C— 








Treatment 






Ct 


1 
1 






0.000540 




1 


A AAAAAK. 


0.000820 


0.000826 


2 






0.000560 






A AAAnA< 


0.000836 


0.000842 




0 oooocvi 




0.Q0O572 








0.000^2 


0.000728 


4 


0000005 




0.000576 




A 


A AAAAAA 
U.UUUUUO 


0.000822 


0.000828 




0000004. 




0.000521 






A nAAAAK 


A AAAO a£ 

0.000846 


0.000852 




0 000005 




0.000602 




O 


A AAAAAiC 
U.UUUUUO 


0.000837 


OX)00843 


7 


0.000004 


V.VVUwU't 


0.000608 




7 


A AAAAAil 
V.UUUUU4 


0.000884 


0.000888 


o 


0 000006 




0.000633 




o 
o 


A AAnAA< 
U.UUUUUO 


0.000867 


O.b00872 


9 


OOOOOQS 




0.000590 




o 


A AAAAA< 


0.000884 


0.000890 




0-000004 




0.000526 




in 


A AAAAAiC 


0.000824 


0.000830 


1 1 
A 1 






0000489 




i 1 


A AAAAnr 


0.00081 1 


0,000818 




V.W/V/VA/J 




0.000559 




iZ 


A AIWV^ 
U.UUUUUO 


0.000847 


0.000853 




0 000005 




0.000514 






U.UUUUUO 


0.000892 


0000898 


14 




n Am<^D 


0.000551 




1** 


A AAAAA< 

U.UUUUU5 


0.000814 


0.000819 


IS 




U.UUUOjy 


0.000622 






A AAAAA< 


V.000899 


0.000904 


16 


OOOOOOd 




0.000557 




16 


0.000006 


0.000812 


0.000818 


17 


0000005 


A fiAACOT 
If. UUU^Z/ 


0.000532 




17 


0.000006 


0.000793 


0.000799 


18 


0.000005 


0.000586 


0.000591 




18 


0.000005 


0.000875 


0.000880 


19 


O.O0OOOS 


0.000549 


0.000554 




19 


0.000005 


0.000786 


0.000791 


20 


0.000003 


0.000539 


0.000542 




20 


0.000006 


0.000860 


04X)0866 


21 


0.000005 


0.000476 


0.000481 




21 


0.000005 


0.000813 


0.000818 


22 


0.000004 


0.000602 


0.000606 




22 


0.000005 


0.000721 


0.000726 


23 


0.000004 


0.000578 


0.000582 




23 


0.000006 


0.000830 


0.000836 


24 


0-000005 


0.000549 


0.000554 




24 


0.000006 


O.00O86O 


0.000866 


25 


0.000006 


0.000567 


0.000573 




25 


0.000005 


0.000885 


0.000890 


26 


0.000006 


0.000574 


0.000580 




26 


0.000006 


0.000858 


0.000864 


27 


0.000004 


0.000551 


0.0005S5 




27 


0000006 


0.000822 


O.0OO828' 


28 


0.000005 


0,000575 


0.000580 




28 


0.000006 


0.000788 


0.000794 


29 


0.000006 


0.000529 


0.000535 




29 


0.000006 


0.000776 


0.000782 


30 


0.000005 


0.000553 


0.000558 




30 


0.000004 


0.000739 


0.000743 



(Cont.) 



65 





k. AT^ = 1 


5,iV, =3 




I. =21,7V, =3 




Treatment 






C, 




Treatment 




Ctt 


Ct 


1 


0.000006 


0 00174^ 


0,001751 




1 


0.000008 


0X)02769 


0.002777 


2 


0.000008 


0 001 696 


0.001704 




2 


0.000010 


0.002433 


0.002443 


3 


0.000008 


0.00)752 


0.001760 




3 


aoooou 


0.002743 


0X)Q27S4 


4 


0.000007 


0 OOII^A 


0.001647 




4 


0.000009 


0.002646 


0.002655 


5 


0000008 


v.vvto/j 


O.Oa}88l 




5 


0.000011 


0.002777 


0.002788 


< 


oooooos 




0.001781 




6 


0.000010 


0.002700 


0.0(tt710 


7 


0 000006 




0.001739 




7 


0.000009 


0.002552 


0.002561 


Q 
O 


0000000 


U.UUloil 


0.001820 




8 


0.000009 


O.0O2452 


0.002461 


o 


0 OOOOOR 


U.UUi tot 


0,001795 




9 


o.oooon 


0.002550 


a00256] 


Xv 


OOOOOOR 




0:0ftl759 




10 


0.000010 


0.002813 


&002823 


1 1 
1 1 


OOOOOOf? 


n AAi7ni 


0,001710 




11 


0.000013 


0.002616 


0.002629 




0 000008 


n nfniiR 

\J.\AjX / I o 


O.0O1726 




12 


0.000010 


0.002627 


0.002637 




0 00(K)07 


U.W 1 o*y 


0.001636 




13 


0.000006 


0.002510 


0.002516 


14 


0.000009 


A 0A1 ^Q4( 


0.001605 




14 


0.000009 


0.002825 


O.0O2834 


15 


0 OOOOOfi 




0.001595 




15 


0.000013 


0.002550 


O.0O2563 


16 


0 000007 


A AAISIO 


0.001846 




16 


aoooou 


O.002544 


0.002555 


17 


0000009 


A AA1 Kd^ 


0.001552 




17 


0,000009 


0.002608 


0.002617 


18 


o.ooooos 


0.001574 


0.001579 




18 


0.000010 


0.002872 




19 


0.000010 


0.001559 


0,001569 




19 


0.000010 


0.002635 


0.002643 


20 


0.000008 


0.001737 


0.001745 




20 


0.000010 


0.002742 


O.0O2752 


21 


0.000010 


0.001685 


0.001695 




21 


0.000014 


0.002538 


0,002552 


22 


0.000010 


0.001602 


0.001612 




22 


0.000014 


0.002484 


O.0Q2498 


23 


0.000007 


0.001571 


0.001578 




23 


0.000011 


0.002486 


O.0O2497 


24 


0.000007 


0.001922 


0.001929 




24 


0.000013 


0.002365 


0.002378 


25 


0.000006 


0.001832 


0.001838 




25 


0.000009 


0.002533 


0.002542 


26 


0.00Ci006 


0.001782 


0.001788 




26 


0.000011 


0.002365 


0,002376 


27 


0.000007 


0.001789 


0.001796 




27 


0.000009 


0.002433 


O.0O2442 


28 


0.000010 


0.001740 


0.001750 




28 


0.000010 


0.002283 


O.0O2293 


29 


0.000006 


0.001862 


0.001868 




29 


0.000013 


0.002587 


O.0O260O 


30 


0.000010 


0.001813 


0.001823 




30 


0.000024 


0.002570 


0.002594 



11^ 



Table A3 Communication Time (Seconds) of the FTTP with Smart Base Station vs. No, 
of Sensor Nodes in the Cluster and No. of Failed Links 





a. = 4, 


=1 








Treatment 




Oft 


Ct 




Treatment 




Cn 




I 


0.000003 


0.000021 


0.000024 




I 


0.000007 


0.000029 


U.O(XX}36 


2 


O.0G00O2 


0.000016 


0.000018 




2 


0.000006 


0.000029 


0.000035 


3 


0.000006 


0.000021 


0.000027 




3 


0.000007 


0.000027 


0.000034 


4 


0.000003 


0.000019 


0.000022 




4 


0.000009 


0.000029 


0.000038 


5 


0.000005 


0.000025 


0.000030 




5 


0.000007 


0.000031 


A AAAA<%B 

0.000038 


6 


0.000005 


0.000021 


0.000026 




6 


0.000006 


0.000026 


0.000032 


7 


O.0000O5 


0.000017 


0.000022 




7 


0.000005 


0.000(Q7 


0.000032 


8 


0.000(»)3 


0:000021 


0.000024 




8 


0.000005 


0000031 


U.0U00J6 


9 


0.000006 


O.QOQ022 


O.00OQ28 




9 


0.000006 


0.000028 


o.puoa34 


10 


0.000004 


0.000038 


0.000042 




10 


0.000D07 


0.000032 


U.U00U39 


11 


0:000002 


0000024 


0.000026 




11 


0.000006 


0.000029 


□.000035 


12 


0.000006 


0.000030 


0.000036 




12 


0.000006 


0.000029 


0.000035 


13 


0.000004 


0.000026 


0.000030 




13 


0.000006 


0.000026 


0.000032 


14 


0.000003 


0.000023 


0.000026 




14 


0.000006 


0.000033 


0.000039 


15 


0.000005 


0.000033 


O.O0OO38 




15 


0.000006 


0.000026 


0.000032 


16 


0.000006 


0.000021 


0.000027 




16 


0.000003 


0.000029 


A AAAA'tt 

0.000032 


17 


0.000006 


0.000020 


0.000026 




17 


0.000006 


0.000026 


0000032 


18 


0.000006 


0.000016 


0.000022 




18 


0.000005 


0.000027 


0.000032 


19 


0.000002 


0.000024 


0.000026 




19 


0.000004 


0.000031 


0.000035 


20 


0.000004 


0.000021 


0.000025 




20 


0.000007 


0.000028 


0.000035 


21 


0.000004 


0.000021 


0.000025 




21 


0.000005 


0.00002S 


0.000030 


22 


0.000002 


0.000020 


0.000022 




22 


0.000005 


0,000029 


0.000034 


23 


O.OOOOC^ 


0.000016 


0.000022 




23 


0.000006 


0000028 


0.000034 


24 


0.000003 


0.000021 


0.000024 




24 


0.000007 


0.000029 


0.000036 


25 


0.000004 


0.000023 


0.000027 




25 


0.000007 


0.000035 


0.000042 


26 


0.000004 


0.000022 


0.000026 




26 


0000005 


0.000027 


0.000032 


27 


0.000005 


0.000022 


0.000027 




27 


0.000004 


0.000030 


0.000034 


28 


0.O00O06 


0000027 


0.000033 




28 


0000004 


0.000030 


0.000034 


29 


0.000007 


0000022 


O0G0029 




29 


0000003 


0.000029 


0.000032 


30 


0-000003 


0.000020 


0.000023 




30 


O.00O0O4 


0.000028 


0.000032 



HQ) 



67 



(Cont.) 



Treatment 




Cn 


Ct 




Treatment 






C| 


1 


0.000006 


0.000040 


0.000046 




i 


U.0UU0U9 


0,000032 


0.000041 


2 


0.000d09 


04)00043 


0^000052 




Z 


U.UUU0U7 


0.000048 


0J00O0S5 


3 


0.000008 


0.000045 


0.000053 




t 


U.UIIQ0U7 


O.O0Q035 


0.000042 


4 


0.000006 


0.000031 


0.000037 




A 
** 


v.OUOUII/ 


0.000047 


O.O0OQS4 


5 


0000006 


0.000030 


0.000036 






U.UUUOOo 


0.000067 


aoooo^ 


6 


O.Q00008 


O.Q00045 


0.000053 




0 


O.U0II0Q8 


O.O00029 


0.000037 


7 


0.000009 


0.000037 


O.000046 




/ 


A AAAAAO 


0.000021 


0.000029 


8 


O.OOOOIO 


0.000033 


0.000043 




Q 

o 


A Afwmo 


0.000042 


0.000051 


9 


0.000008 


0:000026 


0.000034 




0 




0.000034 


0.000042 


10 


0.000007 


0.000024 


O.O00031 








0:000038 


0.000047 


11 


0.000005 


0.000021 


0.000026 




11 

Ml 


A AfWtfkin 
l/.UUUUlU 


0.000032 


0.000042 


12 


0.000006 


0.000024 


0.000030 








0.000035 


0.000043 


13 


0.000009 


O.O00O38 


0.000047 




13 


0.000009 


0000018 




14 


0.000004 


O.000042 


0.000046 




14 


0.000009 


0.000033 


0.000042 


15 


0.000009 


0.000041 


0:000(»0 




15 


0.000005 


0.000056 


0.000061 


16 


0.000007 


0.000048 


0:000055 




16 


O.O00008 


0.000033 


o.om)4i 


17 


0,000010 


0.000024 


0.000034 




17 


aooobo9 


0.000030 


0-000039 


18 


0.000007 


0.000042 


0.000049 




18 


0.000010 


0.000035 


0.000045 


19 


0,000005 


0.000042 


0.000047 




19 


0.000006 


0.000047 


0.000053 


20 


0.000006 


0.000029 


0.000035 




20 


0.000006 


0.000039 


0.000045 


21 


0.000006 


0.000046 


0.000052 




21 


0.000009 


O.0O003O 


0.000039 


22 


O.00O0O7 


0.000031 


0.000038 




22 


0.000009 


0.000034 


0.000043 


23 


0.000006 


0.000038 


0.000044 




23 


O.OOOOIO 


0.000032 




24 


0.000006 


0.000015 


0.000021 




24 


0-000006 


0.000036 


0.000042 


25 


0.000006 


0.000027 


0.000033 




25 


0.000009 


0.000036 


0.000045 


26 


0,000010 


0,000028 


0.000038 




26 


0.000007 


0.000041 


0.000048 


27 


0.000007 


0.000024 


0.000031 




27 


0.000010 


0.000030 


0.000040 


28 


0,000007 


0,000034 


0.000041 




28 


0.000009 


0.000049 


0.000058 


29 


0.000008 


0.000017 


0.000025 




29 


0.000006 


0.000039 


0.000045 


30 


0.000006 


0.000026 


0.000032 




30 


0.000009 


0.000040 


0.000049 



117 



68 



(Com.) 



e. N, =4,N, =2 f. AT, =8.A^, =2 



Treatment 


Pel 


Cn 






& 1 VAUIlCIii 


Co 


Cn 


Ct 




0.000012 


0 000045 






1 


0.000011 


0.000063 


0.000074 


2 


0.000008 


OOOOQSI 






2 


0.000010 


&.000056 


0.000066 




0000007 

W» Www 9 








3 


o.oooou 


0.000054 


0.000065 


A 


0000008 

w- VWWQ 


0 0000^7 






4 


0.000009 


0.000057 


0.000066 


5 


0.000007 


6000(K4. 






5 


0.000010 


0.000042 


a000Q52 


A 


0000006 








6 


0.000009 


0.000054 


0.000063 


7 


0 000006 


ooo(vn7 


AIVWU.1 




•* 

7 


0.000010 


0.000056 


0.000066 


g 


0.000006 


0.000066 






8 


0.00001 1 


0.000065 


O.00OQ76 


9 


0.000007 


0.000040 






9 


0.000009 


0.000042 


0.000051 


10 


0.000007 


0.000041 






10 


o.oooou 


0.000051 


0.000062 


1 1 


0.000006 


0000057 






11 


0.000009 


0.000054 


0.000063 


12 


0.000005 








12 


0.000009 


0.000045 


O.000054 


13 


0.000006 


0000054 






13 


0.000013 


0.000070 


0.000083 


14 


0.000007 


O.00O038 


0.000045 




14 


0.000010 


0.000049 


O.00OOS9 


15 


0.000006 


0.000043 


0.000049 




15 


O.OOQ007 


0.000043 


0.000050 


16 


0,000009 


0.000038 


0.000047 




16 


0.000009 


0.000046 


0.000055 


17 


0.000009 


0.000043 


0.00(»)52 




17 


O.O000Q9 


0.000054 


a000063 


18 


0.000007 


0.000052 


O.O00QS9 




18 


0.000009 


0.000051 


0.000060 


19 


0.000008 


0.000041 


0.000049 




19 


0.000010 


0.000054 


0.000064 


20 


0.000009 


0.000049 


0.000058 




20 


0.000009 


0.000071 


0.000080 


21 


0.000007 


0.000045 


0.000052 




21 


0.000009 


0.000076 


0.000085 


22 


0.000008 


aoooos6 


0.000064 




22 


0.000011 


0.000049 


0.000060 


23 


0.000010 


0.000051 


0.000061 




23 


0.000011 


0.000045 


0.000056 


24 


O.OCK)008 


0.000044 


0.000052 




24 


o.oooou 


0.000058 


0.000069 


25 


0.000008 


0.000044 


0.000052' 




25 


o.oooou 


0.000051 


0.000062 


26 


0.000008 


0.000038 


0.000046 




26 


0.000009 


0.000061 


0.000070 


27 


0.000008 


0.000048 


0.000056 




27 


0.000010 


0.000054 


0.000064 


28 


0.000008 


0.000044 


0.000052 




28 


0.000010 


0.000056 


0.000066 


29 


0.000009 


0.000036 


0.000045 




29 


0.000009 


0.000056 


0.000065 


30 


0.000008 


0.000054 


0.000062 




30 


0.000010 


0.000055 


0.000065 



IIS 



69 

(Cont.) 



Treatment 


c« 


Ci, 


Ct 




Treatpoent 






Ct 


1 


0,000009 


0.000060 


0.000069 




1 


0.000019 


0,000053 


0.000072 


2 


0.000017 


0.000043 


0,000065 




2 


0.000020 


0.000054 


0,000074 


3 


0.000012 


0.00006O 


04)00072 




3 


0.000026 


0.000066 


0.000092 


4 


0.000010 


0.000067 


a000077 




4 


0.000017 


0.000071 


0,000(^8 


5 


0.000018 


O.000OS3 


0.000071 




5 




O.OOQ070 


0000090 


6 


0.000017 


0.000072 


0.000089 




6 


0.000014 


0.000063 


0.000077 


7 


0.000022 


0.000062 


0.000084 




7 


0.000(ti!l 


0,000046 


0.000067 


8 


0.000009 


0.000051 


0.000060 




8 


0.000016 


0.000076 


O.000O92 


9 


0-00001 1 


0-000069 


0.000080 




9 


0.000021 


0.000053 


0.000074 


10 


0.000011 


0.000052 


0.000063 




10 


0.000010 


0.000055 


0.000065 


11 


0.000010 


0-000051 


0.000061 




11 


0.000013 


0.000083 


0.000096 


12 


0.000017 


0.000048 


0.000065 




12 


0.000023 


0.000050 


0-000073 


13 


0.000012 


0.000068 


0.000080 




13 


0.000019 


0.000066 


0.000085 


14 


0,000008 


0.000039 


0.000047 




14 


0.000018 


0.000040 


0.000058 


15 


0.000012 


0.000048 


0.000060 




15 


0.000014 


0.000051 


0-000065 


16 


0.000018 


0.000067 


0.000085 




16 


0.000011 


0.000069 


0.000080 


17 


0,000012 


0,000052 


0.000064 




17 


O.OO0O2O 


0.000058 


0.000078 


18 


0.000020 


0.000049 


0.000069 




18 


0.000018 


O.O0OO45 


O.O0OO63 


19 


O.OOOOH 


0.000071 


0-000082 




19 


0.000019 


0.000061 


0.000080 


20 


0.000022 


0.000051 


0.000073 




20 


0.000014 


0.000073 


0.000087 


21 


0.000018 


0.000064 


0.000082 




21 


0.000020 


0.000048 


0,000068 


22 


0-000013 


0.000051 


0.000064 




22 


0.000020 


O.O0QO78 


0.000098 


23 


o.oooon 


0.000059 


0.000072 




23 


0.000022 


0.000075 


0.000097 


24 


0.000010 


0.000060 


0.000070 




24 


0.000020 


0.000063 


0.000083 


25 


0.000018 


0.000056 


0.000074 




25 


0.000020 


0.000047 


0.000067 


26 


O.0OOOU 


0.000055 


0.000069 




26 


0.000018 


0.000054 


0.000072 


27 


0.000010 


0.000072 


0.000082 




27 


0.000025 


0.000062 


O.0O0O87 


28 


0.000016 


0.000066 


O.0O0O82 




28 


0.000019 


O.O0OO69 


0.000088 


29 


0,000017 


0.000060 


0.000077 




29 


0.000018 


0.000064 


0.000082 


30 


0.000010 


0.000074 


0.000084 




30 


OX)00020 


0.000057 


0.000077 



II? 



(Cont.) 



70 





i. iV, = 4 










TreatB^ent 


c« 




c, 




Treatment 








1 


0.000011 


0.000051 


0.000062 




1 


0.000011 


0.OQOO44 


0.000055 


? 


0.000010 


0.000056 


0.000066 




2 


0J0O0O\6 


boooOfSff 


0.000072 


3 


0.0000/1 


0.000070 


0.000081 




3 


0.000012 


OO0QO5S 


0.000067 


4 


0.000007 


0.000070 


0.000077 




4 


0.000012 




0.000064 


5 


0.O0OOU 


0.000057 


0.000068 




5 


0.000010 


0.000048 


0.000058 


6 


0.000014 


0.000078 


0-000092 




6 


0.000014 


0000058 


O.000072 


7 


0.000010 


0.000081 


0.d0009l 




7 


0.000014 


0.000074 


0:000088 


8 


0.000011 


0.000061 


0.000072 




g 


0.000014 


0000066 


0.000080 


9 


D.oooon 


0.000054 


0.000065 




9 


O.0O0OU 


0000045 


0.000056 


10 


0.000009 


0.000056 


0.000065 




10 


0.000013 


0.000048 


0.000061 


11 


0.000011 


0.000055 


0.000066 




11 


0.000012 


0000056 


0.000068 


12 


0.000009 


0.000059 


0.000068 




12 


0.000014 


0.000053 


O.(K)0067 


13 


0.00001 1 


0.000059 


0.000070 




13 


0.000014 


0.000038 


0.000052 


14 


0.000009 


O.000O74 


0.000083 




14 


0.000014 


0.000062 


0.000076 


15 


0.000013 


0000053 


0.000066 




15 


0.000014 


0.000065 


0.000079 


16 


0.000010 


0.000045 


0.000055 




16 


0.000014 


0.000079 


0.000093 


17 


0.000014 


0.000058 


0.000072 




17 


OJ00O0i3 


0.000079 


0.000092 


18 


0.006010 


0:000054 


0000064 




18 


0.000014 


0.000064 


0.000078 


19 


0.000009 


0.000084 


0.000093 




19 


0.000016 


0.000064 


0.000080 


20 


0.000012 


0.000071 


O.OO0O83 




20 


0.000017 


0.000050 


0.000067 


21 


0.000010 


0.000050 


0.000060 




21 


0-000016 


0.000052 


0.000068 


22 


0.000009 


0.000080 


0.000089 




22 


o.oooon 


0.000059 


0.000076 


23 


0.000011 


0.000053 


0.000064 




23 


0.000014 


0.000046 


0.000060 


24 


0.000012 


0.000072 


0.000084 




24 


o.oooon 


0.000065 


0.000076 


25 


o.oooon 


0.000053 


0.000064 




25 


0.000010 


0.000052 


0.000062 


26 


0.000012 


0.000049 


0.000061 




26 


0.000016 


0.000052 


0.000068 


27 


0.000011 


0-000059 


0.000070 




27 


0.000016 


0.000073 


0.000089 


28 


0.000009 


0.000068 


0.000077 




28 


0.000014 


0.000052 


0.000066 


29 


0.000010 


0.000055 


0.000065 




29 


O.000OM 


O.0000S8 


0.000072 


30 


Qjmon 


0.000061 


O.000O73 


i 


30 


0.000015 


oimns 


0.000090 



(Com.) 



k. =15, N, =3 



Treatment 


c« 


Ct, 


Ct 


1 


0.000018 


0.000085 


0.000103 


2 


0.000013 


0:000060 


0.000073 


3 


O.00O012 


O.O0Q076 


0.000088 


4 


0.000018 


0.000066 


0.000084 


5 


0.000012 


0.000061 


0.000073 


6 


O.O0OOU 


0.000082 


0.000096 


1 


0.060023 


0.000076 


0000099 


8 


D.000019 


O.OO0060 


0.000079 


9 


O.000OI4 


0.000077 


0.000091 


10 


0.000017 


0.0(KX)58 


0.000075 


11 


0.000016 


0.000039 


0.000055 


12 


0.000018 


O.OO0078 


0.Q00096 


13 


ojoooon 


0.000066 


0.000079 


14 


0.000026 


0.000081 


0.000107 


15 


0.000022 


0.000062 


0.000(»4 


16 


0.000011 


0.000068 


0.000079 


17 


0.000614 


0.000062 


0.000076 


18 


0.000013 


0.000067 


0.000080 


19 


0.000014 


0.000048 


0.000062 


20 


0.000022 


0.000077 


0.000099 


21 


0.000014 


0.000083 


O.0O0097 


22 


0.000014 


0.000055 


0.000069 


23 


0.000021 


0.000080 


0,000101 


24 


0.000025 


O.O00O50 


0.000075 


25 


O.OOOOIS 


0.000060 


0.000075 


26 


0.000016 


0.000067 


0.000083 


27 


0.000018 


0.000049 


0.000067 


28 


0.000026 


0.000077 


0.000103 


29 


0.000014 


0.000084 


0.000098 


30 


0.000016 


0.000059 


0.000075 



1. N, =21,Ar, =3 



Treatment 




Ca 


Ct 


1 


0.000018 


0.000083 


0.000101 


2 


0.000022 


0.000069 


O.0O0O91 


3 


0.000018 


0.000073 


0.000091 


4 


0.000622 


0.000077 


0.000099 


5 


0.000023 


0.000083 


0.000106 


6 


0.000020 


0.G0OO51 


0.000071 


7 


aO00022 


a000092 


0:0001 14 


8 


0.000026 


0.000039 


O.0O0065 


9 


0.000019 


0.000147 


0.000166 


10 


0.000030 


0.000070 


0.000100 


11 


0.000023 


0.000085 


0.000108 


12 


0.000021 


0.000044 


0.000065 


13 


aoooon 


0.000076 


0.000093 


14 


0.000022 


0.000085 


0.000107 


15 


0.000026 


0.000055 


O.0O0081 


16 


0.000026 


0.000063 


0.000089 


17 


0.000019 


0.000045 


0.000064 


18 


0.000022 


O.O00056 


0.000078 


19 


O.000028 


0.000063 


0.000091 


20 


0.000026 


0.000067 


0.000093 


21 


0.000036 


O.00006O 


0.000096 


22 


0.000025 


O.0OO068 


0.000093 


23 


0.000024 


0.000076 


oxnoloo 


24 


0.000024 


0.000054 


0.000078 


25 


0.000020 


0.000063 


0.000083 


26 


0.000028 


O.00O049 


o.oooon 


27 


0.000022 


0.000047 


a000069 


28 


0.000029 


0.000056 


0.000085 


29 


0.000022 


0.000051 


0.000073 


30 


0.000024 


O.00007O 


0.000094 




APPENDIX B Examples of Execution of FTSIA 



Trae Value 



I3 

I5 

I7 
Is 



I I 

I! 

I I 



I I 
I I 
I I 
I I 
I I 
I f 



Steps: 

Step L Find out distinct S-icterval 
frcwn abstract sensors Ii to Ig AB 
Step 2. Since there is no other abstract 
sensors, AB is the final output interval 



Step 1 Distinh S-loterva) *A 
Step 2 Final output interval 

Figure B. 1 Example of Execution of FTSIA When iV, = 8, / = 0 



True Value 



U 
Is 

Is 



Stepl 
Step 2 
Steps 



Step 5 Final output interval 




DistiscrT-inUr^ 

IMstfnctT-interval 

Distinct 7*intervaJis 
Distinct 8-interv9l 



Steps: 

Step I. Find out the distinct 7-Interval AB 
from first seven abstract sensors AB 
Step 2. Ignore Ij because its higher boui^ is 
the smallest among the first seven abstract 
sensors and find out the distinct 7*interval 
from abstract sensors t to 1% CD 
Step 3. Combine distinct intervals and 
classify different levels of distinct intervals 
(AC and BD are distinct 7-interval$; CB is 
a distinct 8-interval} 

Step 4. Compare the distinct intervals with 
all the abstract sensors and find out the 
possibly faulty sensors: if AC contains the 
true value, is tamely faulty; if CD 
contains the true value, all abstract sensors 
are correct; if BC contains the true value. 1| 
is tamely faulty. Therefore, I, and I, are 
possibly tamely faulty. 
Step 5. Ignore readings from all the 
possibly lamely faulty sensors, and find out 
final output interval fix)m the rest correct 
sensors. 



Figure B.2 Example of Execution of FTSIA When A^, = 8./, = 1 



Trae Value 



I| 

I2 
I3 
U 
Is 
h 
h 
h 



I 



Step2 Dbtfnti^in&rratl 
Step 5 Final eatpiit interval 



u 

ll I 
It I 

!|i 

t 



I I 



Steps: 

Step 1. There is no distinct 6-rnterval 
from first six abstract sensors 
Stq> 2. Ignore Ii because its hi^ier 
bound is die smallest among the first 
six abstract sensors and find out the 
distinct 6-intcrval from abstract 
sensonl2tol7 AB 
Step 3. Ignore reading from Ij, and 
there is no distinct 6-interval from 
abstract I) sensors to Ig 
St^ 4. Compare the distinct intervals 
with all the abstract sensors and find 
out the possibly foulty sensm. Since 
AB is the only distinct 64^tervai, I| 
and Igimist be faulty. 
Step S. AB is the final output interval. 



Figure B3 Example of Execution of FTSIA When A^, = 8, / = 2 



True Value 



Iz 
h 
U 
h 
h 
I7 
Is 



ir 

ii 

1 1 



Step 1 ]Mstiiic«74iileHal * 



Step 5 Ffnal outpat interval ^ 



B 



Steps: 
N,-f=7 

Step 1. End out the distinct 7-intcrval 
from abstract sensors Ii to I7 AB 
Step 2. Ignore It because its higher 
bound is the smallest among the first 
seven abstract sensors and there is no 
distinct 7-interval from abstract 
sensors I2 to Ig 

Step 3. Compare the distinct intervals 
widi an the abstract sensors and find 
out the possibly faulty sensors. Since 
AB is only distimrt 7-interval, Ig must 
be faulty 

Step 4. AB is the final output interval. 



Figure B.4 Example of Execution of FTSIA When N^^S.f^^l 
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True V^ie 



Ii 

h 
I4 
I5 
h 
h 
h 

Stepl 
Step! 
Steps 



IM$t$nc|!6-inier/a! *Ar 



DisliiiclMiterval JC 
flnaicm^t interval Lm 



Steps: 

Step 1. Fmd out the distinct 6-interval 
frotn abstract sensors Ii to I« AB 
Step 2. Igm>re Ig because its higher 
bound is the smallest among the first six 
abstract sonsons and find the distinct 6- 
interval fiom abstract sensors latol? 
CD 

Step 3. Combine distinct intervals and 
classify different levels of distiiKrt 
intervals (bodi AC and BD are distinct 6- 
intervals) 

Step 4. Compare the distinct intervals 
with ail the abstract sensors and find out 
the possibly faulty srasors: if AB 
contains ^ true value, I7 to Ig are faulty; 
if CD contains the true value» I] to Ig are 
faulty. Therefore, I| and I7 are possibly 
tan»ty faulty, and Ig must be feuhy. 
Step 5. Ignore readings fiom possibly 
faulty and faulty sensors, find out the 
final output interval from the rest correct 
sensors 



Figure B^5 Example of Execution of FTSIA When N, = 8. /, = 1, = 1 



i 



75 



-APPENDIX C Experiment Results of FTSIA and Statistical Calculations 

(x*=0.5, e=±0.005) 

In the statistical calculations, the formula used for calculating mean is -j^ (2^ r,. ) and the 



formula used for calculating standard deviation (Std.) is V — ^'^ I' ' — 



Table C.l Experiment Results of FTSIA and Statistical Calculations (Case 1) 



a. AT, =8,/ =0 



TFeatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


I 


0.49796252 


0J0599591 


0.5019792 


0.0019792 


03020511 


0.0020511 


0.0007325 


2 


0.49979105 


0.5O5089S2 


0.5024404 


0.0024404 


03020586 


0.0020586 


0.0018668 


3 


0.49926955 


0.50551097 


0.5023903 


00023903 


05024965 


0.0024965 


0.0012247 


4 


0-49980358 


0-50579084 


0.5027972 


0.0027972 


0.5029328 


0.002932S 


0.0014704 


5 


0-49995181 


0JOS18419 


0502568 


0.002568 


0.5021969 


0.0021969 


0.0015634 


6 


0.49812857 


0J0567376 


0.5019012 


0.0019012 


05021544 


0.0021544 


0.0008269 


7 


0.4999855 


0.50513723 


0.5025614 


0.0025614 


0.5030277 


0.0030277 


O.0O15744 


8 


0.49980421 


0,50565663 


0.5027304 


0.0027304 


05026188 


0.0026188 


0.0011699 


9 


0.49987774 


0.5O538874 


05026332 


0.0026332 


03022356 


0.0022356 


0.0017834 


10 


049960093 


050522843 


0.5024147 


0.0024147 


03023244 


0.0023244 


0.0014075 


11 


0.49991374 


0.50523443 


O502574J 


00025741 


03030171 


0.0030171 


0.0014693 


12 


0.49874722 


0.5050154 


0.5018813 


0.0018813 


05022298 


00022298 


0.0012337 


Average 


0.49940304 


0.505408863 


0.502406 


0002406 


03024453 


0.0024453 








b. 










Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std 


1 


0.49928084 


0.50534114 


0.502311 


0.002311 


03022704 


0.0022704 


0.0009459 


2 


0.49910! 86 


0.50500907 


0.SO20555 


0.0020555 


03019179 


0.0019179 


0.0012485 


3 


0.49954836 


0.5053165 


Q3024324 


0.0024324 


03030002 


0.0030002 


0.0014989 


4 


0.49939947 


050512739 


OJ022634 


0.0022634 


03022939 


0.0022939 


O00I5699 


5 


0.49936191 


0.SOS43349 


03023977 


0.0023977 


03028966 


0.0028966 


0.0013283 


6 


0.49995425 


050508228 


03025183 


0-0025183 


O50Z2597 


0.0022597 


0.0014881 


7 


0.4996912 


0.5O504O13 


03023657 


0.0023657 


05027383 


0.0027383 


0001482 


8 


0.49930976 


050525148 


0.5022806 


0.0022806 


03023218 


0.0023218 


0.0012336 


9 


049980147 


050599519 


0.5028983 


00028983 


05029607 


00029607 


0.0011435 


10 


04997875I 


050534796 


05025677 


0.0025677 


0302793 


0.002793 


0.0012274 


11 


0.49958979 


0.50560502 


05025974 


0.0025974 


05026188 


0.0026188 


0.0012875 


12 


0.49992005 


0.50501645 


03024683 


0.0024683 


03026511 


0.0026511 


0.0016061 


Average 


049956221 


0505297175 


OS024297 


0.0024297 


03025602 


0.0025602 





(Cont.) c. N, = 21./ = 0 



Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Eaor 


Mean 


Error 


Std 


I 


0,49940835 


0.50521986 


0.5023141 


0.0023141 


0.5023434 


0.0023434 


0.001434 


2 


0.49991704 


0.50527108 


05025941 


0X)025941 


05025089 


0.0025089 


0.0015061 


3 


0.49951541 


040501712 


05022663 


0.0022663 


05024916 


0.0024916 


OJ0013459 


4 










0.5029658 


0.0029658 


0.0015172 


5 


0.4998058! 


OJ05006I1 


0502406 


0.002406 


05021951 


0.0021951 


0.0016855 


6 


0.49947866 


0.50550706 


05024929 


0.0024929 


0502355 


0.002355 


0.001274 


7 


0.4998U17 


0.50501519 


05024132 


0.0024132 


05022662 


0.0022662 


0.0015969 


8 


0.4998184 


0.50530213 


0.5025603 


0.0025603 


0.5025368 


0.0025368 


0.0012457 


9 


0,49898706 


0J0522O8 


0502T039 


OJQ(Q1039 


0.5022881 


0.0022881 


0.0011116 


10 


0.49967244 


0.50510212 


05023873 


0.0023873 


05025692 


0.0025692 


0.0015299 


11 


0.49940494 


0.50500673 


050220158 


0.0022058 


05022872 


0.0O22872 


0.0012085 


12 


0.499$8991 


050524911 


0-^195 


0:0026195 


0.5029167 


0.0029167 


0.0015229 


Average 


0.49964909 


0505193728 


05024214 


0:0024214 


0502477 


0.002477 





Treatment 


FTSIA 


Statistic^ Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49983935 


0.50505835 


05024489 


0.0024489 


0.5023792 


0.0023792 


0.0013202 


2 


0.49962378 


050526936 


05024466 


0.0024466 


0.5024579 


0.0024579 


0.0013043 


3 


0.49933206 


0.50502746 


0.5021798 


0.0021798 


05020785 


0.0020785 


0.001371 


4 


0.49998649 


05050075 


0.502497 


0.002497 


0.5026989 


0.0026989 


0.0014821 


5 


0.49971876 


0.50512297 


0.5024209 


0.0024209 


05023309 


0.0023309 


0.0011845 


6 


0.49996842 


0.5050165 


05024925 


0.0024925 


0.5024155 


0.0024155 


0.0016019 


7 


0.49995932 


0.50506687 


05025131 


0.0025131 


05029196 


0.0029196 


0.0014238 


8 


0.49998378 


0.50516355 


05025737 


00025737 


0.5025107 


0.0025107 


0.0015578 


9 


0.49995396 


0,5050368 


05G24954 


0.0024954 


0502648 


0.002648 


0.0014875 


10 


0.49976186 


050504613 


0502404 


0.002404 


05023729 


0,0023729 


0.0015783 


11 


0.49979355 


050517038 


0502482 


0.002482 


0.5027683 


0.0027683 


0.0014278 


12 


0.49958318 


0.50504776 


0.5023155 


0.0023155 


0.502562 


0.002562 


0.0013961 


Average 


0.49979204 


0.505086136 


05024391 


0.0024391 


05025119 


0.0025119 








e. 


^, = 36./ = 0 








Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49981836 


050506372 


0502441 


0.002441 


0.5(03291 


0.0023291 


0.0013947 


2 


0.49993868 


0.50506525 


0502502 


0.002502 


0,5025447 


0.0025447 


0.0014281 


3 


0.49992337 


0.50508638 


05025049 


0.0025049 


05023862 


0.0023862 


0.0014936 


4 


0.49995917 


05050574 


05025083 


0.0025083 


05028427 


0,0028427 


0.0015979 


5 


0.49992327 


0.50500978 


0.5024665 


0.0024665 


0.5026589 


0-0026589 


0.001589 


6 


0.49982858 


050513719 


05024829 


0.0024829 


05023423 


0.0023423 


0.0014108 


7 


0.49990614 


0.50525254 


0.5025793 


0.0025793 


0.5025357 


0.0023357 


0.0013777 


8 


0-49987083 


050517544 


0.5025231 


O.0O25231 


05022054 


0.0022054 


0.0015986 


9 


0.49983627 


050563913 


05027377 


0.0027377 


0.5025759 


0.0025759 


0.0012234 


10 


0.4998221 


05052905 


0.5025563 


0^25563 


0.5028588 


0.0028588 


0.0013412 


11 


0.49996308 


0.5051605 


0.5025618 


0.0025618 


0.5025942 


0.0025942 


0.0014594 


12 


0.49987818 


050517756 


0.5025279 


0.0025279 


0.5026436 


0.0026436 


0.001426 


Average 


0.499889 


0.505176283 


0.5025326 


0.0025326 


0.5025431 


0.0025431 





Table C.2 Experiment Results of FTSIA and St^stical Calculaticnis (Case 2) 



a- JV, =8./, «l 



Treatment 




FTSIA 


Statistical Calculatioss 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0,49936145 


0.5065619 


0.5029617 


0.0029617 


03038056 


0.0038056 


0.0027807 


2 


0.49934696 


0-50627515 


03028111 


0,0028111 


0303102 


0.003102 


0..»22272 


3 


0.4989259 


0.50669498 


0.5028104 


0.0028104 


03032258 


0.0032258 


0.CO2438I 


4 


0.49747933 


0.50552288 


03015011 


Q.0015011 


0.5017775 


0.0017775 


0.COI88I9 


5 


0.497396 


0.50590588 


0.5016509 


0.0016509 


03022772 


0.0022772 


0.C022951 


6 


0.49807041 


0.50S41893 


03017447 


0.0017447 


03020745 


0.0020745 


0.CC23422 


7 


0.49977289 


0.50549349 


0.5026332 


0.0026332 


03025069 


0.0025069 


0.GQZ762 


8 


0.499332 U 


0.50597289 


0.5026525 


0.0026525 


03030647 


0.0030647 


0.CC30047 


9 


0.49952402 


0.50649&48 


0.5030113 


0.00301 13 


03033534 


0.0033534 


0.0020561 


10 


0.49864137 


0.50733809 


03029897 


0.0029897 


03033138 


0.0033138 




ii 


0.49894644 


0J0695258 


03029495 


0.0029495 


03029276 


0.0029276 


0.0020002 


12 


0.49914953 


. 0:505^568 


03024676 


0^24676 


0.5032767 


0.0032767 


0.jD02967I 


Average 


0.49982887 


I O.50i^l744 


ojmsm 


9.0025153 


0.^)28921 


0.0028921 










Treatment 




FTSIA 


Statistical CalculatioriS 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49914953 


0.50679545 


0.5029725 


0.0029725 


03039017 


0.0039017 


0.CO35406 


2 


0.49951886 


0.50719015 


03033545 


0.0033545 


0.5036211 


0.0036211 


0.>D026768 


3 


0.49991709 


0.5064845 


0.5032008 


0.0032008 


03039179 


0.0039179 


0-C032475 


4 


0.49927258 


0,50763448 


03034535 


0.(0034535 


03039449 


0.0039449 


0.C022883 


5 


0.49886323 


0.50583714 


03023502 


0.0023502 


0.5030506 


0.0030506 


0.0027044 


6 


0.49858124 


0.50610477 


0302343 


0.002343 


0.5035919 


0.0035919 


0.CQ3103T 


7 


0.49912138 


0-50636803 


03027447 


0.0027447 


0.5033371 


0.003337] 


0.0030249 


8 


0.49856593 


0-50667595 


03026209 


0.0026209 


0.5033062 


0.0033062 


0.0032311 


9 


0.49770695 


0.50635045 


0.5020287 


0.0020287 


0.5032618 


0.0032618 


O.0O36O4 


iO 


0.49798286 


0.50603726 


0.5(^0101 


0.0020101 


030337 


0.00337 


0.G033615 


H 


0.49966803 


0-50701709 


a50334i6 


0.0033426 


0.5038428 


0J0O38428 


0X028248 


12 


0.49948746 


0.50736066 


03034241 


0.0034241 


03038157 


0.0038157 


O.0O32612 


Average 


0.49898626 


0.506654661 


03028205 


0.0028205 


03035801 


0.0035801 





(^7 



Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


EnOT 


Std. 


1 


0.49891307 


0.50824385 


0.5035785 


0.0035785 


0.5O4157 


0.004157 


O.0O32255 


2 


0.499 182 1 4 


0.50884185 


0.504012 


0,004012 


0.504896 


0.004896 


0.0032092 


3 


0.49981912 


0S0966I4 




U.VU4/403 


03058456 


0.0058456 


0.0026614 


4 


0.49703591 


0.50699213 


0302014 


0.002014 


03037965 


O0037565 


0.0036467 


5 


0.497S694 


0.50710928 


0.5023393 


0.0023393 


03033403 


00033403 


0.0028803 


6 


0.49897776 


0J0863352 


03038056 


0.0038056 


03043195 


0.0043195 


0.0034355 


7 


0.4976622 


0,5075267 


03025945 


0.0025945 


03043287 


00043287 


00034448 


8 


0.49901937 


050787969 


0.5034495 


0.00344195 


03043575 


0.0043575 


0.0029468 


9 


0.49961377 


0.50791225 


0.503763 


0.003763 


03044677 


0.0044677 


0.0029U7 


10 


0.49885962 


0.50877049 


0.5038151 


0,0038151 


05045478 


0:6045478 


O0029251 


11 


0.49999813 


0,50798293 


03039905 


00039905 


03048647 


00048647 


0.0033054 


12 


0.49895114 


0.50727876 


0.503115 


0.003115 


03039462 


0.0039462 


0.0032327 


Average 


0.49880014 


0.508069404 


03034348 


0.0034348 


03044056 


0.0044056 





d. N, = 15,/, = I 



Treainient 




FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49909356 


0.50572427 


O5024089 


O0024G89 


03026127 


0.0026127 


0.0022849 


2 


049932539 


0.50542086 


0.5023731 


0,0023731 


03026977 


00026977 


0.0021969 


3 


0.4996285 


0,5055161 


0,5025723 


0,0)25723 


03028187 


0.0028187 


0.0018977 


4 


049950 1 7 


0.50541406 


05024579 


0.0024579 


0502484 


0.002484 


0.002241 


5 


0.4997 13 H 


0305589 


0.5026511 


0.0026511 


03028167 


00028167 


0.0023155 


6 


049950075 


0.50542465 


05024627 


00024627 


03025368 


0.0025368 


0.0020929 


7 


0.49865443 


0.50570372 


0.5021791 


0.0021791 


05025509 


0.0025509 


0.0017688 


8 


O49991073 


030520296 


05025569 


00025568 


03030292 


00030292 


0.0026083 


9 


0.49959115 


030573746 


03026643 


0.0026643 


03028845 


00028845 


0.0020643 


10 


0.49980359 


0.50537178 


05025877 


00025877 


0.5025324 


0.0025324 


0.0021395 


11 


049984095 


0.50575773 


05027993 


0-0027993 


03029165 


0.0029165 


0.0024803 


12 


0.49986476 


0.50569074 


03027778 


0.0027778 


03029558 


00029558 


0X)01957 


Average 


0.49953572 


0.505546111 


0.5025409 


0.0025409 


03027363 


0.0027363 








e. 










Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


I 


0.49970868 


O50682646 


03032676 


0.0032676 


05035006 


0-0035006 


O002680I 


2 


0.49942348 


030629661 


O502860I 


Q.002860I 


05034681 


0.0034681 


0.0026445 


3 


0.49976561 


0.50742117 


0.5035934 


0.0035934 


0.5041834 


0.0041834 


0.0024974 


4 


0.49888065 


030723518 


03030579 


0.0030579 


05033432 


0.0033432 


0.0018143 


5 


0.49882218 


0.50593075 


03023765 


0.0023765 


0.5032316 


0.0032316 


00024873 


6 


049937553 


0.50575139 


0.5025635 


00025635 


03033534 


0.0033534 


0.0028113 


7 


0.49982171 


0.50657344 


O5031976 


O003I976 


03033475 


00033475 


00024491 


8 


0.4995684 


0.50616925 


03028688 


0.0028688 


05035307 


0-0035307 


0.0029544 


9 


049942638 


0.50642014 


03029233 


0.0029233 


0.5032869 


O0032S69 


0.0026499 


10 


049959378 


0.50661158 


05031027 


0-0031027 


O5036319 


0.0036319 


0.0022404 


11 


049970677 


O5062J373 


05029603 


00029603 


O50340I7 


O00340I7 


000236 


12 


049945674 


030623522 


0302846 


0002846 


03032195 


0-0032 1 95 


00026946 


Average 


0.49946249 


0.506473743 


03029681 


0.0029681 


0.5034582 


0.0034582 
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(Cont) f. ^^^,5y^«5 



1 iCAtllJCIll 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49987083 


0.50746842 


0.5036696 


a0036696 


03043337 


0,0043337 


0.003777 


2 


0.4993184 


0.5084S492 


0.5039017 


a.0Q39O17 


03043518 


0.0043518 


0.0028738 


3 


0,49963695 


0.50687051 


0^032537 


0.0032537 


03043026 


0.0043026 


0.0035155 


4 


0.49941426 


0.50755314 


0.5034837 


0.0O34837 


03039954 


0.0039954 


0.0027635 


5 


0-4995663 


0.50722234 


0.5033943 


0.0033943 


03044185 


0.0044185 


0.0028373 


6 


0.49984025 


0.507 J0237 


0.5034713 


0,0034713 


0.5042335 


00042335 


0.0032507 


7 


0.49998819 


0.50636643 


0.5031773 


0.0031773 


03041108 


0.0041108 


0.0034147 


S 


0,49974628 


0.50749325 


0.5036198 


0.0036198 


0.5039659 


0.0039659 


0.0022068 


9 


0.499I9T71 


0.50736576 


0.5032817 


0.0032817 


0304167 


0-004167 


0.0029304 


10 


0.49998293 


0.5083734 


0.5<W1782 


0.0041782 


G3046221 


0.0046221 


0.0028168 


11 


0.49980767 


0.50705643 


0.5034321 


0.0034321 


0304032 


0.004032 


0.002995 


12 


0.49928838 


0.50718763 


0.503238 


0.003238 


0.5044149 


0.0044149 


0.003177 


Average . 


0.49963818 


0.507378717 


0.5035085 


0.0035085 


03042457 


0.0042457 








& 










Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Errcff 


Mean 


Error 


Std. 


1 


0.49908182 


0.50908182 


OJ040818 


0.0040818 


03046128 


0.0046128 


0.0035794 


2 


0,49880045 


0.50880045 


0.5038;00S 


0.0038005 


03042373 


0.0042373 


0^29527 


3 


0,49890159 


0.50890159 


0.5039016 


0.0O39O16 


0.5055221 


0.0055221 


0,0030501 


4 


0.49854646 


0.50854647 


0.5035465 


0.0035465 


03043975 


0.0043975 


0.0035102 


5 


0.4992234 


0.5092234 


0.5042234 


0.0042234 


0.5052727 


0.0052727 


0.0026531 


6 


0.49899324 


0.50899324 


0.5039932 


0.0039932 


03045131 


0.0045131 


0.0033109 


7 


0-49931998 


0.50931998 


0.50432 


0.00432 


03050514 


0.0050514 


0.0029778 


8 


0.49892765 


0.50892765 


a50392T7 


0.0039277 


0.5051088 


0.0051088 


0.0030651 


9 


0,49942244 


0,50942244 


0.5044224 


0.0044224 


03049713 


0.0049713 


0.0030557 


10 


0.49986701 


0.50986701 


0.504867 


0.004867 


0.5053539 


0.0053539 


0.0029991 


11 


0.49936981 


0.50936981 


O.S043698 


0.0043698 


030S1178 


0.0051178 


0.0029175 


12 


0.49950319 


0.50950319 


03045032 


0.0045032 


0.5059093 


0-0059093 


0.0030103 


Average 


0.49916309 


0.509163088 


0.5041631 


0.0041631 


03050057 


0.0050057 








h. 


/V,=21,/, =1 








Treatment 


FrSL\ 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49920187 


0.5050874 


0.5021446 


0.0021446 


03021843 


0.0021843 


0.0015967 


2 


0.49907606 


0.50541744 


0.5022468 


0.0022468 


03025643 


0.0025643 


0.0018831 


'X 
D 


0.49967122 


0.50554719 


0.5026092 


0.0026092 


0.5027374 


0.0027374 


0.0020205 


4 


0.49967618 


0.50569237 


0.5026843 


0.0O26843 


03027957 


0.0027957 


0.0019476 


5 


0.49978336 


0.50503005 


0.5024067 


0.0024067 


03024649 


0.0024649 


0.0020494 


6 


0.49944217 


0.505446 


0,5024441 


0.0024441 


0502667 


0.002667 


0:0016343 


7 


0,49978497 


0.5058003 


0.5027926 


0.Q027926 


0.5028382 


0.0028382 


0.002O723 


8 


0.49990121 


0.5052 1825 


0.5025597 


0.0025597 


03027)88 


0.0027188 


0.0017539 


9 


0.49995088 


0.50534043 


0.5026457 


0.0026457 


0.5028989 


0-0028989 


0.0019951 


10 


0.49989302 


0.505159 


0.502526 


0.002526 


0.5024061 


00024061 


0.0018107 


11 


0.49987116 


0.50503286 


0.502452 


0.002452 


0.5029285 


00029285 


00019622 


12 


0.49969013 


0.50592841 


0.5028093 


0.0028093 


0.5028085 


0.0028085 


0,0019928 


Average 


0.49966185 


0.505391642 


03025268 


0.0025268 


03026677 


0.0026677 





Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


EiTor 


Mean 


Error 


Std. 


i 


0.4996278 


0.50539909 


0.5025135 


0.0025135 


0.5026392 


0.0026392 


0.0025234 


2 


0.49946925 


0.50674555 


0.5031074 


0.0031074 


0.5033072 


0.0033072 


0.0022555 


3 


0,4999579 


0 50581314 








0.0030096 , 


0.0023899 


4 


0.49912807 


0.50530749 


0.5022178 


0.0022178 


0.5026453 


0.0026453 


0.0027413 


5 


0.49938447 


0.50596042 


0J026725 


0.0026724 


OJ029878 


0.0029878 


0.0021412 


6 


0.499651 


0,50544031 


0.5025457 


0.0025457 


0 J029696 


0.0029696 


0.0024379 


7 


0.49978823 


0.50629452 


0^030414 


0.0030414 


0.5034586 


0.0034586 


0.0024426 


8 


0.499!431S 


0.50565307 


a50239Sl 


0.0023981 


0.51^6251 


0J(mS'2Sl 


0.0022132 


9 


0,49925369 


0,50548495 


OJ023693 


0.0023693 


0JO27O15 


0.0027015 


0.0023041 


10 


0.49986755 


0.5060F061 


03029391 


a0029391 


0.503€^ 


0.0630544 


0.0022048 


11 


0.49947339 


0.50626276 


05028681 


0.0028^1 


05(^809 


0.0029809 




12 


0.49918695 


0.50567518 


0.5O243U 


0,0024311 


0.5024963 


0.0024963 


I 0.0021883 


Average 


0.49949429 


0.505837258 


0.5026658 


0.0026658 


0.5029063 


0.0029063 





j. AT, = 21./, =5 



Treatment 




FTSIA 


Statistical Ca:lculation$ 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std 


1 


0.49989728 


0.50573381 


0.5028156 


0.0028156 


0.5037183 


0.0037183 


0.0027985 


2 


0.49923914 


0.50620662 


0.5027229 


0.0027229 


0.503336 


0.003336 


0.0029422 


3 


0-49909696 


0.50669485 


0.5028959 


0.0028959 


0.5038178 


a0038178 


0.0030921 


4 


0.49931139 


0-50651872 


OJ502915I 


0.0029151 


0.5036555 


0.0036555 


0.0023377 


5 


0.49904724 


0 50688365 


0.5029655 


0.0029655 


0.5035252 


0.0035252 


0.0023126 


6 


0.49866755 


0.50642247 


0.502545 


0.002545 


0.5035636 


0.0035636 


0.0028635 


7 


0.49937126 


0,50651502 


OJ029431 


0.0029431 


0.5035302 


0.0035302 


0.0027819 


8 


0.49891775 


0.50678096 


0.5028494 


0.0O28494 


0.5033285 


0.0033285 


0.0Qn699 


9 


0.4991 188 


0.50583317 


0.502476 


0.002476 


0.5030478 


0:0030478 


0.002568 


10 


0.49939535 


0.50577655 


0.502586 


0.002586 


0.50^75 


0.0034675 


0.0031506 


11 


0.49930806 


0.50609451 


O.5027O13 


0.0027013 


03032721 


0.0032721 


0.0025768 


12 


0.49987071 


0.50631505 


0.5030929 


0.0030929 


0.5035321 


0.0035321 


0.0027859 


Average 


0.49927012 


0.506314615 


0.5027924 


0.0027924 


0.5034829 


0.0034829 








k. 


AT, = 21./, =7 








Treatment 




FTSIA 


Sutistical Calculations 


Lx>w bound 


High bound 


Mean 


EiTor 


Mean 


Error 


Std. 


1 


0.49981846 


0.50787374 


0.5038461 


0.0038461 


0.5041158 


0.0041158 


0.0028689 


2 


0.49981248 


0.50843261 


0.5041226 


0.0041225 


0.5046395 


0.0046395 


0.0027877 


3 


0.49991319 


0.50794918 


0.5039312 


0.0039312 


0304658 


0.004658 


0.0030199 


4 


0.49980047 


0.50826195 


04040312 


0.0040312 


0.5043475 


0,0043475 


0.0024763 


5 


0.49942243 


0.50777226 


0.5035974 


0.0035973 


0.5041457 


0.0041457 


0.0029134 


6 


0.49939799 


0 50814742 


0JQ37727 


0.0037727 


05045666 


0.0045666 


0.0027338 


7 


0.4988308 


0.50750597 


05031684 


0.0031684 


0.5043382 


0.0043382 


0.0033891 


8 


0.49930107 


0.50811252 


0.5037068 


0.0037068 


0.50425 


0.00425 


0.0027992 


9 


0.49905259 


0.50728429 


0.5031684 


0.0031684 


0.503979 


0.003979 


0.0028866 


10 


0.49998888 


0.50734449 


0.5036667 


0.0036667 


0.5043319 


0.0043319 


0.0030015 


11 


0.49990333 


0.50706028 


OJ0348I8 


0.0034818 


0.5043224 


0.0043224 


0.0033459 


12 


0.4991441 


0-50753506 


0.5033396 


0.0033396 


0.5040821 


0.0040821 


0.0029881 


Average 


0.499532.15 


0.507773314 


0J036527 


0.0036527 


0.5043147 


0.0043147 





n 



Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 
i 


0.49968267 


0.50968267 


0J046827 


0.0046827 


03Q51SQ5 


0.0051505 


0.002984 




0.49937428 


0.50937428 


03043743 


0.0043743 


03056529 


0.0056529 


0.0030123 


'1 

J 


0.49937428 


0.50937428 


0J043743 


0.0043743 


03052322 


0.0052322 


0.0027451 


4 


0.49988151 


0.50988151 


0.5048815 


0.0048815 


0.5055707 


0.0055707 


0.0028647 


5 


0.49876603 


0.50876603 


0J03766 


0.003766 


0.5047361 


0.0047361 


0.0025604 


6 


0.49889554 


0.50889554 


0.5038955 


0.0038955 


03049002 


0.0049002 


a003453 


7 


0.49950986 


OJ0950986 


0.5045099 


0.0045099 


03048014 


0.0048014 


0.003317 


8 


0.49948509 


0J0948509 


0.5044851 


0.0044851 


0.5049487 


0.0049487 


OiX>32447 


9 


0.49922582 


0-50922582 


03042258 


0.0042258 


03049345 


0.0049345 


0.0031369 


10 


0.49947976 


0.50947976 


0.5044798 


0.0044798 


0.5048989 


0.0048989 


0.0028323 


11 


0.49933787 


0.50933787 


03043379 


0.0043379 


03051951 


0.0051951 


0.0032867 


12 


0-49953019 


0.50953019 


03045302 


0.0045302 


0305419 


0.005419 


0.0029089 


Average 


0.49937858 


0.509378575 


03043786 


0.0043786 


030512 


0.00512 
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Table C.3 Experiment Results of FTSIA and Statistical Calculations (Case 3) 



Treatment 


FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 




Mean 


Error 


Std. 


1 
1 


0.49830274 


0.50525146 


0.5017771 


0.0017771 


0.5075333 


0.007S333 


0.016518 


2 


0.49987353 


0.50501107 


0.5024423 


0.0024423 


0.5080563 


Q.0080563 


0.0163823 


3 


0.49945923 


0^513215 


0J022957 


0.0022957 


0.507739 


0.007739 


00134797 


4 


0.49842174 


0.50540274 


0.5019122 


0.0019122 


0.5058399 


0.0058399 


0.0123811 


5 


0.49915461 


0.50562071 


0.5023877 


0.0023877 


0.5069586 


0.0069586 


0.0121969 


6 


0.49909057 


0.50560486 


0.5023479 


0.P023479 


0.5055159 


0.0055159 


Q.0099107 


1 


0.49924392 


0.50507159 


0:5021578 


0.0021578 


0.5060495 


0.0060495 


0.0095769 


8 


0.49881164 


0.5052022 


OJ020069 


0.0020069 


03049012 


0.0049012 


0.0074582 


9 


0.49980127 


0.50507573 


0J024385 


0.0024385 


0.5050885 


0.0(^0885 


0.006S384 


10 


0.49921703 


0.50522586 


05022215 


0.0022214 


0.5070597 


0.0070597 


0.0137379 


11 


0.49977707 


. 04P503Q31 


.0.5024037 


0.0024037 


0JO7041 


a007041 


0.013193 


12 


0.49901493 


0.5067047 


0J028598 


0.0028598 


0.5061998 


0.0061998 


0.0108978 


Average 


0.49918072 


0.505361115 


0.5022709 


0.0022709 


0.5064986 


0.0064986 





b. AT, =8./. =2 



Treatment 




FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49870636 


05055557 1 


0.502131 


0.002131 


0.5110072 


0.0110072 


0.0176853 


2 


0.49819933 


0.50599439 


0.5020969 


0.0020969 


0.5084628 


0.0084628 


0.0125514 


3 


0.49888227 


0.5062481 1 


0.5025652 


00025652 


0.5109883 


0.0109883 


0.0157251 


4 


0-49912602 


0.50571385 


0.5024199 


0:0024199 


0.5111247 


0.0111247 


0.0163524 


5 


0.4996572 


0.5051419 


0:5023996 


00023996 


OJU0628 


0.0110628 


00163108 


6 


0.49878538 


0.50665322 


0.5027193 


0.0027193 


0.5119823 


0.0119823 


0.0172604 


7 


0.49903478 


0.50550432 


0.5022696 


00022696 


0.5102407 


0.0102407 


0.0153966 


8 


0.4986346 


0.50515685 


0.5QI8957 


0.0018957 


0.5119636 


0.0119636 


00186969 


9 


0.49884751 


0.50540094 


0.5021242 


O0021242 


0.5081833 


0.0081833 


0.0112803 


10 


0.4994794 


0.50549011 


0.5024848 


00024848 


0.5103879 


0.0103879 


0.0162172 


11 


0.49802428 


0.50648732 


0-5022558 


0.0022558 


0.508039 


0.008039 


00105507 


12 


0.49984907 


0.50514029 


0.5024947 


0.0024947 


0.5134565 


00134565 


0.0192871 


Average 


049893552 


0.505707251 


05023214 


0.0023214 


0.5105749 


0.0105749 





13^ 



(Com.) 



Treatment 


FTSIA i Statistical Calculations 


Low bound 


High bound 


Mean 


Error Mean 


Error 


Std. 


I 


0.49944543 


0.50501754 


05022315 


0.0022315 05131801 


0.0131801 


0,0165447 




0.49923084 


0.50618961 


0.5027102 


0.002710:- 05163098 


0.0163098 


0,0195086 


J 


0.49948915 


0.50595748 


0.5027:53 


0.0027233 : 05164627 


0.0164627 


0.0188129 


4 


0.49791586 


0.50627818 


0.502097 


0.002097 05146064 


0.0146064 


0^)174466 


5 


0.49954359 


0.50502469 


a5022S41 


0.0022841 ' 05135391 


0.0135391 


0.0154008 


6 


0.49913236 


0.50500461 


0.5020685 


0.0020685 05119023 


0.0119023 


0.0156807 


7 


0.49901769 


0^0609031 


0.502554 


0.002554 0512231 


0.012231 


0.0138531 


8 


Q.4997669 


0.50593613 


0.5028515 


0.0028515: 05137721 


0.0137721 


0.0169249 


9 


0,49918897 


0.50589985 


0.5025^44 


0.0025444 05140257 


0.0140257 


O.0168013 


10 


0-49748344 


0.50518901 


0.5013362 


0.0013362 05116034 


0.0116034 


0.013842 


11 


0.49916516 


0.50598235 


0.5025-38 


0.Q02S738; 0513901 


0.013901 


0.0176S 


12 


0.49942671 


0.50573576 


0.5025S12 


0:0025812? 05170181 


0.0170181 


0.0203459 


Average 


0.49906718 


0.505692127 


0.5023'96 


0.0023796 - 0514046 


0.014046 








Treatment 


FTSIA Statistical Calculations 


Low bound 


High bound 


Mean 


Error Mean 


Error 


Std. 


1 


0.49918766 


0,50630671 


0.5027472 


0.0027472 ; 05043665 


0.0043665 


0.0093878 


2 


0.49948799 


0.50517088 


05023294 


0.0023294 : 05043665 


0.0043665 


0.0075608 


3 


0.499933 


0.50504101 


0.502487 


0.002487 ; 0.504:1741 


0.0041741 


0.0066408 


4 


0.49973235 


0.50567204 


0.5027022 


0.0027022. 05048323 


0.0048323 


0.0083196 


5 


0.49961401 


OJ052569S 


0.5024355 


0.0024355 ! 05045129 


0.0045129 


0.0085679 


6 


0.49974381 


0.50503725 


0.5023905 


0.0023905 : 0.5034205 


0.0034205 


0.0053537 


7 


0.49937253 


0.50517716 


05022T49 


0.0022748 0.5050986 


0.0050986 


0.0109576 


8 


0.4988081 


0.50502817 


0.5019131 


0.0019181 05041521 


0.0041521 


0.0068359 


9 


0.49881568 


0.50502514 


05019:(M 


0.0019204 ^ 05046925 


0.0046925 


0.0119113 


10 


0.49947265 


0.50502651 


05022496 


0.0022496 0.5040509 


0.0040509 


0.0074641 


11 


0.49930384 


0.505237 


05022-04 


0.0022704 0 5057298 


0.0057^8 


0.011275 


12 


0.49934828 


0.50505565 


0502202 


0.002202 05041576 


0.0041576 


0.008513 


Average 


0.49940166 


0.505252873 


0.50232T3 


0:0023275 0.5044629 


0.0044629 








Treatment 


FTSIA Statistical Calculations 


Low bound 


High bound 


Mean 


Error : Mean 


Error 


Std. 


1 


0.49900849 


0.50560675 


0.50230^6 


0.0023076 05091692 


0.0091692 


0.0237996 


2 


0.49949217 


0.50514396 


0.5023181 


0.0023181 0.5098424 


0.0098424 


0.0159538 


3 


0.49939763 


0.5050199 


05022088 


0.0022088 0.5077875 


0.0077875 


0.0129313 


4 


0.49886258 


0.50511814 


05019904 


0-0019904 0.5080301 


0.0080301 


0.0137491 


5 


0.49901183 


0.50516398 


050208T9 


0.0020879 05087581 


0.0087581 


0.0135335 


6 


0.4992263 


0^50507799 


0.5021522 


0.0021521 0.508966 


0.008966 


0.0145932 


7 


0.499079 


0.5050282 


05020536 


0.0020536 0.5097317 


0.0097317 


0.0170256 


8 


0.49912469 


0.50500875 


0.5020667 


0.0020667 0.5080259 


0.0080259 


0.0128308 


9 


0.49921088 


0.50532749 


0.5022692 


0.0022692 0.5092227 


0.0092227 


0.0152136 


10 


0.49896978 


0.50518432 


0.5020^1 


0.0020771 05079268 


0.0079268 


0,013086 


11 


0.49934466 


0.50507624 


0.5022105 


0.0022105 0.5090934 


0.0090934 


0.0155553 


12 


0.49958023 


0.50506043 


05023203 


0.0023203 0.5098407 


0.0098407 


0.015221 


Average 


0.49919235 


0.505151346 


05021 "19 


0.0021719 0.5088662 


0.0088662 
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Treatment 




FTSLA 


Statistical Calculations 


Low bound 


rljgn bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49917667 


0.505 1 1 148 


0.5021441 


0.0021441 


0.5158601 


0.0158601 


0.019877 


2 


0.49928281 


. 0^0566513 


0.502474 


0.002474 


03145221 


0.0145221 


0.0177471 


3 


0.49959884 


0.50543568 


05025173 


0.0025173 


03122182 


0.0122182 


0.015674 


4 


0,49934423 


0.50626777 


0.502806 


0.002806 


03134236 


0.0134236 


O.O170806 


5 


0.49918115 


0.5055491 


0J023651 


0.0023651 


03108567 


0.0108567 


0.01309 


6 


0.4994214 


0.505051 


0.5022362 


0.0022362 


0311652 


0.011652 


0.0143348 


7 


0.4995079 


0.50506496 


0.5022864 


0.0022864 


03105026 


0.0105026 


0.013805 


8 


0.49886887 


0.505S4219 


0.5023555 


0.0023555 


03121111 


0.0121111 


0.0151241 


9 


0.49961189 


0.50533472 


0.5024733 


0.0024733 


03129126 


0.0129126 


0.0159438 


10 


0.4998642 


0.50500773 


O.502436 


0.002436 


03120777 


omiom 


0.0161682 


11 


0.49814443 


0.505244 


0.5016942 


0.0016942 


03101859 


0.0101859 


0.0134742 


12 


0.4994554 


0150613233 


03027939 


0.0027939 


03133779 


0.0133779 


0.0175827 


Average 


0.49928815 


0^05475508 


03023818 


0.0023818 


0312475 


aOl2475 














Treatment 




FTSU 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49921 147 


0.50538233 


0.5022969 


0.0022969 


03134034 


0.0134034 


0.0132934 


2 


0.49953193 


0J0501715 


0.5022745 


0.0022745 


03175264 


0.0175264 


0.0181207 


3 


0.49906269 


0.50608634 


0.5025745 


0.0025745 


03140893 


0.0140893 


0.015053 


4 


0.49908672 


0.50518782 


0.5021373 


0.0021373 


03142892 


0.0142892 


0.0159702 


5 


0.49928674 


0.50512343 


0-5022051 


0.0022051 


0.5149995 


0.0149995 


0.0159002 


6 


0.49862235 


0.50561659 


0.5021195 


0.0021 195 


0315003 


0.015003 


0.0146906 


7 


0.49888998 


0.50509682 


0.5019934 


0.0019934 


03168192 


0.0168192 


0.0173196 


8 


0.49914178 


0.50505716 


0.5020995 


0.0020995 


03161403 


0.0161403 


0.0165769 


9 


0.49955719 


0i;0555008 


0.5023536 


0-0025536 


03170762 


aOl70762 


0.0183007 


10 


0.49903783 


0.50568651 


0.5023622 


0.0023622 


03153519 


0.0153519 


0.0159533 


11 


0.49961057 


0.50528602 


0.5024483 


0.0024483 


03149022 


0.0149022 


0.0171499 


12 


0.49966269 


0.50500601 


0-5023344 


0.0023344 


03173278 


0.0173278 


0.0181771 


Average 


0.49922516 


0.505341355 


0.5022833 


0.0022833 


03155774 


0.0155774 








h. 










Treatment 




FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49931695 


0-50538657 


0-5023518 


0.0O23518 


03041177 


0.0041177 


00065764 


2 


0.49944147 


0.50334463 


0.5023931 


0.0023931 


03045873 


0.0045873 


0.0087526 


3 


0.49977857 


0.5051036 


0.502441 1 


0.0024411 


0.5043139 


0.0043139 


0.0092816 


4 


0.49929358 


0-5052798 


0.5022867 


0.0O22867 


0,5043502 


0.0043502 


0.0077W4 


5 


0.49942866 


0.50511409 


0.5022714 


0.0022714 


03037087 


0,0037087 


00060528 


6 


0.49906262 


0.50505645 


0.5020595 


0.0020595 


03042384 


0.0042384 


O.0094674 


7 


0.499474 


0.50517891 1 


0.5023265 


0.0023265 


03046626 


0.0046626 


00087634 


8 


0.49933806 


0.50515041 1 


0.5022442 


0.0022442 


0.503702 


0.003702 


0.006492 


9 


0.49967052 


0.50521446 ! 0.5024425 


0.0024425 


03044678 


0.0044678 


0.0092793 


10 


0.4990245 


0.50535402 \ 


0.5021893 


0.0O21893 


0.5039407 


0.0039407 


00093746 


11 


0.49925803 


0.50530834 j 0.5022832 


0.0022832 


0.5041243 


0.0041243 


0.005212 


12 


0.4994118 


0.5050493 ! 


0.5022306 


0.0O22306 


03042305 


0.0042305 


0.0074928 


Average 


0.4993749 


0.505211715 ! 


0.5022933 


0.0022933 


0.5042037 


0.0042037 
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i. A^, =21,/, =3 







FTSIA 


Statistical Calculations 


Low bound 


xxign oouno 


Mea.n 


Error 


Mean 


Error 


Sid. 


1 


0.499268S9 


0.505432 


0.5023503 


0.0023503 


0.5077672 


0.0077672 


0.0143445 


2 


0.49946094 


0.dU9U7853 


, 0 J022697 


0.0022697 


0.5061471 


0.0061471 


O.0095461 


3 


0.49945736 


0.50508226 


0.5022698 


00022698 


0.5077436 


0 0077436 




4 


0.49877968 


050508729 


05019335 


0.0019335 


0.505906 


0.005906 


0.0108534 


5 


0.49949521 


0.50516978 


0J023325 


0.0023325 


0.5074487 


0.0074487 


0.0139683 


6 


0.49875457 


0.50524475 


05019996 


0.0019996 


0.5057807 


0.0057807 


0.0094563 


7 


0.49942347 


0.50507701 


0.5022502 


0.0022502 


0.5063741 


0.0063741 


0.0105449 


8 


0.49917074 


0.50544583 


0..5023083 


0.0023083 


0.5064334 


0.0064334 


0.0113358 


9 


0.49930134 


0.50543647 


05023689 


0.0023689 


0.5078356 


0.0078356 


0.0136095 


10 


0.49986855 


050509117 


0.5024799 


0.GO24799 


0.5070241 


0.0070241 


0.0122775 


11 


0.49910589 


0.50510477 


0.5021053 


0.0021053 


05067803 


0.0067803 


0.0123945 


12 


0.49924919 


O.505O7O19 


05021597 


0.0021597 


0.508541 1 


0.008541 1 


0.0155197 


Average 


0.49927794 


0.505193338 


0.5022356 


0.0022356 


0.5069818 


0.0069818 





J- ^, =21,/„=5 



Treatment 




FTSIA 


Statistical Calculations 


Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.4989793 


0.50524254 


0.502U09 


0.0021109 


0.5097801 


0.0097801 


0.013795 


2 


0.49941933 


OJ0533127 


0.5023753 


00023753 


0.5106901 


0.0106901 


0.0149554 


3 


0.4990898 


0.50609337 


05025916 


00025916 


0.5094831 


0.0094831 


a01343a7 


4 


0.49944274 


0.50516118 


0.502302 


0.002302 


0J09493 


0.009493 


0.0132579 


5 


0.49923013 


0.5051598 


0.502195 


0002195 


0.5083682 


0.0083682 


0.013U9 


6 


0.49915662 


0.50502213 


0.5020894 


0.0020894 


0.5102857 


0.0102857 


0.0149801 


7 


0.49910243 


0.50595729 


0.5025299 


0.0025299 


0.510803 


0.010803 


0.0157481 


8 


0.49951605 


0,50510798 


0.502312 


0.002312 


OJ10I004 


0.0101004 


0.0142097 


9 


0.49949946 


05052413 


05023704 


00023704 


0.5097258 


0.0097258 


0.0141988 


10 


0.49948672 


0.50514715 


0.5023169 


0.0023169 


0.509692 


0.009692 


0.01531 


11 


0.49923988 


050548895 


05023644 


00023644 


0.5108094 


0.0108094 


0.0162715 


12 


0.49943569 


0.50532238 


0502379 


0002379 


0.5104157 


0.0104157 


0.0151396 


Average 


0.49929985 


0.505356278 


0.5023281 


0.0023281 


0.5099705 


0.0099705 







-J — ' 


Treatment 




FTSL\ 


Statistical Calculations 




Low bound 


High bound 


Mean 


Error 


Mean 


Error 


Std. 


1 


0.49893459 


0.50549386 


05022142 


0.0022142 


0.5127394 


0.0127394 


0.0152902 


2 


0.49926369 


0.50504884 


0.5021563 


0.0021563 


0,512814 


0.012814 


0.0172981 


3 


0.49937003 


0.50512409 


05022471 


0002247 1 


0.5150132 


0.0150132 


0.0186415 


4 


0.49918837 


0.50503421 


0.5021113 


0.0021113 


0.513678 


0.013678 


0.0174939 


5 


0.49932718 


0.50504101 


0.5021841 


O.0O2184) 


0J;i41l52 


0.0141152 


0.0175016 


6 


0.49931434 


0.50532103 


0.5023177 


0.0023177 


OJ 132846 


0.0132846 


0.0164223 


7 


04992437 


0.50500753 


0.5021256 


0.0021256 


0.512574 


0.012574 


0.0169267 


8 


0.49934622 


0.50508527 


0.5022158 


0.0022157 


0.5134157 


0.0134157 


0.0177348 


9 


0.49897656 


0.50542936 


0.502203 


0.002203 


0.5128614 


0.0128614 


0.0163791 


10 


0.49939036 


0.50513913 


0.5022648 


0.0022647 


0.5119821 


0.0119821 


0.015084] 


11 


0.49897199/ 


O505 13231 


05020522 


00020522 


0.5144925 


0.0144925 


0.0186373 


12 


0.49954425' 


050502574 


0.502285 


0.002285 


0.5099386 


0.0099386 


0.0121739 1 


Average 


0.49923927 


0.505156865 


05021981 


00021981 


0.5130757 


0.0130757 1 





i 



86 



Treatment 




FTSIA 


Statistical Calcul^itioiis 


Low bound 


ffigh bound 


Mean 


Error 


Mean 


Error 


Std. 


t 
1 


0.4997133 


0.50534933 


0.5025313 


0.0025313 


0.5136752 


0:0136752 




4 


0.4997523 


0.5O542793 


0.5025901 


0.00259Q] 


0-5149228 


0.0149228 


0.0150769 


1 

o 


0.49949133 


0.50504049 


0,5022659 


0.0022659 


OJ 146835 


0.0146835 


0.0155 115 


4 


0.49932866 


0.50505773 


0J021932 


0.0021932 


05203112 


0.0203112 


0.0198055 


5 


0.49885099 


0.50566457 


0.5022578 


0,0022578 


0.5153462 


0.0153462 


0.0157767 


6 


0.49940298 


0.50595539 


0.5026792 


0.0026792 


0.5147693 


0.0147693 


0.0147205 


7 


0.49924746 


0.50519481 


05022211 


.0.0022211 


0.5156894 


0.0156894 


Q.Q151724 


8 


0.49966613 


0.50504135 


0.5023537 


0.0023537 


0,5150388 


0.0150388 


0.0)52035 


9 


0,49941179 


. Q.505.7349 


03025734 


0.0025734 


0315244 


0015244 


0.0143989 


10 


0.49897576 


0.50547452 


0.5022251 


0.0022251 


03144501 


0.0144501 


0.015345 


11 


0.49882502 


0.50539256 


0.5021088 


O.002t088 


03162709 


0.0162709 


0,0177374 


12 


0.49930065 


0.50501257 


0.5021566 


0.0021566 


0313945 


0.013945 


0.0146133 


Average - 


0.49933053 


0.505362179 


0,5023464 


0.0(^3464 


03153622 


0.0153622 





^ <^Hrror of statistical mean- Error of mean from FTSIA) 



0.012 



0.01 



0,008 



0.006 



0.004 



0.002 




Experiment 



Figure C. 1 Differences of Error of Statistical Mean and Error of Mean from FTSIA in 29 
Experiments with 12 Treatments each (Experiments 1 to 29 correspond to 
Tables C. I .a to C.3.1 respectively) 

* the point surrounded by a rectangular means the difference is significant (a=95%)in the corresoondine 
experiment ^ * 
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1. Introduction 

The measurement of fluid flow is an important technology used in nearly every industry. 
Flow-sensors are used to measure the movement of fluid flow (liquid or gas flow). They 
are widely used in a variety of applications, such as automotive applications, 
environmental monitoring, biomedical diagnostics, and industrial testing. 
Before the advent of micro minimization, sensor systems tended to be bulky and 
expensive, so single sensor systems have typically been designed to save space and 
manufacturing cost. Although single sensor systems are relatively easy to construct and 
analyze, they have limitations. For instance, if we need to measure several variables, say, 
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temperature, pressure and flow density, all at the same time, single sensor systems cannot 
be used. In addition, a single sensor cannot guarantee to deliver accurate information all 
the time, because it is inevitably affected by noise, refresh delay, and other uncertain 
disruptions. Therefore, single sensor systems are not suitable for critical applications. 
Development of MEMS has enabled production of large amounts of micro sensors at low 
cost. Microelectromechanical system (MEMS) technology deals with design and 
fabrication of entire electrical and mechanical systems, usually on a single silicon chip. 
Micro flow-sensors are one of the most common MEMS devices. Due to their small 
sizes, micro flow-sensors have been developed in recent years (Berberig et al., 1998; 
Okulan et al., 2000; Wu and Sansen, 2002) with the following advantages: 

1 . Interfere less with the environment they are measuring 

2. Lower manufacturing cost 

3. Can be applied in narrow spaces, such as inside living organisms, small pipes, and 
automobile engines. 

4. Multiple sensors can be designed with redundancy to improve information 
accuracy, responsiveness, and fault tolerance of the system. 

Distributed micro flow-sensor arrays and networks (DMFSA/N) are built from 
collections of spatially scattered, intelligent micro flow-sensor nodes. Each node has the 
ability to measure the local flow within its accuracy limits, process the raw sensor data, 
and cooperate with its neighboring nodes. It enables more reliable and timely monitoring 
and control. Sensors incorporated with dedicated signal processing functions are called 
intelligent sensors, or smart sensors. The main roles of dedicated signal processing are to 
enhance design flexibility and realize new sensing functions. Additional roles are to 
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reduce loads on central processing units and signal transmission lines by distributing 
information processing to the lower layers of the system (Yamasaki, 1995). Here, the 
micro flow-sensor array refers to a group of the same type of micro flow-sensors which 
are placed close together to measure the same variable of interest. Alternatively, sensors 
with different measurement functions can be combined in an array. Different groups of 
micro flow-sensor arrays can be distributed widely in the flow environment to form a 
micro flow-sensor network. 

Developing an effective DMFSA/N, however, also faces great challenges: 

1. First, since sensor nodes are distributed spatially, appropriate sensor network 
architecture has to be designed to support the efficient communication among the 
sensor nodes 

2. Second, minimizing communication cost, C, which is often related to 
communication time, C/, and energy consumption, Cg, is important in online 
control and energy constrained applications. Therefore, cost-efficient 
communication protocols need to be developed 

3. In addition, it is unrealistic to expect all the micro sensor nodes and 
communication links to be fault-free and available all the time, therefore, 
developing fault-tolerant communication protocols is another important issue in 
the design of DMFSA/N. 

The objective of this paper is to address the design of the network architectures, NA, and 
communication protocols, CP, for DMFSA/N. A new cluster network architecture, two 
alternatives of fault-tolerant time-out communication protocol (FTTP) - one with blind 
stations, the other with smart base stations, are proposed in this paper (Nof and Liu, 



3 



2003). TIE (Teamwork Integration Evaluator)/MEMS, a N4EMS sensor network 
simulation tool, has been developed to evaluate the effectiveness of the network 
architectures and communication protocols. Experimental results show that the proposed 
cluster architecture can support efficient communication in the sensor network, and the 
FTTP with smart base stations requires less communication time than the FTTP with 
blind stations, although smart base stations have more complex structures than blind base 
stations. In this paper, some current research on the design of sensor network 
architectures and communication protocols is reviewed in section II; the new cluster 
architecture and FTTP for DMFSA/N are presented in section III; in section IV, 
TIE/MEMS is described. Two sets of simulation experiments were conducted using 
TIE/MEMS to evaluate the five combinations of sensor network architectures and 
communication protocols, and the two alternatives of FTTP. The experimental results and 
discoveries are presented and discussed. Finally, conclusions, significance of the 
discoveries, and future challenges are discussed in section V. 

2. Sensor Network Architecture and Communication Protocol Review 

The abstraction of a sensor network can be represented as NA = {V^E), where V\s the set 
of sensor node and E is the set of communication links connecting the sensor nodes. 
V = = 1,2,..., ATJ, where Ns is the number of sensor nodes in the network and is 
the sensor node /. £ = {e^}(/ = l,2,...,iV,;y = l,2,...,iVJ where e,-,.is the link between 

V . and v^. . 
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The two most conventional network architectures that can be used for sensor networks 
are the committee architecture (COA), and the hierarchical architecture (HIA), first 
discussed by Wesson et al (1981). In COA, each node in the network is autonomous and 
can send information to all other nodes in the network. Thus \fNs is equal to «, it needs 
interconnections that are of the order 0(«^). With such a high number of 
interconnections, this architecture is very robust to link failures because the messages 
from a source node to a destination node can be easily rerouted by bypassing the failed 
links. However, COA is not suitable for the sensor network with a large value of Ns, 
because it is too expensive with respect to the communication cost, C, and not easily 
extensible (figure 1(a)). In HIA, nodes are connected to form strict hierarchies. Each node 
(except the lowest-level leaf node) receives information fi-om the lower level node; 
integrates the information received according to its position in the hierarchy, and then 
sends information to the node at the upper level. The root node is the node at the highest 
level and responsible for interpretation of data at the global level. A special case of HIA 
is a complete binary tree (figure 1(b)), Compared with COA, HIA is extensible more 
easily, but it is not tolerant to link failures. Two other sensor network architectures, 
proposed by Iyengar et al (1990) and Iyengar et al (1994) are flat Tree architecture 
(FTA) and Multilevel binary de Bruijn architecture (MBDA). In FTA, the network is 
partitioned into clusters and each cluster is internally organized as a binary tree. All the 
root nodes in the network are fully interconntected (figure 1(c)). FTA is not sufficiently 
robust to link failures. MBDA is a multilevel network (MLN) with the top level fully 
connected and with each of the other levels interconnected as a de Bruijn network (figure 
1(d)). MBDA is not suitable for the DMFSA/N because of its complex structure. 
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In recent years, there has been increasing interest in the development of communication 
protocols for the distributed sensor network, especially the wireless sensor network 
(Gossink et al., 1998; Agre et al, 1999; Chandrakasan et ah, 1999; Heinzelman et al, 
1999; Bhushan and Rengaswamy, 2000; Min et al., 2002). The sensor network can 
improve the accuracy of information obtained via collaboration among sensor nodes. 
The simplest communication protocol is the point-to-point protocol (PTP). In this 
protocol, sensor nodes A and B communicate exclusively with each other without 
interfering with other nodes (figure 2(a)). Another commonly used communication 
protocol is broadcasting protocol (BP). In this protocol, the information sent out by the 
sender node can be received by all its neighbors using only one transmission (figure 
2(b)). Since a node always sends data to its neighbors, regardless of whether or not they 
have already received the data from another source, BP can lead to the implosion 
problem and waste resources by sending requested copies of data to the same node. In the 
Gossiping protocol (GP), instead of indiscriminately sending information to all its 
neighboring nodes, each sensor node only forwards the data to one randomly selected 
neighbor (Heinzelman, 1999) (see figure 2(c)). While GP requires less communication 
than BP, it is not as robust to link failures as BP because a node can only rely on one 
other node to reroute information for it if its link fails. 



[Insert figure 1 here] 
[Insert figure 2 here] 



3. DMFSA/N Architecture and Communication Protocol 

In this section, the proposed cluster network architecture and fault-tolerant time-out 
communication protocol for DMFSA/N are presented. 

3.1 Cluster Network Architecture for DMFSA/N 

The new cluster network architecture for DMFSA/N, Cluster architecture (CLA) is 
shown in figure 3. It is preferred because it overcomes the major limitations of both COA 
and HIA by limiting the amount of communication, while at the same time tolerating 
failures of the sensor nodes and links. In this structure, DMFSA/N are divided into sensor 
cluster units (SCUs), each of which consists of a set of intelligent micro flow-sensor 
nodes and a base station. The intelligent sensor nodes within the same SCU are also 
called sibling nodes. An intelligent micro flow-sensor node consists of a dedicated 
processing unit and an associated micro flow-sensor that measures the variables of 
interest. A base station is needed for each cluster to integrate information from the 
individual sensor nodes in the cluster. Considering the complexity of computation and 
energy constraints of the micro sensor node, a more powerful processor is selected for the 
base station. The base station acts as the control node of the respective SCU, controlling 
the signal transmission of the sensor nodes, managing information rerouting in case of 
link failures, and integrating the information from all the active sensor nodes in the SCU. 
To enable the system to be tolerant to link failures, a sensor node sends information not 
only to the base station, but also to one or more other sibling nodes that are also called its 
backup nodes. Two communication protocols can function among the sibling nodes: (1) 
BP, in which each sensor node broadcasts its information to all its sibling nodes; (2) 
Evenly wide gossiping protocol (EWGP), a revised GP (Liu and Nof, 2001a), in which 
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instead of sending information to only one other sibling node by random selection 
according to GP, the backup nodes are chosen in the way that the roles of backup nodes 
are evenly distributed among the sibling nodes, so it will not happen that some nodes 
assume too many 'duties', while others are 'starved'. 

The two communication protocols, GP and EWGP, can also be used for the 
communications among the base stations. The external processor is the commander of the 
entire network, in which information from all the SCUs is integrated for decision making 
as required. The clusters can be deployed along the flow path, depending on the flow 
situation and applications. 

[Insert figure 3 here] 

3.2 Fault-tolerant Time-out Communication Protocol 

The design of fault-tolerant time-out protocol (FTTP) is motivated by the combination of 
time-out and task coordination protocols (Esfarjani and Nof, 1998; Williams et al., 2002, 
Peralta et al., 2003), and the need for fault-tolerant integration (Liu and Nof, 2001b). At 
the beginning of each data transmission cycle, the base station broadcasts the data 
transmission request to all the sensor nodes in the SCU. Receiving the request, each 
sensor node sends its data to the buffer inside the base station. Each sensor node has a 
corresponding buffer in the base station to store its data so that the base station can trace 
the source of the data. Then the data will be retrieved into the processing unit in the base 
station. According to FTTP, the base station will stop waiting for the information from 
the sensor node if a certain amount of time, J, has passed, where T is the pre-determined 
time-out period and its value, which can also be adaptable, depends on the application. 
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The base station will then announce the rerouting task message to the sibling nodes of the 
node with failed link. Two alternatives exist for where the base station should send the 
task announcement: (1) with a Blind Base Station, the message is sent to all the sibling 
nodes; (2) with a Smart Base Station, the message is only sent to the sibling node that can 
finish the task earliest, based on the current status of the nodes in the system. 

3.2.1 FTTP with Blind Base Station 

In this protocol, when the link of a sensor node, say node v^^ , fails, the base station will 
broadcast the rerouting task message (RTM) for v^to all its sibling nodes. Receiving 
RTM, its backup node bki (the ith backup node of v^^) generates a bid value bvki and sends 
it to the base station, e BKk , where BKk is the set of backup nodes of . bvki can be 
calculated as follows: 

K =('Wa/+1K (1) 

where yia is the transmission rate of node bki (time spent in transmitting one message) and 
mja is the number of tasks that have been offered to node bkh Because of the computation 
limit of the intelligent sensor node, it is unrealistic for it to keep track of its transmission 
rate over time, yki is assumed to be constant for node bju* The base station then evaluates 
the priorities of the backup nodes based on their bid values (the lower the bid value, the 
higher the priority), and sends the task offer message (TOM) to the backup node with the 
highest priority, bk(i). Receiving the task offer message, bk(i) sends the required backup 
data message (BDM) to the base station. 



The structures of the base station and intelligent nodes, and the logic chart of the FTTP 
with blind base station are shown in figure 4 and figure 5 respectively. 

[Insert figure 4 here] 
[Insert figure 5 here] 

3.2.2 FTTP with Smart Base Station 

In this protocol, in order to limit the unnecessary communication and reduce the 
communication traffic, instead of broadcasting RTM for Vj^ , whose link fails, to all its 

sibling nodes, the base station first evaluates the priorities of the its backup nodes and 
then sends TOM only to the node that can finish the task the earliest. The time when the 
backup node bki finishes the task, can be calculated as, assuming the transmission 
starts at time 0: 

^i.=('w^+lK (2) 
where is the number of tasks that have been offered to the backup node bki, and /f^ is 

the transmission rate of bki (time spent in transmitting one message). The base station 
then sends TOM to the node bk(i) whose tk(j) is the smallest. In order to find out fe, the 
base station needs to have a knowledge base which keeps track of the information for 
each node - its backup nodes, transmission rate, and the number of tasks that have been 
offered. Figure 6 shows the structures of base station and intelligent sensor nodes for the 
FTTP with smart base station. The structure of the knowledge base is shown in figure 7, 
in which n is the number of sensor nodes in the cluster and the backup nodes are 
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hypothetical for illustration only. Figure 8 shows the logic chart of the FTTP with smart 
base station. 

Compared with the FTTP with blind base station, the FTTP with smart base station 
requires less communication, and has the advantage of being able to keep track of the 
updated transmission rate of each sensor node, but it requires a base station with more 
complex structure (with an additional knowledge base). Therefore, a tradeoff has to be 
considered in order to determine when to apply each protocol. 

[Insert figure 6 here] 
[Insert figure 7 here] 
[Insert figure 8 here] 

4. TIE (Teamwork Integration Evaluator)/MEMS 

TIE (Teamwork Integration Evaluator)/MEMS is developed in this research as a new 
version of the TIE simulators developed by the PRISM group (Khanna and Nof, 1994; 
Huang and Nof, 1996; Anussornnitisam and Nof, 2000; Anussomnitisam, 2003; Liu and 
Nof, 2001c) at the School of Industrial Engineering in Purdue University. The objective 
of TIE/MEMS is to simulate and evaluate the architectures and communication protocols 
developed for the distributed micro sensor network. TIE/MEMS is programmed with 
MPI (Message Passing Interface). MPI is not a new programming language; rather it is a 
library of subprograms that can be called fi-om C and Fortran 77 programs. Since its 
completion in June 1994, MPI has received widespread acceptance because it is based on 
message passing, one of the powerful and widely used paradigms of programming 
parallel systems (Snir et al., 1996). MPI can be used to simulate different communication 
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schemes among the sensor nodes, such as point-to-point communication, collective 
communication, group communication, and so on (Liu and Nof, 2001c, Liu and Nof, 
200 Id; Huang and Nof, 2002). 

4.1 Structure of TIE/ MEMS 

The overall structure of TIE/MEMS is depicted in figure 9. The inputs of TIE/MEMS 
are: 

(1) micro sensor network architecture, NA\ 

(2) communication protocol, CP; 

(3) sensor network parameters (such as the number of micro sensor nodes, iVy, the number 
of sensor nodes in each cluster, Nc, the number of failed links, etc.); 

(4) constraints unique to the given MEMS sensors, e.g., physical and chemical constraints 
causing micro-stress, proximity noise, and other disturbances to measurement and 
communication. 

TIE/MEMS consists of a collection of libraries, each of which is built to simulate one 
type of sensor network architecture and communication protocol. After the required 
inputs are entered, the corresponding library will be called to execute the program. The 
output of TIE/MEMS is communication cost, C, which is represented by the 
communication time C/, or/and energy consumption Cg. If C is too high and can be 
reduced, NA or CP, or both need to be revised. TIE/MEMS can be extended by adding 
libraries for the newly developed NA and CP. Table 1 lists some examples of constraints 
related to MEMS sensors. 

[Insert figure 9 here] 
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[Insert table 1 here] 

4,2 Simulation Experiments and Results 

Two sets of experiments were conducted to evaluate the sensor network architectures and 
communication protocols described above using TIE/MEMS. The experiments were 
implemented on the SGI Origin 2000 workstation in the Department of Computer 
Science at Purdue University. 

4.2.1 Experiment with Network Architectures and Protocols 

In this experiment set, five combinations of sensor network architectures and 
communication protocols were evaluated. They are HIA & PTP, COA & PTP, COA & 
BP, CLA & BP and CLA & EWGP with the number of backup nodes, Nt, equal to one. A 
total of thirty experiments were run for each of the five combinations and the outputs of 
the experiments are Ct and Ce per communication round. The result of Ct (sample mean) 
vs. Ns is shown in table 2. Table 3 shows the results of Ce vs. iV^ (assuming the 
transmitting cost is 600mv per message transmission and the receivmg cost is 200mw per 
message reception). 

Statistical analysis results (a=95%) show that when Ns is small {Ns =3), there is no 
significant difference of Ct between HIA & PTP and CLA & EWGP, and they have the 
least Ct. There is no significant difference between COA & BP and CLA & BP, and COA 
& PTP has the greatest Ct among the five combinations. When Ns increases, however, Ct 
of the COA network, especially, the COA & PTP, increases dramatically, and CLA & 
EWGP has the least C/. HIA networks have the least Cg, followed by CLA networks, and 
COA networks have the highest Therefore, COA and CLA are preferred in relatively 
small sensor networks because of their fault-tolerant characteristics, and HIA and CLA 
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are favored in large sensor networks due to their cost eflBciency. CLA is a good choice in 
either case. 

[Insert table 2 here] 
[Insert table 3 here] 

4.2.2 Experiments with FTTP 

In this experiment set, FTTP with blind base station and FTTP with smart base station 
were evaluated. Nc and Ni are inputs of the experiment. Thirty experiments were run and 
the output is q per communication round. Table 4 summarizes the experiment results, 
which show, as expected, that the FTTP with blind base station requires significantly 
(a=95%) less evaluating time, Cet, than the FTTP with smart base station. The reason for 
this is that in the FTTP with blind base station, the bid values are calculated by the 
intelligent sensor nodes locally before they are sent to the base station, and the base 
station only evaluates the bid values and chooses the sensor node with the minimal bid 
value. Much more information transmission between the sensor nodes and the base 
station, however, is involved in the FTTP with blind base station. Because the 
information transmission time, C//, is much longer than Cet that requires only computation, 
the total communication time, which is the sum of Cet and of the FTTP with blind 
base station is significantly longer than that of the FTTP with smart base station. The 
difference of C/ becomes more evident when Nc and Ni increase. On the other hand, the 
structure of the smart base station is more complex than that of the blind base station, as 
shown in figure 4 and figure 6. Figure 10 shows the ratio of Ct of the FTTP with blind 
base station to the FTTP with smart base station, rts, vs. Nc when Ni equal to 1, 2 and 3. 
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[Insert table 4 here] 
[Insert figure 10 here] 

5. Conclusions and Contributions of the Research 

MEMS technology is the merger of the integrated circuit (IC) world with the mechanical 
world. It allows the integration of the sensors and actuators with IC at an unprecedented 
level. MEMS flow-sensor can be used in the applications in which conventional sensor 
systems cannot or are not suitable. DMFSA/N, built from a collection of cooperating 
intelligent micro flow-sensor nodes, can improve the reliability and fault-tolerance of the 
system. 

5.1 Array/Network Architecture Design Observation 

The appropriate sensor architecture for DMFSA/N depends on the number of sensor 
nodes needed in the system to provide accurate and timely information. When the number 
of sensor nodes in the system is relatively small, COA and CLA are better architectures 
because of their fault-tolerant characteristics. Otherwise, HIA and CLA are preferred in 
order to limit the communication cost. For example, table 2 shows that when Ns is greater 
than 15, Q of COA & BP is more than 6 times of Ct of HIA & PTP, and Q of COA & 
PTP is more than 30 times of C, of HIA and PTP. The proposed CLA can support 
efficient information communication and integration in the system, because it not only 
requires less communication cost, but also tolerates a certain level of sensor and link 
failures. For example, table 2 shows that when A^^ is greater than 15, G of CLA & BP is 
less than 1/4 of C/ of HIA & PTP, and Q of CLA and EWGP (with equal to 1) is less 
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than 1/8 of Ct of HIA & PTP. In addition, if the flow environment is not steady or non- 
homogeneous along the path, which is very common in the real flow measurement, CLA 
is better than the other alternatives because the clusters can be easily deployed along the 
path. 

5.2 Protocol Design Observation 

The best communication protocols depend on the number of sensor nodes in the system, 
sensor constraints and requirement of response time of the system. Although BP has 
better fault tolerance than EWGP, it requires relatively more communication. For 
example, Table 3 shows that when Ns is equal to 15 {Nc is equal to 5), Q of CLA and BP 
is almost twice Q of CLA and EWGP, and the difference increases when Ns increases. 
Therefore, BP is suitable for the system in which the sensors or other IC components are 
not very robust because it has better fault-tolerance. Compared with the FTTP with blind 
base station, the FTTP with smart base station requires less communication time, but it 
demands the base station to have an additional knowledge base that keeps track of the 
information of the sensor nodes. Therefore, FTTP with smart base station is preferred in 
the on-line control applications, in which response time is critical. 

In order to evaluate the effectiveness of micro sensor network architectures and 
communication protocols, TIE/MEMS has been developed using MPI in this research. 
TIE/MEMS can be expanded by adding libraries for the new developed sensor network 
architectures and communication protocols. 
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This research has been undertaken in response to actual problems encountered in the flow 
measurement area. Given the widespread use and importance of flow-sensors, the 
limitations of the current flow measurement technology are surprising. DMFSA/N opens 
original new directions for this area. Some of the main contributions of this research are: 

1. Intelligent sensors manufactured using advanced IC technology will 
fundamentally change the nature of sensing and control systems, and modeling 
computational tool, such as TIE/MEMS, to evaluate their performance are useful, 

2. The ability to place and activate the micro flow-sensors where and when they are 
needed enables better control of the flow process. 

3. Redundancy enables more reliable monitoring and control because occurrence of 
failures of some sensors will not inhibit the sensing ability of the system. 

For contributions 2 and 3, a FTTP is necessary to enable timely integration of the 
distributed sensors' readings. For future research, TIE/MEMS is being further developed. 
MEMS flow-sensor physical constraints (as shown in Table 1) are being quantified and 
incorporated into the simulation evaluation of the sensor network architectures and 
communication protocols. 
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Table 1 . Examples of Constraints Related to MEMS Sensors 







Scaling issues 


Viscosity force; Surface effect 


Problems in 
micromachining 


Squeeze film damping; Particular contamination; Static stiction 


Response 


Decrease of response time due to protective gels, coatings or 
diaphragms 


Reliability of 
microelectronic device 


Breakdown of a wire-bond to a junction; Defective encapsulation; 
Defective semiconductor material; Thermal mnaway 


Life time 


Aging; Contamination of the device; Limited power supply 


Noise/interferences 


Degraded signal due to the variance of temperature, humidity, etc. 



Table 2. Communication Time (ms) of Five Combinations of Network Architectures, NA, 
and Communication Protocols, CP, vs. No. of Sensor Nodes, Ns 



NA 


CP 




3 


8 


15 


20 


30 


HIA 


PTP 


0.109 


0.441 


0.782 


0.902 


1.674 


COA 


PTP 


0.725 


6.153 


22.920 


40.894 


93.200 


COA 


BP 


0.471 


2.233 


4.835 


9.217 


20.555 


CLA 


BP 


0.438^ 


0.161^ 


0.192' 


0.205* 


0.242' 


CLA 


EWGP (JVi=1) 


0.112' 


0.118^ 


0.115' 


0.115* 


0.121' 



'Nc=3; ^Nc=4; ^Nc=5; '*Nc=5; ^Nc=6 



Table 3. Energy Consumption (mw) of Five Combinations of Network Architectures, 
NA, and Communication Protocols, CP, vs. No. of Sensor Nodes, Ns 



NA 


CP 


Ns 


3 


8 


15 


20 


30 


HIA 


PTP 


600* 


600* 


600* 


600* 


600* 


N/A 


800** 


1200** 


1400** 


1400** 


400*** 


800*** 


800*** 


800*** 


1600*** 


COA 


PTP 


4800 


44800 


168000 


304000 


696000 


COA 


BP 


2400 


15400 


50400 


87400 


191400 


CLA 


BP 


1000^ 


1200^ 


1400^ 


1400" 


1600' 


CLA 


EWGP (A^A=1) 


800^ 


800^ 


800' 


800" 


800* 



(* Ce for the leaf nodes ; ** Cg for the intermediate nodes ; *** Ce for the root node; no 
intermediate nodes for HIA with when Ns is equal to 3; ' Nc=3; ^Nc=4; ^Nc=5; '*Nc=5; ' 
Nc=6) 
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Table 4. Communication Time (|lis), Ct, of the FTTP with Blind Base Station and FTTP 
with Smart Base Station vs. No. of Sensor Nodes in the Cluster, iVc, and No. of 
Failed Links, Ni 



CTTD 

r 1 1 r 


A/ 


hi 


Cet 










4 


2.17 


142.77 


144.93 




1 


8 


2.70 


261.13 


263.83 




15 


2.90 


477.37 


480.27 






21 


3.93 


773.07 


777.00 






4 


3.13 


346.60 


349.73 


FTTP with Blind 


2 


8 


4.20 


581.90 


586.10 


Base Station 


15 


5.87 


1146.50 


1152.37 






21 


7.40 


1706.13 


1713.53 






4 


4.63 


556.80 


561.43 




3 


8 


5.65 


827.13 


832.77 




15 


7.70 


1721.53 


1729.23 






21 


10.45 


2587.52 


2597.97 






4 


4.33 


22.43 


26.77 




1 


8 


5.64 


28.91 


34.55 




15 


7.16 


33.06 


40.23 






21 


8.17 


40.49 


48.66 






4 


7.90 


45.33 


53.23 


FTTP with Smart 


2 


8 


9.89 


54.48 


64.37 


Base Station 


15 


13.90 


58.50 


72.40 






21 


18.46 


60.68 


79.14 






4 


10.60 


62.13 


72.73 




3 


8 


13.30 


66.96 


80.26 




15 


16.83 


70.69 


87.51 






21 


23.86 


72.91 


96.77 
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Figures: 

Figure 1. Sensor Network Architectures: (a) Committee Architecture (b) Hierarchical 

Architecture (c) Flat Tree Architecture (d) Multilevel Binary de Bruijn Architecture 

Figure 2. Communication Protocols: (a) Point-to-Point Protocol (b) Broadcasting 
Protocol (c) Gossiping Protocol 

Figure 3.New Cluster Architecture for DMFSA/DMFSN (dash line means the link may 
or may not exist depending on the communication protocol used) 

Figure 4. Structures of the Base Station and Intelligent Sensor Nodes for the FTTP with 
Blind Base Station (dash line means the link may or may not exist depending on 
the protocol used.) 

Figure 5. Logic Chart of the FTTP with Blind Base Station 

Figure 6. Structures of the Base Station and Intelligent Sensor Nodes for the FTTP with 
Smart Base Station (dash line means the link may or may not exist depending 
on the protocol used.) 

Figure 7. Logic Chart of the FTTP with Smart Base Station 

Figure 8. Logic Chart of the FTTP with Smart Base Station 

Figure 9. Structure of TIE/MEMS (dash line- under development) 

Figure 10. Ratio of Communication Time of the FTTP with Blind Base Station to the 

FTTP with Smart Base Station, r^j, vs. No. of Sensor Nodes in the Cluster, A/c, 
when the No. of Failed Links, M, Equal to 1, 2 and 3 
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Figure 1. Sensor Network Architectures: (a) Committee Architecture (b) Hierarchical 

Architecture (c) Flat Tree Architecture (d) Multilevel Binary de Bruijn Architecture 
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Figure 4. Structures of the Base Station and Intelligent Sensor Nodes for the FTTP with 
Blind Base Station (dash line means the link may or may not exist depending on 

the protocol used.) 
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Figure 5. Logic Chart of the FTTP with Blind Base Station 
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Figure 6. Structures of the Base Station and Intelligent Sensor Nodes for the FTTP with 
Smart Base Station (dash line means the link may or may not exist depending 
on the protocol used.) 
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Abstract: Distributed micro flow-sensor arrays and networks (DMFSA/N), built from collections of spatially 
scattered, cooperating intelligent micro flow-sensor nodes, can improve the accuracy and reliability of system. 
However, it is unrealistic to expect all the sensor nodes and communication links in the system to function properly 
all the time. A new fault-tolerant time-out Protocol (FTTP) with two alternatives - one with blind station, and the 
other with smart base station, and a fault-tolerant sensor integration algorithm (FTSIA) were developed and 
implemented in this research. TIE (Teamwork Integration Evaluator)/MEMS was developed to evaluate the 
communication protocols for the distributed micro sensor network. Experiment results show that although FTTP 
with smart base station requires more complex structures, compared with FTTP with blind base station, it involves 
less communication. Experiment results also show that FTSIA can always give reliable results even when some 
sensors provide faulty information, and the results from FTSIA are significantly more accurate than the statistical 
mean if some sensors give faulty readings. 

Keywords: MEMS sensors, time-out, teamwork integration evaluator, sensor integration, communication protocol 
1. Problem Background 

Flow-sensors are used to measure the movement of fluid flow (liquid flow or gas flow). Conventionally, single 
sensor systems have been used to save space and manufacturing cost. But they are not suitable for critical 
applications because a single sensor is inevitably affected by noise, refresh delay and other uncertainty issues and 
cannot guarantee to deliver accurate information all the time. Development of MEMS (Micro Electro Mechanical 
Systems) has enabled production of large amounts of micro sensors at low cost. Micro flow-sensors are one of tiie 
most common MEMS devices which have been developed in recent years [1, 2, 3]. Distributed micro flow-sensor 
arrays and networks (DMFSA/N), built from collections of spatially scattered, cooperating intelligent micro flow- 
sensor nodes, can improve the accuracy, timeliness, and reliability of systems. Sensors incorporated with dedicated 
signal processing functions are called intelligent sensors, or smart sensors. The main roles of dedicated signal 
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processing are to enhance design flexibility and realize new sensing functions. Additional roles are to reduce loads 
on central processing units and signal transmission lines by distributing information processing to the lower layers 
of the system [4], Here, the micro flow-sensor array refers to a group of the same type of micro flow-sensors which 
are placed close together to measure the same variable of interest. Another possibility is to combine hybrid sensors, 
measuring different variables, in the same array [5, 6]. Different groups of micro flow-sensor arrays can be 
distributed widely and selectively in the flow environment to form a micro flow-sensor network. 
Due to the inherent constraints of sensors and other electronic components, and the effects of flow environment, 
however, it is unrealistic to expect sensor nodes and the communication links to function properly all the time. The 
objective of this paper is to address the design of fault-tolerant communication protocol and information integration 
algorithm for DMFSA/N [7]. A system is said to be fault-tolerant if its performance is not affected by faults in the 
system, such as sensors measure the variable of interest not accurately or timely, communication links break down, 
etc. A fault-tolerant time-out Protocol (FTTP) with two alternatives - one with blind stations, and the other with 
smart base stations, and a fault-tolerant sensor integration algorithm (FTSIA) have been proposed and implemented 
in this research. TIE (Teamwork Integration Evaluator)/MEMS has been developed to evaluate the communication 
protocols for the distributed micro sensor network. 

This paper is organized as follows. In section II, the cluster architecture proposed by Liu and Nof [8], which can 
support fault-tolerant communication and integration for DMFSA/N, is reviewed. The new FTTP with blind base 
station and FTTP with smart base station are described in Section III. In section IV, TIE/MEMS is described. Two 
alternatives of FTTP are evaluated using TIE/MEMS. In section V, the proposed fault-tolerant sensor integration 
algorithm (FTSIA) is presented. Three cases of simulation experiments are conducted to test FTSIA and the results 
from FTSIA are compared with the statistical mean. A case study of measuring the airflow using a pressure sensor 
array is described in section VI to illustrate the application of FTSIA in real flow measurement. The experiment 
setup and results are also described. Finally, conclusions and future challenges are given in section VII. 

2. A Cluster Architecture for DMFSA/N (CLA) 

Since micro flow-sensor nodes are distributed spatially in DMFSA/N, a suitable architecture has to be designed to 
support the efficient and fault-tolerant communication and integration in the system. CLA, a cluster architecture for 
DMFSA/N, was proposed and developed by Liu and Nof [8]. In this structure, DMFSA/N are divided into sensor 
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cluster units (SCUs), each of which consists of a set of intelligent micro flow-sensor nodes and a base station. The 
intelligent sensor nodes within the same SCU are also called sibling nodes. An intelligent micro flow-sensor node 
consists of a dedicated processing unit and an associated micro flow-sensor that measures the variable of interest. A 
more powerful processing unit is used as the base station. The base station acts as the control node of the respective 
SCU, controlling the signal transmission of the sensor node, managing information rerouting in case of link failures 
and integrating the information from all the sensor nodes in the SCU. In order to design the system with tolerance to 
link failures, the sensor node sends information not only to the base station, but also to one or more other sibling 
nodes that can be considered as its backup nodes. Two communication protocols can be used among the sibling 
nodes: (1) Broadcasting Protocol (BP), in which each sensor node broadcasts its information to all its sibling nodes; 
(2) Evenly wide gossiping protocol (EWGP), a revised gossiping protocol (GP), in which instead of sending 
information to only one other sibling node by random selection according to GP [9], the backup nodes are chosen in 
a way that the roles of backup nodes are evenly distributed among the sibling nodes. The advantage of even 
distribution is to avoid a situation where some nodes are assigned too many 'duties', while others are 'starved'. 
These communication protocols can also be applied for the communication among base stations. The external 
processor is the commander of the entire network, in which information from all the SCUs is integrated for decision 
making as required. The clusters can be deployed along the flow path, depending on the particular flow situation and 
applications (figure 1). Comparison of the CLA and other architectures, including committee, hierarchical, flat tree, 
and multilevel binary de Bruijn architectures is described in detail in Liu and Nof [8]. 
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Figure 1. Proposed Cluster Architecture for DMFSA/DMFSN (dash line means the link may or may not exist 
depending on the communication protocol used) 

III. Fault-tolerant Time-out Communication Protocol 

The design of fault-tolerant time-out protocol (FTTP) is motivated by the combination of time-out and task 
coordination protocols [10, 11,12, 13], and the need for fault-tolerant integration [14]. At the beginning of each data 
transmission cycle, the base station broadcasts the data transmission request to all the sensor nodes in the SCU. 
Receiving the request, each sensor node sends its data to the buffer inside the base station. Each sensor node has a 
corresponding buffer in the base station to store its data so that the base station can trace the source of the data. Then 
the data will be retrieved into the processing unit in the base station. According to FTTP logic, the base station will 
stop waiting for the information from a sensor node if a certain amount of time, T, has passed, where 7 is a pre- 
determined time-out period and its value which can be adaptable depends on the application. The base station will 
then announce the rerouting task message to its sibling nodes. Two alternatives exist for where the base station 
should send the rerouting task announcement: (1) with a Blind Base Station, the message is sent to all the sibling 
nodes; (2) with a Smart Base Station, the message is only sent to the sibling node that can finish the task earliest, 
based on the current status of the nodes in the system. 

3.1 FTTP with Blind Base Station 
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In this protocol, when the link of a sensor node, say node v^^ , fails, the base station will broadcast the rerouting task 
message (RTM) for V/^ to all its sibling nodes. Receiving RTM, its backup node bti (the iih backup node of ) 
generates a bid value bvki and sends it to the base station, bu £ BKk,f where BKk is the set of backup nodes of . 
bvki can be calculated as follows: 



(1) 



where yki is the transmission rate of node bid (time spent in transmitting one message) and njki is the number of tasks 
that have been offered to node bkh Because it is unrealistic for the intelligent sensor node to keep track of its 
transmission rate over time due to its computation limit, yk, is assumed to be constant for node bkh The base station 
then evaluates the priorities of the backup nodes based on their bid values (the lower the bid value, the higher the 
priority), and sends the task ofTer message (TOM) to the backup node with the highest priority, bk(i)' Receiving the 
task offer message, bk(i) sends the required backup data message (BDM) to the base station. The structures of the 
base station and intelligent nodes, and the logic chart of the FTTP with blind base station are shown in figures 2 and 
3 respectively. 




Intelligent 
Sensor Node 



Figure 2. Structures of the Base Station and Intelligent Sensor Nodes for the FTTP with Blind Base Station (Dash 
line means the link may or may not exist depending on the protocol used.) 
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Figure 3. Logic Chart of the FTTP with Blind Base Station 
3.2 FTTP with Smart Base Station 

In this protocol, in order to limit the unnecessary communication and reduce the communication traffic, instead of 
broadcasting RTM for v^^ , whose link fails, to alt its sibling nodes, the base station first evaluates the priorities of 

the its backup nodes and then sends TOM only to the node that can finish the task the earliest. The time when the 
backup node bid finishes the task, tki, can be calculated as (assuming the transmission starts at time 0): 

/^=(/w,,+lK (2) 
where is the number of tasks that have been offered to the backup node 6^ and y^^ is the transmission rate of 

bid (time spent in transmitting one message). The base station then sends TOM to the node bk(\) whose t^^y^ is the 

smallest. In order to find out , the base station needs to have a knowledge base which keeps track of the 
information for each node - its backup nodes, transmission rate, and the number of tasks that have been offered. 
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Figure 4 shows the structures of base station and intelligent sensor nodes for the FTTP with smart base station. The 
mode! of the knowledge base is shown in figure 5, in which n is the number of sensor nodes in the cluster and the 
backup nodes are hypothetical for illustration only. Figure 6 shows the logic chart of the FTTP witii smart base 
station. 

Compared with the FTTP with blind base station, the FTTP with smart base station requires less communication, 
and has the advantage of being able to keep track of the updated transmission rate of each sensor node, but it 
requires a base station with more complex structure (with an additional knowledge base). Therefore, a tradeoff has 
to be considered in order to determine when to apply each protocol. 
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Figure 4. Structures of the Base Station and Intelligent Sensor Nodes for the FTTP with Smart Base Station (Dash 
line means the link may or may not exist depending on the protocol used.) 
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Figure 5. Logic Chart of the FTTP with Smart Base Station 
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Figure 6. Logic Chart of the FTTP with Smart Base Station 

4. TIE (Teamwork Integration £valuator)/ME]V[S 

TIE (Teamwork Integration Evaluator)/MEMS has been developed in this research as a new version of the TIE 
simulators developed by the PRISM group [15, 16, 17, 18] in School of Industrial Engineering at Purdue University. 
The objective of TIE/MEMS is to simulate and evaluate the architectures and communication protocols developed 
for the distributed micro sensor network. 

The inputs of TIE/MEMS include micro sensor network architecture, NA, communication protocol, CP, sensor 
network parameters (such as the number of micro sensor nodes, N^, the number of sensor nodes in each cluster, Nc, 
the number of failed links, N/, etc.), and constraints unique to the given MEMS sensors, e.g., physical and chemical 
constraints causing micro-stress, proximity noise, and other disturbances to measurement and communication. The 
output of TIE/MEMS is the communication cost, C, which is represented by the communication time c,, or/and 
energy consumption Ce. If C is too high and can be reduced, NA or CP, or both need to be revised. TIE/MEMS can 
be extended by adding libraries for the newly developed NA and CP. More details on TIE/MEMS can be found in 
Liu andNof[18]. 

4.1 Simulation Experiments with FTTP 

Simulation experiments were conducted to evaluate the FTTP with blind base station and FTTP with smart base 
station using TIE/MEMS. The experiments were implemented on the SGI Origin 2000 workstation at the 
Department of Computer Science in Purdue University. Nc and Ni are inputs of the experiment. Thirty experiments 
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were run and the output is c, per communication round. Table 1 summarizes the experiment results, which show, as 
expected, that the FTTP with blind base station requires less evaluating time, Ceh than the FTTP with smart base 
station. The reason for this is that in the FTTP with blind base station, the bid values are calculated by the intelligent 
sensor nodes locally before they are sent to the base station, and the base station only evaluates the bid values and 
chooses the sensor node with the minimal bid value. Much more information transmission between the sensor nodes 
and the base station, however, is involved in the FTTP with blind base station. Because the information transmission 
time, C//, is much longer than c^/ that requires only computation, the total communication time, C/, which is the sum 
of Cet and c,/, of the FTTP with blind base station is longer than that of the FTTP with smart base station. The 
difference of C/ becomes more evident when Nc and Ni increase. On the other hand, the structure of the smart base 
station is more complex than that of the blind base station, as shown in figure 2 and figure 4. Figure 7 shows the 
ratio of c, of the FTTP with blind base station to the FTTP with smart base station, r^s, vs. Nc when Ni equal to 1, 2 
and 3. 



Table 1. Communication Time, C/, of the FTTP with Blind Base Station and FTTP with Smart Base Station vs. No. 
of Sensor Nodes in the Cluster, N^ and No. of Failed Links, N 
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Figure 8. Ratio of Communication Time of the FTTP with Blind Base Station to the FTTP with Smart Base Station, 
rts, vs. No. of Sensor Nodes in the Cluster, A/^, when the No. of Failed Links, Equal to 1, 2 and 3 



5. Fault-tolerant Sensor Integration Algorithm 

Fault-tolerant sensor integration algorithm (FTSIA) refers to the algorithm used to combine information from 
different sensors in the system and produce reliable results even if some sensors yield faulty information. The 
algorithm developed in this research follows the ideas originated from Marzullo [19] and later improved by 
Jayasimha [20]. It deals with the competitive integration of sensor information, in which each sensor ideally 
measures identical information, but in reality, is subject to noise and other disruptions. This algorithm does not deal 
with the situation when sensors do not respond at all. This case is handled by the FTTP functions described earlier. 
First, the logic of FTSIA is described, followed by experiments demonstrating FTSIA advantages, and then 
explaining how FTSIA can be combined with FTTP to gain fault-tolerance advantages under time-out integration 
regime. 



5.1 Some Definitions 

A concrete sensor (CS) is a device that measures the variable of interest in the environment. An intelligent micro 
flow-sensor node is an example of CS. An abstract sensor (AS) is represented as a dense interval that contains the 
value of the physical variable read by the CS. The reason for using AS is that it is more reasonable to consider the 
reading from a sensor as a continuous set of values instead of a point value, because of the variation of the data due 
to noise, manufacturing, and other uncertainty issues. AS is derived from the point value received from CS. If a CS, 
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say , reads a value, say r,. , and its maximum tolerable variation can be ± ff^ , then AS of can be defined as 

the interval [L^ ? ] = - ^/ 5 ^/ + ^, ] » where L/ is the lower bound of v,. , and Ui is the higher bound of v,- . 
An AS is correct if its interval contains the true value of the physical variable being measured. Otherwise, it is 
faulty. If AS /4/ and A 2 are correct, they must intersect and the intersection contains the true value of physical 
variable. An AS is said to be tamely faulty if it is not correct but its interval overlaps with that of a correct AS. A 
faulty AS that is not tamely faulty is said to be widely faulty. 
A distinct m-interval is an interval in which no more than m intervals intersect. 

5.2. The Logic of Algorithm 

Input: Readings from Ns CSs, with at most / (either tamely faulty or widely faulty of them are faulty. 

Output: A single output interval and a list of all the possibly tamely faulty sensors and surely faulty sensors (either 

tamely faulty or widely faulty). 

Assumption: f <jN^ 

Pseudo Program of FTSIA: 

1. Import readings from A^^ CSs and generate the respective ASs; 
LB -< LB^.LB^ LB , LB ^ > , where LB,=t,-€,\ 

MB =< HB^,HB^,,.,, MB,,,.,, HB > , where HB , = + f , ; 

2. Sort the ASs by their lower bounds, and this sorted list is i =< , ^2 Lff > and the corresponding upper 
boundsetis j/ * =< C/,',^// f/^ ' > ; 

3. Sort the first {Ns -/-I) intervals of U and this sorted list is U ; 

4. Initialize the variables 

/ := O ; /* Initialize the interval set, in which each interval has three attributes, lower bound /, higher bound u 

and level of interval m */ 
last := 0 ; /* Initialize the number of intervals in / */ 

5. Main part of the algorithm 
for j := / down to 0 do 

5.1 Insert u ^ .'in its appropriate place in U to keep V sorted ; 

5.2 Find the intersection oi{N ^ - f) intervals 

if (L„ ,^t/,„),then /* j/^, is the first element in (/ ♦/ 
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5.3 Insert TV := [L ^ _y , t/ ] ,which is a ( A^^ - / )-interval, into I and make appropriate changes to I 
5.3.1Locate the connected interval sequence 5 = [C,,C2,..., CJ in / in which each C/ intersects with ^ 

5.3.2 Condition 1 : If S is empty, insert to the tail / 
if (5 = O ) 

I := N ; 
last := 1 ; 

else 

5.3.3 Condition 2: If there is only one interval C, is S, find the intersection of N and Ci 
if(5 = [C,] ) 

else /*if (t/^,) >«,,)*/ 

/ := (/ u N, u A^2)- C, ; 

to/ := last + 1 ; 
end if 

else /Mf (/^^ <I^,-y)*/ 

Wat, := iV, - / ; 

to/ := te/ + 1 ; 
else /*if t/^,) > w,^ */ 
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last := te/ + 2 ; 
end if 
end if 
end if 

else 

5.3.4. Condition 3: intervals C2, Cp.i completely overlaps N 
for ifc := 2 to p - 1 do 

end for 

if( =1^ _.)then 

else /* if /^^ <L^ _j , Split [/^^ , u^^ ] into two intervals */ 
'"at, ^c, ' 

'"yv, •= 'w,, +1; 

/ := (/ u TV, u A^2)- C, ; 

ter := last + 1 ; 
end if 

if{(/(,j < w,^ )then 



/:=:(/uA^3u7V,)-C,; 

last := last + 1 ; 
else /♦ if (7^,^ > u,^ */ 

if(i/^,j > i/,^ )then 

'Wat, •= - /J 

/:=(/u JV3)-C,; 

end if 
end if 
end if 

5.3.5 Remove the AS [ L^j^ , ] from the list of ASs 

Z := I - L^,^',U = U - 
end for 

6. Find out all the possibly tamely faulty and surely faulty sensors 
for / := 0 up to last - 1 do 
for j := 0 up to N ^ do 

ifiiHBj >- ) or {lb J <= )) /* sensory is faulty if interval i contains the true value */ 
Ff = Ff Kj j \ /* insert sensor j into the set of faulty sensors of interval i; F; is the set of faulty sensors of 
interval i */ 

end if 
end for 
end for 

fori:=Oupto 

forj:=Otolast-l do 
if (/eF^ then 

:= rt^ + 1 ; 

end if 
end for 
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sensor i surely correct; 
else 

= last ) 

sensor i is surely faulty; 
else 

sensor i is possibly tamely faulty; 
end if 

CLB := LB - {LB .}\ /*CLB is set of the lower bound of the surely correct sensors */ 
CHB := HB - {HB /*CHB is set of the higher bound of the surely correct sensors ♦/ 
end if 
end for 

7. Find out the output data 

sort CLB and the sorted list is CL =< CL^,CL CL ^ >\ 
sort CHB and the sorted list is QH =< CH ^,CH ^ ,'-CH „ > ; 
= [CjL „ , C// , ]; /* lo is the output interval */ 



5.3 Simulation Experiments and Results 

Three experiment cases were simulated to test FTSIA. The mean from FTSIA, which is the average of the higher 
bound and lower bound of the output interval of FTSIA, were compared with the statistical mean of the same 
measurement values. The results from the FTSIA and statistical calculations are also affected by the distribution of 
sensor readings. In order to evaluate the reliability of both methods, the worst situation was simulated in all three 
cases as described below. All the sensor nodes generate readings above the true value of the measured physical 
variable, i.e., Ar,. =r^-r*>0 (/ = 1,2,..., A^J, where r* is the true value of the measured physical 

variable, is the reading received from the sensor node v,. , and Ar, is the error of the reading from v,. . The 
sensor readings were assumed to follow the uniform distribution. In all three cases, assuming r* = 0.5 , and the 
maximum tolerable variation of the sensor reading, 5* , is the same for all the sensor nodes and € = ±0.005 . Shift of 
values from 0.5025 to 0.505 needs to be detected because the correct readings should not be greater than 0.505. 
1) Case 1. In this case, all the sensors are assumed be correct, generating readings with errors within ±1% of r* and 
their ASs contain r*. Five experiment sets were conducted in this case: Ns is equal to 8, 15, 21, 27 and 36. 
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Assuming both the type I risk and type II risk are 0.025, the sample size required to detect the shift of values firom 
0.5025 to 0.505 is no less than seven. 

2) Case 2. In this case, some sensors are assumed to be tamely faulty, i.e., they generate readings with errors beyond 
±1% but within ±2% of r* . Since the tamely faulty ASs intersect with some correct ASs, they can form (N^-f)- 

intervals that do not contain the true value. Because it is impossible to know in advance which interval is the correct 
interval, according to this algorithm, all the possibly faulty ASs are removed, and only those surely correct ASs will 
be kept. Following this logic, however, some correct ASs will also be considered as possibly faulty and thus 
removed. In order for the algorithm to be effective in this case,/ must satisfy that < 7 A^, , because there will be 

no correct ASs left after removing all the possibly faulty ASs when >^N^, Twelve experiment sets were 

conducted in this case: N,. is equal to 8 with / equal to 1, 2 and 3; Ns is equal to 15 with f, equal to 1, 3, 5 and 7; and 
Ns is equal to 21 with ft equal to I, 3, 5, 7 and 10. Assuming both the type I risk and type II risk are 0.025, the 
sample size required to detect the shift of value from 0.5025 to 0,505 is no less than seven. 

3) Case 3. In this case, some sensors are assumed to be widely faulty, i.e., they generate readings with errors beyond 
±2% of r*. In order to control the experiments, however, their readings are assumed to be within ±10% of r*. 
Because when the readings from the widely faulty sensors are very close, it is possible that the widely faulty ASs 
intersect and form ( - f )-intervals that deviate the results from the correct one. In order to prevent this 

situation, the number of widely faulty sensors,/^, needs to satisfy <2^s' Twelve experiment sets were 

conducted in this case: A^^ is equal to 8 with fy, equal to 1, 2 and 3; A^Js equal to 15 with equal to 1, 3, 5 and 7; and 
Ns is equal to 21 with equal to 1, 3, 5, 7 and 10. Assuming both the type I risk and type II risk are 0.025, the 
sample size required to detect the shift of value from 0.5025 to 0.505 is no less than one. 

Each experiment set was conducted twelve times. Figure 9 shows the result from FTSIA and the statistical mean vs. 
Ns when/=0. The results from FTSIA and statistical mean vs.^I when Ns is equal to 8 are shown in figure 10. The 
results from FTSIA and statistical mean vs.jC when A^^ is equal to 8 are shown in figure 1 1. 
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Figure 8. Ratio of Communication Time of the FTTP with Blind Base Station to the FTTP with Smart Base Station, 
rts, vs. No. of Sensor Nodes in the Cluster, Nc^ when the No. of Failed Links, Equal to 1 , 2 and 3 
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Figure 9.Result from FTSIA and Statistical Mean vs. Number of Sensor Nodes (r*=0.5, £=±0.005) 
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Figure 10. Result from FTSIA and Statistical Mean vs. Number of Sensor Nodes (r*=0.5, «=±0.005, JV,=8) 



S.4 Summary and Conclusions of Experiments 
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Analysis of Variance (ANOVA) was used to analyze the results from the experiments in section 5.3 (a=95%). The 
mean from FTSIA and statistical mean were compared. Table 2 summarizes the analysis results, which show that 
when /=0 or / is small relative to Ns, there is no significant difference between the mean from FTSIA and the 
statistical mean. However, when ft increases or/w>0, the mean from FTSIA is significantly more accurate than the 
statistical mean. Table 3 shows the effects of N, and / on the mean from FTSIA in case 2. The reason for such 
effects is that tamely faulty ASs can intersect with some correct ASs and form ( A^, - / )-intervals that do not 
contain the true value, and the correct ASs that do not intersect with the inconect (N^ -f )-intervals will also be 
removed as possibly faulty ASs. The correct ASs left usually have relatively higher errors because they not only 
intersect with the correct ASs, but also with some tamely faulty ASs. For example, in case 2 shown in section 5.3, 
when Ns^lS and f, =7, and when M=21 and f, =10, only one correct ASs is left after removing all the possibly faulty 
ASs and this sensor has the highest sensor reading, i.e., the largest error, although still within the tolerable limits. 
Therefore, the performance of FTSIA degrades gracefiiUy when / increases, but the output interval it generates still 
contains the true value. Table 3 also indicates that in order to enable FTSIA generate more accurate results, the 
number of correct sensors left after deletion, «c, should be at least half of A/,, i.e., n^>jN^, «ccan be calculated as 



Therefore, when all the faulty sensors are tamely faulty in the system, in order to generate more accurate result using 
the FTSIA, iV^ should satisfy 



For example, when f,=\ and A^, > 4, when / =3 and A^, >12, and when/r=5 and Ns ^ 20, the FTSIA will yield more 
accurate results. 

The widely faulty ASs (shown in case 3), however, do not cause the same problem as tamely friulty ASs. They do 
not intersect with the correct ASs and thus are easier to be detected and removed, without affecting the correct ASs 
(see table 4). Therefore, /w do not have significant effects on the results from FTSIA, as long as the/n. is less than 
half of A^j, i.e., < jN^ . Therefore, when all the faulty sensors are widely faulty in the system, A^^ should 



(3) 



(4) 



satisfy 



w 



(5) 
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Since equations (4) and (5) are linearly independent, it can be concluded tiiat when both tamely faulty and widely 
faulty sensors exist in the system, in order for FTSIA to generate more accurate results, should satisfy 

iV,>4*/+2*A (6) 
The outputs of FTSIA also indicate which sensor nodes are possibly tamely faulty, and which are surely faulty 
(tamely or widely). If some sensor node exists in the all the possibly faulty sensor sets, it is surely faulty. On the 
other hand, if some sensor node appears in some faulty sensor sets, but not all of them, it is possibly tamely faulty. 

5.5 FTSIA Algorithmic Characteristics 

The simulation experiment results show that the proposed FTSIA is effective because although in the simulated 
worst situation, when f, is large relative to N„ the FTSIA may degrade gracefully, it always generates the output 
interval that contains the true value. This algorithm has a complexity of 0(Ns\og(Ns) because of the sorting in step 2 
and 3 shown in the logic of the algorithm (the complexity of computing statistical mean is of 0(Ns)), This algorithm 
is simple and easy to be implemented. 

The algorithm proposed by [19] only provides a single interval from all the sensors without fault detection. The 
algorithm proposed by Jayasimha [20] provides all the (A';j-i)-interval (0<i<J)y but in the flow measurement 
applications, single output results are preferred. The improvement of the FTSIA proposed in this thesis over the two 
algorithms described above is that it not only detects the possibly faulty sensors and surely faulty sensors but also 
generates a single output interval from the surely correct sensors after removing those faulty sensors. 
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Table 2. ANOVA Results (a=95%) of Effects of No. of Total Sensor Nodes and No. of Faulty Sensors on the Mean 
from FTSIA and Statistical Mean (From three experimental cases) 
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Table 3. ANOVA Results (a=95%) of Effects of No. of Total Sensor Nodes and No. of Tamely Faulty Sensors on 
the Mean from FTSIA (From Case 2) 
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Table 4. ANOVA Results (a=95%) of Effects of No. of Total Sensor Nodes and No. of Widely Faulty Sensors on 

the Mean from FTSIA (From Case 3) 
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6. Case Study 

This section illustrates the application of the FTSIA described in section V through a case study - measuring the 
flow using a pressure sensor array. 

6.1 Decision on the number of sensors, Ns 

Before the experiment was conducted, Ns had to be decided to generate more accurate result using the FTSIA. The 
initial experiments showed that among the 8 available pressure sensors, / is equal to 2. Therefore, according to 
equation (6), all 8 available pressure sensors had to be used to form a pressure sensor array in the experiment. 

6.2 Description of the experiment setup 

This experiment was conducted in the PRISM lab at Purdue University. The experiment setup is shown in figure 12. 
It consists of a vacuum cleaner, which was used to suck air to go through the pipe connected to it, and a pressure 
sensor array, which had 8 pressure sensors. 8 holes were drilled on the pipe and signals were then input into the 
pressure sensors through the 1/8 Tygon Tubing. Because the input signals were very low (less than 0.0 Iv), 
differential inputs were used to reduce the noise effect. Data Acquisition board (DAQ) was used to change the 
analog signals into the digital signals which were input into the computer for data analysis. In this experiment, the 
DAQ used was KPCI-3107, manufactured by Keithley Instruments Inc. KPCI-3107 which is compatible with 
LabView, etc. Finally, the data retrieved from LabView were input into the FTSIA to generate the output interval. 
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Figure 12. Experiment Setup of Flow Measurement with a Pressure Sensor Array 



6.3 Experiment results 

Table 5 shows an example of the data, t},, received from the pressure sensor array and the lower bounds, Li, and 
higher bounds, Uj, of the data intervals, assuming £=±0.0005. After the Li and Ui are input into the FSTIA, the 
outputs generated are as shown in table 6. The output results show that there is an 8-interval, which means all ASs 
intersect, two 7-intervals and one 6- interval. When there are tamely faulty sensors in the system, and if the lower 
bounds (higher bounds) of their ASs are lower (higher) than the maximum (minimum) of the higher bounds (lower 
bounds) of the correct ASs, it is possible that all ASs intersect, so it cannot be concluded tiiat there is no faulty 
sensor in the system. Therefore, all three possibly tamely faulty ASs, 4, 5, and 7, have to be removed, and the output 
interval (the intersection of all the surely correct ASs) is [0.006747, 0.007685]. 
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Table 5. An Example of the Data Received from 8 Pressure Sensors (volt) and the Lower Bounds and Higher 
Bounds of the Data Interval (Assuming £=±0.0005) 
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Table 6. Summary of the Outputs from FTSIA (Data from table Sare inputs) 
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Calculation 


Mean 
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Standard Deviation 
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7. Conclusions and Future Tasks 

This paper proposes a fault-tolerant time-out protocol (FTTP) with two alternatives and a fault-tolerant sensor 
integration algorithm (FTSIA) to address the fault-tolerance issues in DMFSA/N. TIE/MEMS was developed to 
evaluate the two alternatives of FTTP. Simulation experiments show that the FTTP with smart base station requires 
less communication time than FTTP with blind base station, but it demands the base station to have an additional 
knowledge base that keeps track of the information of the sensor nodes. Therefore, FTTP with smart base station is 
preferred in the on-line control applications, in which the response time is critical. Three experiment cases were 
simulated to test the reliability of the proposed FTSIA and the results from FTSIA were compared with the 
statistical mean. The experiment results show that when / is large relative to N^, and when >0, the FTSIA gives 
more accurate results than the statistical mean. The experiment results also indicate that A^^ should satisfy equation 
(6) in order for FTSIA to give more accurate results. 

The proposed FTSIA is used for integrating 1-dimensionaI data from the sensors that ideally should provide the 
same information. In the flow environment, however, the variables of interest can be a vector with more than one 
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dimensions. In such case, multi-dimensional fault-tolerant sensor integration algorithm needs to be developed. 
Moreover, it is very common that in the flow sensor network, each sensor provides partial and overlapping 
(complementary), or non-overlapping (cooperating) information [21] about the flow. Therefore, fault-tolerant 
complementary integration and cooperating integration are two other issues that will be also addressed in our future 
research. 
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AH references cited herein are incorporated in their entirety as if each were incorporated 
separately. This invention has been described with reference to illustrative embodiments and 
is not meant to be construed in a limiting sense. Various modifications of the illustrative 
embodiments* as well as additional embodiments of the inventi(xi» will be apparent to persons 
skilled in die ait upon reference to this description. 



